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ABSTRACT 


Mental  workload  has  been  very  difficult  to  describe  quantitatively.  An  exces¬ 
sive  workload  can  lead  to  a  decrease  in  accuracy  &  performance,  while  a 
sustained  high  level  of  workload  can  lead  to  mental  exhaustion.  Previous 
research  has  indicated  that  heart  rate  variability  and  evoked  potentials  in 
the  EEG  (electroencephalogram)  may  be  linked  to  mental  workload.  Unfor¬ 
tunately  most  of  the  work  to  date  has  examined  these  two  biocybernetic 
variables  independently  rather  than  jointly.  Recent  advances  now  allow  one 
to  measure  the  magnetic  fields  produced  by  the  brain  (MEG)  using  a  bQUID 
magnetometer  (Superconducting  Quantum  Interference  Device).  Much  of  the 
MEG  research  to  date  has  concentrated  on  lower  order  brain  processes 
rather  than  the  higher  cognitive  processes  associated  with  workload. 

The  current  research  examines  all  these  biocybernetic  variables  jointly 
in  an  effort  to  quantify  mental  workload.  A  paradigm  was  developed  to  vary 
several  aspects  of  mental  workload  and  verify  the  "hybrid  capacity"  model 
of  human  information  processing  that  was  developed  at  Purdue.  This 
research  included  constructing  a  data  acquisition  system  to  implement  this 
paradigm  and  simultaneously  record  heart  rate,  respiration,  EEG  &  MEG 
data.  It  is  believed  that  this  is  the  first  time  that  such  a  varied  data  set  was 
recorded  simultaneously. 
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Behavioral  measures  of  response  did  not  match  the  loading  levels 
predicted  by  the  model  for  all  conditions.  It  did  match  though  for  the 
higher  levels  of  loading.  Power  spectral  estimates  of  the  prestimulus  interval 
and  simple  averaging  did  not  appear  to  give  useful  measures  of  the  loading. 
The  MEG  data  was  highly  corrupted  due  to  head  movement.  The  use  of 
latency  corrected  averaging  (LCA)  allowed  the  identification  of  peaks  (P300 
&  N200)  which  were  capable  of  distinguishing  the  loading  levels  that 
showed  significant  differences  in  the  behavioral  measures  of  performance. 
The  heart  rate  variability  was  able  to  distinguish  between  the  same  loading 
levels  that  the  behavioral  variables  and  EEG  peaks  could. 
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CHAPTER  1 


INTRODUCTION 


1.1.  What  is  Loading 

As  engineering  systems  become  more  complex,  the  man-machine  inter¬ 
face  may  become  the  limiting  factor  to  increased  performance.  This  is  espe¬ 
cially  evident  in  complex  real  time  interfaces  such  as  between  a  pilot  and  a 
modern  airplane.  Cognitive  psychologists  have  developed  many  models  of 
human  information  processing  (review  in  Kantowitz,  1985a,  1974;  Taylor, 
1976;  Townsend,  1974;  Broadbent,  1971;  Sternberg,  1969).  In  order  to  jus¬ 
tify  the  model,  they  must  design  experiments  with  human  subjects  geared  to 
test  hypotheses  predicted  by  that  model.  The  measurement  of  workload  is 
of  central  importance  to  these  applications. 

The  term  “load”  or  “workload”  has  many  meanings  associated  with  it 
and  it  is  necessary  for  us  to  first  clarify  our  definition  of  it  before  we  can 
attempt  to  measure  it.  Weiner  (1982)  describes  workload  in  the  context  of 
mechanical  stress  and  strain.  As  long  as  the  body’s  regulatory  systems  are 
able  to  compensate  for  the  load  (e.g.  increased  cardiac  output),  the  subject 
is  said  to  be  stressed.  When  the  body  is  overloaded  and  can  no  longer  fully 
compensate  (e.g.  build  up  of  lactic  acid),  the  subject  is  said  to  be  under 
strain.  This  distinction  becomes  harder  to  make  for  psychological  systems. 
Overload  is  generally  perceived  as  a  change  of  emotional  state  (e.g.  anxiety). 
A  decrease  in  performance  accuracy  has  been  used  in  the  context  of  over¬ 
load,  however  it  has  also  been  associated  with  a  “high  workload”. 

1.2.  Measures  of  Mental  Load 

While  cardiac  output  and  oxygen  uptake  may  provide  a  good  measure 
of  physical  workload  there  is  no  corresponding  measure  for  mental  workload. 
This  is  primarily  because  there  is  no  adequate  model  for  human  information 
processing.  It  is  difficult  to  come  up  with  a  measure  for  something,  when  we 


are  not  even  sure  as  to  what  we  wish  to  measure.  Three  different  approaches 
are  generally  used  to  detect  mental  load  (Kalsbeek,  1973;  Gopher  &  Braune, 
1984): 

1)  personal  experience. 

2)  objective  description  of  the  task  parameters 

3)  measures  of  response  (either  behavioral  or  physiological) 

This  is  a  starting  point  for  clarifying  the  concept  of  mental  workload,  but 
there  are  still  many  shortcomings. 

Personal  experience  is  probably  more  a  measure  of  perceived  cognitive 
difficulty  than  actual  loading.  Hart  &  Bortolussi  (1984)  used  a  pilot  opinion 
survey  to  create  a  database  of  flight  scenarios  with  predetermined  levels  of 
pilot  workload.  Pilots  were  asked  to  rate  their  workload  and  other  parame¬ 
ters  as  they  went  through  the  predetermined  flight  scenarios.  There  was  a 
high  degree  of  correlation  between  the  perceived  workload,  stress  and  effort 
ratings  as  one  would  expect  in  such  a  circular  approach  to  measuring  work¬ 
load.  What  is  of  interest  to  note  though,  is  that  the  performance  ratings 
were  not  correlated  with  the  other  measures.  Gopher  &  Braune  (1984)  report 
similar  findings  in  a  different  experiment  and  conclude  with  “It  was  hinted 
that  subjective  measures  may  be  more  related  to  the  formal  properties  of  the 
task  than  to  the  details  of  performance.>> 

Kalsbeek  (1973)  also  discusses  “whether  the  concept  of  mental  load  is 
indivisible.”  It  may  be  that  a  single  measure  of  loading  is  not  sufficient  to 
describe  the  total  cognitive  processing  involved  with  a  given  task.  Some 
argue  (Moray  et  al.,  1979,  p.105)  that  the  only  true  measure  of  workload  is 
the  subject’s  own  interpretation  of  his  loading,  but  this  has  many  limita¬ 
tions.  Subjects  do  not  have  difficulty  assigning  a  number  to  a  given  task,  but 
how  does  one  select  the  range  of  numbers  to  chose  from.  Is  the  range  to  be 
fixed  and  invariant  with  time?  A  student  may  rate  the  first  exam  as  a  10  on 
a  difficulty  scale  of  1  to  10.  When  the  second  and  more  difficult  exam  is 
given,  is  it  then  ranked  as  a  15  on  the  same  1  to  10  scale  or  do  we  rescale  all 
previous  interpretations  of  loading  so  that  we  do  not  exceed  our  allowable 
range?  How  do  we  compare  ratings  between  students  who  may  base  their 
evaluations  on  completely  different  aspects  of  the  same  task.  One  student 
may  base  his  rating  on  the  time  aspects  of  the  test,  while  another  may  base 
his  on  the  style  of  the  exam,  while  yet  another  may  base  his  on  the  actual 
questions  asked. 
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This  problem  of  multiple  workload  aspects  for  a  given  task  is  not 
unique  to  the  subjective  method  of  load  description.  The  load  can  also  be 
described  in  terms  of  the  task  parameters,  but  how  does  one  compare  the 
workload  of  responding  to  a  pair  of  four  choice  stimuli  given  250  msec  apart 
with  a  pair  of  two  choice  stimuli  given  62.5  msec  apart.  Krol  &  Opmeer 
(1971)  report  different  trends  in  cardiovascular  activity  for  cockpit  workload 
during  simulated  landings  and  the  mental  load  of  parachute  jumpers.  They 
differentiate  this  as  informational  and  emotional  load.  A  single  measure  of 
load  appears  inadequate  to  describe  the  entire  spectrum  of  cognitive  process¬ 
ing,  so  how  then  do  we  proceed  in  our  quest  to  develop  metrics  for  workload 
classification. 

One  possible  solution  is  to  develop  a  simplified  model  of  cognitive  pro¬ 
cessing  and  define  our  workload  metric  in  light  of  this  model.  Chapter  2 
reviews  some  of  the  earlier  information  processing  models  leading  to  the 
development  of  the  “Hybrid  Capacity  Model”  (Kantowitz  &  Knight,  1976) 
which  is  the  model  used  for  the  current  research. 

1.3.  Biocybernetic  Variables  and  Experimental  Procedure 

The  current  research  utilizes  the  heart  inter-beat-interval  (R  to  R  inter¬ 
vals  of  the  electrocardiagram),  the  respiratory  rate,  event  related  potentials 
(electroencephalograms)  and  evoked  fields  (magnetoencephalograms)  to  gain 
an  understanding  of  mental  workload.  Chapter  3  reviews  the  measurement 
of  these  biocybernetic  variables  and  their  application  to  workload  studies. 

In  Chapter  4  the  experimental  setup  and  paradigm  are  discussed.  A 
double  stimulation  task  which  also  varied  the  number  of  choices  was  used 
for  the  current  study.  This  allows  us  to  vary  the  workload  by  increasing  the 
rate  that  information  is  presented  as  well  as  varying  the  degree  of  uncer¬ 
tainty  in  the  stimulus.  Response  times  and  accuracy  are  recorded  through 
out  the  experiment  so  that  we  have  converging  measures  (Kantowitz,  1985b); 
i.e.  workload  is  described  by  both  task  description  and  behavioral  response. 

An  extremely  powerful,  yet  flexible  data  acquisition  system  was 
designed  and  constructed  to  carry  out  this  research.  The  system  is  discussed 
in  Chapter  4  and  in  Appendix  B.  The  experiment  was  performed  in  the 
summer  of  1986  at  the  Aerospace  Medical  Research  Laboratory  (AMRL)  at 
Wright  Patterson  Air  Force  Base  (WPAFB). 
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1.4.  Data  Analysis  and  Results 

The  data  analysis  is  divided  into  two  chapters  in  this  report.  Chapter 
5,  entitled  “Preliminary  Data  Analysis”  and  Chapter  6,  entitled  “Advanced 
Analysis”.  The  results  and  comparisons  to  previous  work  are  discussed  in 
Chapter  7. 

The  preliminary  analysis  comprised  of  prescreening  the  data  to  reject 
eyeblink  contamination,  amplifier  saturation,  ectopic  heart  beats,  etc. 
Power  spectra  were  computed  from  the  prestimulus  interval  for  the  event 
related  potentials  (ERP)  and  the  evoked  fields  (EF).  Average  responses  were 
computed.  Mean  heart  rates  and  variances  were  computed. 

The  advanced  analysis  included  using  latency  corrected  averaging 
(Aunon  &  McGillem,  1977,  1978)  on  the  ERP  to  separate  out  significant 
peaks  in  the  response.  Spectral  techniques  were  used  on  the  respiration  and 
heart  rate.  Linear  regressions  were  performed  with  respect  to  the  response 
time. 
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CHAPTER  2 

HUMAN  INFORMATION  PROCESSING 

2.1.  Introduction 

The  scope  of  human  information  processing  is  so  broad  that  no  single 
document  could  do  it  justice.  The  goal  of  this  research  is  to  develop  metrics 
for  workload  and  so  our  discussion  of  human  information  processing  will  be 
i  focused  in  this  direction.  Several  different  techniques  are  available  for  load¬ 

ing  the  human  processing  capabilities.  The  first  we  shall  refer  to  as  single 
stimulation;  as  the  name  implies,  a  single  stimulation  is  presented  to  a  sub¬ 
ject.  Loading  can  be  varied  by  varying  the  complexity  of  the  stimulus,  i.e. 
j  locating  a  target  element  from  a  large  versus  a  small  set.  Stimulus  presenta- 

|  tion  can  also  be  degraded  to  increase  the  loading.  This  technique  predom- 

i  inantly  varies  the  perceptual  and  coding  aspect  of  information  processing 

and  is  consequently  not  often  used  when  one  is  more  concerned  with  the 

I  total  loading  in  an  information  channel. 

1 

2.2.  Double  Stimulation 

The  second  method  of  loading  is  referred  to  as  double  stimulation  (Kan- 
|  towitz,  1974)  and  is  the  more  dominant  method  since  it  allows  greater  flexi¬ 

bility  in  focusing  the  loading.  Loading  is  accomplished  by  forcing  the  infor¬ 
mation  channel  to  do  two  things  at  once.  This  can  be  accomplished  by 
decreasing  the  interstimulus  interval  (ISI)  between  two  successive  stimuli  to 
the  point  where  the  the  channel  does  not  have  sufficient  time  to  return  to  its 
i  unloaded  state  between  stimulations.  The  workload  can  then  be  increased 

by  decreasing  the  ISI  and  thereby  increasing  the  amount  of  time  that  the 
processing  will  be  dealing  with  both  stimulations  as  opposed  to  a  single 
stimulation.  If  the  ISI  becomes  extremely  short  the  two  stimuli  are  then 
“grouped”  and  perceived  as  a  single  stimulation  with  two  elements.  This 
technique  of  loading  utilizes  a  serial  approach  to  double  stimulation.  A 
parallel  approach  can  also  be  used  by  requiring  the  channel  to  handle  dual 
tasks  (primary  and  secondary)  at  the  same  time.  Response  time  (RT)  and 
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accuracy  are  recorded  as  a  measure  of  the  channel’s  performance.  The 
theory  states  that  channel  performance  ■will  decrease  as  channel  loading 
approaches  and  exceeds  channel  capacity.  Limitations  of  these  measures 
with  regard  to  channel  capacity  will  be  discussed  later.  Note  that  depend¬ 
ing  on  the  task  selected  (e.g.  continuous  reading),  either  or  both  of  these 
measures  may  not  be  available. 

Cognitive  psychologists  construct  behavioral  experiments  to  study 
human  information  processing.  Inferences  are  drawn  from  the  results  of 
these  experiments  about  the  underlying  process  and  models  are  then  con¬ 
structed  to  fit  these  inferences.  Models  constructed  from  loading  due  to  suc¬ 
cessive  stimuli  will  be  discussed  first.  These  models  will  then  be  revaluated 
in  light  of  some  recent  results  from  dual  task  experiments. 

There  are  several  basic  paradigms  that  can  be  used  in  a  serial  double 
stimulation  (Kantowitz,  1974,  p.  86).  The  most  common  of  these  is  the 
stimulus  1  -  response  1 ,  stimulus  2  -  response  2  (Si  -  Rl,  S2  -  R2)  paradigm 
since  it  requires  a  response  to  each  stimuli.  It  is  well  known  that  there  is 
delay  in  responding  to  the  second  stimulus  and  that  this  delay  increases  as 
the  ISI  is  decreased.  This  delay  has  been  called  the  psychological  refractory 
period  (PRP)  since  it  was  first  believed  that  the  human  information  process¬ 
ing  channel  became  refractory  or  disabled  for  a  brief  period,  similar  to 
refractory  period  for  nerve  fibers  (Telford,  1931).  It  is  now  believed  that  this 
effect  has  little  in  common  with  its  name,  but  nonetheless,  the  effect  has 
retained  its  descriptive  name.  This  paradigm  is  also  referred  to  as  the  PRP 
paradigm.  Other  paradigms  Si,  S2  -  R2  and  Si  -  Rl,  S2  which  do  not 
require  responses  to  both  stimuli  have  also  been  used  to  verify  models. 

2.3.  Limited  -  Capacity  Channel  Model 

The  traditional  model  for  double  stimulation  is  the  limited  capacity 
channel  model  (Broadbent,  1958;  1971)  shown  in  Figure  2.3.1.  The  limited 
channel  model  was  able  to  explain  changes  in  reaction  time  due  to  PRP, 
stimulus  modality  and  different  stimuli  probabilities.  This  became  the 
benchmark  for  other  models.  PRP  could  be  explained  by  S2  being  held  in 
short  term  storage  while  the  channel  was  dealing  with  SI,  since  the  limited 
capacity  channel  was  capable  of  handling  only  one  task  at  a  time.  Broad- 
bent  (1971  p.  313)  now  believes  that  delays  in  category  selection  is  a  better 
explanation.  Stimuli  with  higher  probabilities  of  occurrence  would  have  fas¬ 
ter  reaction  times  (RT)  due  to  the  selective  filter.  Delays  for  ISI  longer  than 
the  response  time  for  stimulus  1  (RTl)  could  be  explained  since  the  channel 
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was  updating  memory  after  it  had  finished  its  response  to  Si.  Unfortunately 
many  earlier  researchers  focused  on  PRP,  and  did  not  consider  the  effects  of 
S2  -  R2  on  RTl  (response  time  l).  It  was  found  that  RTl  increased  slightly 
due  to  S2  -  R2  (Herman  &  Kantowitz,  1970)  and  the  limited  channel  capa¬ 
city  can  not  offer  any  explanation  for  this  effect.  Broad  bent’s  review  of  the 
limited  channel  capacity  model  in  1971  does  not  explain  this  issue.  The  lim¬ 
ited  capacity  channel  also  has  difficulties  explaining  some  effects  of  dual  task 
loading  which  will  be  discussed  later. 

2.4.  Stage  Models 

The  stage  approach  to  model  human  processing  (Sternberg,  1969;  Tay¬ 
lor,  1976;  Townsend,  1974)  has  become  a  popular  method.  Figure  2.4.1 
shows  a  simple  four  stage  serial  model  of  the  reaction  process  for  a  binary 
decision.  Models  can  be  classified  as  serial,  parallel  or  hybrid  depending  on 
how  the  stages  are  connected  together.  Kantowitz  (1985b)  gives  an  abbrevi¬ 
ated  definition  of  a  stage  as  a  single  transformation  of  information,  while 
Townsend  (1974)  gives  a  more  precise  mathematical  definition. 

Taylor  (1976)  breaks  down  stages  into  smaller  terms  called  elements. 
An  element  is  identified  with  the  smallest  unit  of  information  processed  in  a 
particular  stage  of  a  particular  model.  The  processing  within  each  stage  can 
then  be  evaluated  along  several  lines  including:  parallel  versus  serial; 
exhaustive  versus  self  terminating;  and  limited  versus  unlimited.  Parallel 
and  serial  refer  to  whether  a  stage  can  process  more  than  a  single  element 
at  a  time.  Figure  2.4.2  gives  an  example  of  the  time  considerations  for  a  sin¬ 
gle  four  element  stage.  Exhaustive  and  self  terminating  refers  to  the  point 
at  which  the  stage  will  terminate  its  processing.  An  exhaustive  stage  will 
always  use  the  same  amount  of  time  since  it  will  process  all  elements  while  a 
self  terminating  process  will  end  after  it  processes  the  “critical  element”.  A 
self  terminating  process  is  a  possible  explanation  to  reaction  time  increasing 
linearly  with  the  number  of  choice  alternatives.  It  will  take  longer  on  aver¬ 
age  to  locate  the  critical  element  if  there  are  more  elements  to  search 
through.  Unlimited  refers  to  whether  the  stage  can  work  with  many  ele¬ 
ments  at  the  same  speed  at  which  it  can  handle  a  single  element.  Note  that 
parallel  processing  does  not  necessarily  imply  unlimited  processing  since  the 
processing  time  for  each  element  may  increase  as  the  number  of  elements 
increase  and  the  resources  of  the  stage  are  stretched. 
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A  diagram  of  the  flow  of  information  within  the  nervous  system,  as  conceived  by 


Figure  2.3.1  Limited  Capacity  Channel  (Broadbent  1058) 
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l*rocc"ing  singes  in  binary  classification.  Abuse  the  broken  line  arc 
•flown  ihe  four  factors  examined  lit  low  the  line  is  shown  the  analysis  of  RT 
inferred  from  aJJiti  c  relations  between  factor  pairs  IA2.  IAJ,  2 &>,  2A4.  and 
TA4.  the  linear  effect  of  factor  2.  and  other  considerations  described  in  the  text 
The  quality  of  the  test  stimulus  influences  the  duration  of  an  encoding  stage  in 
which  a  stimulus  representation  is  formed  This  representation  is  then  used  in  a 
verial-comparoon  stage,  "hose  duration  depends  linearly  on  size  of  positive  set; 
in  each  of  its  substages  the  represen'anon  is  compared  to  a  memory  represen¬ 
tation  of  one  member  of  the  set.  In  the  third  stage  a  binary  decision  is  made 
that  d  pends  on  whether  a  match  has  occurred  during  the  serial-comparison  stage 
that  precedes  it.  its  mean  duration  is  greater  for  negative  than  for  positive  de¬ 
cisions  The  selection  of  a  response  based  on  the  decision  is  accomplished  in  the 
final  stage  whose  duration  depends  on  the  relative  frequency  with  which  a 
response  of  that  type  is  required. 

Figure  2.4.1  Stage  Models  (Sternberg  1868) 
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Figure  2.4.2  Serial  &  Parallel  Element  Processing  (Taylor  1076) 
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Figure  2.6.1  Variable  Allocation  Model  (Kahneman  1073) 


2.5.  Stage  Capacities 

Given  the  complexity  that  even  a  single  stage  can  contain,  most  infor¬ 
mation  processing  models  consist  of  a  serial  channel  of  stages  to  prevent  the 
model  from  getting  unwieldy.  The  concept  of  the  loading  of  an  information 
channel  is  probably  best  described  in  terms  of  the  capacity  of  the  channel 
and  is  developed  by  Broadbent  (1071)  in  his  limited-capacity  channel  model 
of  performance.  Loading  can  be  described  as  to  how  much  of  the  available 
capacity  is  used.  Overloading  takes  place  when  the  available  capacity  is 
exceeded.  Of  central  concern  in  an  information  channel  is  how  the  available 
capacity  is  distributed  throughout  the  channel.  Reaction  time  for  a  process 
is  often  used  as  an  indicator  of  channel  capacity,  though  this  measure  has 
many  shortcomings  (Broadbent,  1971;  Townsend,  1974).  Kantowitz  (1985a) 
stresses  caution  when  using  RT  and  states:  “Ideally,  lag  should  not  be  used 
as  an  index  of  capacity  without  strong  converging  operations.”  A  converging 
criterion  would  be  changes  in  a  behavioral  variable  which  correlated  with 
changes  in  physiological  variables.  If  we  accept  for  the  moment  that  reac¬ 
tion  time  is  a  measure  of  capacity,  then  capacity  allocation  can  be  deduced 
by  decomposing  the  reaction  time  for  individual  stages  (Taylor,  1976). 

The  simpliest  method  of  decomposing  reaction  time  is  by  using  two 
tasks  that  differ  only  by  a  postulated  stage.  The  difference  in  reaction  times 
then  represents  the  time  for  that  stage.  This  method  is  known  as  the 
method  of  subtraction,  and  has  fallen  into  disfavor  though  it  has  never  been 
invalidated  by  modern  experimental  methods  (Taylor,  1976).  Sternberg 
(1969)  proposes  a  method,  the  method  of  additive  factors,  which  uses  a  fac¬ 
torial  experiment  to  decompose  reaction  time.  This  method  is  very  popular 
since  it  allows  one  to  determine  what  factors  affect  what  stages.  The 
rationale  of  the  method  is  that  factors  which  interact  in  a  factorial  experi¬ 
ment  must  influence  the  same  stage,  while  those  which  do  not  interact  prob¬ 
ably  influence  different  stages.  This  method  was  used  to  justify  the  binary 
response  model  shown  in  Figure  2.4.1.  Unfortunately  the  very  power  of  this 
method,  the  ability  to  assign  different  factors  to  different  stages,  will  lead  to 
its  downfall  when  we  discuss  the  dual  task  results.  The  basic  limitation  of 
these  methods  is  that  they  assume  each  stage  is  independent  and  has  its 
own  invariant  source  of  capacity.  Taylor  (1976)  expands  upon  Sternberg’s 
method  to  avoid  some  of  these  assumptions,  but  at  the  cost  of  adding  com¬ 
plexity  to  the  stage  description. 
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Kantowitz  (1985a)  advocates  that  the  concept  of  capacity  is  the  central 
topic  in  information  processing.  “Capacity  is  the  more  important  concept, 
since  while  capacity  can  exist  withovt  precise  definition  of  a  channel,  a  chan¬ 
nel,  at  least  in  psychological  research,  can  be  inferred  only  from  a  measure¬ 
ment  of  capacity.  ”  This  heuristic  solution  of  replacing  an  information  chan¬ 
nel  with  the  notion  of  capacity  sidesteps  all  the  difficulties  inherent  in  trying 
to  define  a  channel.  The  concept  of  a  channel  will  be  retained  in  this  discus¬ 
sion;  however  the  term  will  not  necessarily  imply  a  formal  description  of  the 
channel.  If  we  are  only  concerned  with  the  measurement  of  loading  and  we 
accept  that  capacity  utilization  is  an  acceptable  measure  of  loading,  then  it 
is  not  necessary  to  further  refine  our  model. 

2.6.  Variable  Capacity  Model 

Kahneman  (1973)  also  agrees  that  capacity  is  the  central  issue  and  has 
proposed  a  variable  allocation  model.  The  strength  of  this  model  is  that  it 
is  based  on  a  few  assumptions  and  does  not  attempt  to  provide  a  structure 
for  the  processing.  It  only  provides  an  abstract  model  to  justify  the  effects 
of  loading.  It  is  very  hard  to  disprove  something  that  is  very  unspecific  and 
makes  few  predictions.  The  basic  mechanism  of  the  variable  capacity  model 
is  that  total  capacity  increases  as  a  function  of  demand.  This  increase  is  at 
a  slower  rate  than  the  increase  in  capacity  flowing  to  the  primary  task. 
This  effect  is  illustrated  in  Figure  2.6.1.  The  channel  is  initially  able  to  sup¬ 
ply  the  primary  (or  single  task)  with  all  the  capacity  that  it  demands  so 
that  there  is  no  decrease  in  channel  performance  for  low  loading.  As  the 
loading  increases  the  channel  can  no  longer  keep  up  with  demand  and  per¬ 
formance  decreases.  For  dual  task  stimulations  the  amount  of  spare  capacity 
for  the  secondary  task  is  reduced  as  primary  task  difficulty  increases,  so  that 
both  task  performances  should  decrease.  Capacity  is  not  intended  to  be 
unlimited  in  this  model.  The  undefined  capacity  limit  though  is  only 
reached  in  the  asymptotic  sense.  One  of  the  shortcomings  of  this  model  is 
that  it  does  not  allow  any  way  of  distinguishing  between  task  demands  and 
the  performer’s  intentions  (Kantowitz,  1974,  p.94);  i.e.  how  does  the  model 
reflect  the  relative  importance  of  each  task  to  the  performer? 

2.7.  Dual  Task  Stimulation 

Dual  task  stimulation  (timesharing  paradigm)  involves  the  performance 
of  two  different  task  simultaneously.  Limited  capacity  channels,  stage 
models  and  variable  capacity  models  all  predict  that  performance  should 
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decrease  when  the  second  task  is  added  to  the  paradigm.  This  relationship 
is  presented  in  Figure  2.7.1.  If  the  addition  of  the  secondary  task  causes 
overloading,  as  with  the  addition  of  the  secondary  task  to  the  difficult  pri¬ 
mary  task,  then  primary  task  performance  will  decrease  considerably.  If  the 
channel  does  not  become  overloaded,  as  with  the  addition  of  the  secondary 
task  to  the  easy  primary  task,  then  there  will  only  be  a  slight  decrease  in 
performance.  If  the  two  lines  in  Figure  2.7.1  are  parallel,  it  implies  that  the 
same  decrease  in  performance,  due  to  the  introduction  of  the  secondary 
task,  was  simply  added  to  both  curves.  This  effect  is  called  the  additivity 
effect  (Kantowitz,  1976).  Simple  additivity  can  be  justified  by  all  of  the 
models  discussed  by  either  assuming  that  overload  was  not  reached  or  that 
overload  caused  a  decrease  in  secondary  task  performance. 

In  order  to  fully  understand  timesharing  paradigms,  it  is  necessary  to 
perform  each  component  task  at  two  levels  and  each  component  task  alone 
to  generate  a  single  stimulation  baseline  (Kantowitz,  1974).  These  results 
can  then  be  presented  as  shown  in  Figure  2.7.2.  Note  that  while  this  graph 
has  the  same  shape  as  Figure  2.7.1  it  represents  dual  task  performance, 
while  Figure  2.7.2  represents  the  interaction  of  single  task  and  dual  task 
performance.  Limited  channel  capacity  and  variable  capacity  models  only 
allow  for  the  possibility  of  additivity  as  shown  in  Figure  2.7.2  if  one  assumes 
that  two  simultaneous  difficult  tasks  do  not  strain  capacity.  Stage  analysis 
(Sternberg,  1969)  does  allow  for  additivity  at  this  level  by  then  placing  each 
task  in  a  different  parallel  stage  with  its  own  independent  source  of  capa¬ 
city. 

The  results  of  an  experiment  performed  by  Kantowitz  &  Knight  (1974, 
1976)  demonstrates  the  use  of  these  two  graphs.  They  obtained  additivity 
for  dual  task  interactions  (as  in  Figure  2.7.2),  but  overloading  for  single  task 
versus  dual  task  (as  in  Figure  2.7.1).  If  one  takes  liberties  (i.e.  considering 
no  dual  task  to  be  a  very  easy  task)  with  these  graphs,  the  results  of  Kan¬ 
towitz  &.  Knight’s  experiment  can  be  summarized  in  a  single  graph  as  in  Fig¬ 
ure  2.7.3.  The  inconsistencies  with  the  different  models  are  apparent  in  this 
graph.  The  first  part  of  this  graph  implies  that  capacity  is  strained  in  going 
from  the  single  task  to  the  dual  task,  but  that  capacity  is  no  longer  strained 
as  we  increase  the  difficulty  of  the  dual  task.  The  limited  capacity  channel 
and  the  variable  capacity  channel  are  not  able  to  explain  this  result.  The 
stage  approach  also  fails  since  the  first  part  of  the  graph  predicts  (via  the 
method  of  additive  factors)  that  the  primary  and  secondary  task  are  located 
in  the  same  stage,  while  the  second  part  of  the  graph  predicts  that  they  are 
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Figure  2.7.1  Primary  &  Secondary  Task  Interactions  (Kantowiti&Knight  1976) 
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Figure  2.7.2  Dual  Task  Interactions  on  Performance  (Kantowitz&Knight  1976) 
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located  in  different  stages.  While  the  existing  models  could  be  modified  to 
incorporate  these  results,  an  alternate  approach  would  be  to  construct  a  new 
model  based  on  these  findings. 


2.8.  Hybrid  Capacity  Model 

A  new  model  (Figure  2.8.1),  is  introduced  by  Kantowitz  &  Knight  (1976) 
to  explain  these  findings.  Although  the  authors  refer  to  it  simply  as  a  hybrid 
model,  we  shall  refer  to  it  as  the  hybrid  capacity  model.  This  model  is  basi¬ 
cally  a  conglomeration  of  the  other  models.  The  static  capacity  allocator 
(SCA)  can  be  used  to  divide  capacity  in  accordance  with  the  performers 
intentions  (e.g.  relative  task  importance).  The  entire  channel  has  a  single 
limited  source  of  capacity  which  is  dynamically  allocated  to  some  stages 
directly  and  to  some  stages  through  the  SCA.  This  implies  that  in  stages  1 
and  2  of  Figure  2.8.1  the  ratio  of  capacities  between  these  stages  is  set  at 
the  start  of  the  experiment,  but  that  the  actual  levels  of  capacity  will 
change  dynamically  with  task  demands. 

The  model  behaves  as  a  serial  stage  for  low  workloads  where  the  total 
capacity  required  does  not  exceed  the  available  supply.  It  is  not  necessary 
for  the  limited  capacity  source  to  dynamically  trade  off  capacity,  so  that 
each  stage  appears  to  have  its  own  independent  supply  of  capacity.  As 
workload  is  increased,  capacity  available  to  each  stage  is  not  sullicient  to 
keep  up  with  demand  and  the  model  behaves  as  a  limited  capacity  channel. 

This  model  appears  to  place  the  bottleneck  of  the  processing  channel  at 
the  later  output  stage  since  this  is  where  several  parallel  stages  feed  a  single 
stage.  There  are  many  different  views  as  to  where  capacity  limitations 
cause  bottlenecks  (Kantowitz,  1974),  and  some  even  argue  that  there  are  no 
bottlenecks  (Kahneman,  1973).  This  arrangement  was  probably  a  personal 
bias  of  the  author  as  expressed  by  Kantowitz  (1985a  p.  160):  “Indeed,  the 
present  author  is  guilty  of  this  oversimplification  by  claiming  that  response 
processes  are  a  more  important  locus  of  limitation  than  are  stimultis 
processes.  While  this  may  be  true  (e.g.  I  still  believe  it),  it  ignores  implica¬ 
tions  of  capacity  limitations  from  elements  for  limitations  upon  stages.” 
Note  that  this  statement  was  not  made  in  reference  to  any  model  in  particu¬ 
lar,  but  rather  during  a  discussion  of  capacity  limitations  in  stages. 

This  model  is  applicable  to  both  PRP  and  timesharing  paradigms.  The 
parallel  stages  controlled  by  the  SCA  are  not  necessarily  a  primary  and 
secondary  task  as  in  timesharing  paradigm.  They  can  also  be  used  to 
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Figure  2.7.3  Combined  Dual  Task  Effect  on  Performance 


Schematic  tcndcrinp  of  a  h> brid  model.  Solid  lines  represent  capacity,  allocation  to 
•tapes  while  dotted  lines  represent  information  floss  between  stapes  The  information  flu's 
inputs  to  stapes  I  and  2  have  been  deleted  to  iniprose  the  clarity,  ol  the  diapram  A  limned 
Capacity  source  dynamically  feeds  both  stupe  S3  and  a  static  capacity  allocator  (SCA)  which 
Partitions  capacity  between  stapes  SI  and  S2.  Stapes  SI  and  S2  operate  in  parallel  and  stupe  S3 
operates  in  serial  with  both  preeedinp  stapes.  Ha  eh  stupe  can  potentially  be  broken  down  into 
•tnaller  stapes  depcndtnp  upon  the  level  of  analy  sis.  For  example,  slape  S3  is  a  molar  representa¬ 
tion  of  response  selection,  execution  and  control  processes. 

Figure  2.8.1  Hybrid  Capacity  Mode!  (Kantowitz  &  Knight  1976) 


represent  attention  and  arousal  effects.  It  is  well  known  that  accuracy  can 
be  traded  off  for  speed  (Pachella,  1974).  The  SCA  could  probably  be  used 
with  some  postulated  states  to  explain  this  effect.  It  may  also  be  possible  to 
explain  effects  of  stimulus  quality  on  RT  by  using  the  SCA. 

While  this  model  does  not  give  any  insight  as  to  how  information  is 
actually  processed  in  a  channel,  we  must  remember  that  the  focus  of  this 
research  is  on  workload  (i.e.  capacity)  measurement  and  this  model  will  be 
sufficient  for  our  needs.  We  shall  simply  refer  to  the  stages  in  this  model  as 
early  and  late  stages  of  processing  without  worrying  specifically  about  what 
processing  is  performed  in  them. 


CHAPTER  3 


BIOCYBERNETIC  VARIABLES 


3.1.  Heart  Rate 

3.1.1.  Introduction 

Heart  rate  (HR)  is  normally  variable  over  time  in  a  relaxed  individual. 
HR  is  easily  determined  by  computing  the  R-R  interval  from  an  electrocardi- 
agram  (ECG).  A  lead  II  configuration  (Guyton,  1981)  is  used  in  order  to 
enhance  the  R  wave  for  easy  computer  identification.  One  electrode  is 
placed  on  the  upper  right  chest  and  the  other  on  the  left  ankle  with  the 
right  ankle  used  as  a  ground.  The  signal  is  amplified  and  then  threshold 
detected  by  the  computer  to  determine  the  exact  time  of  occurence  for  each 
R  wave. 

Numerous  regulatory  mechanisms  are  constantly  competing  amongst 
themselves  to  maintain  the  circulatory  system  in  what  can  best  be  described 
as  a  state  of  dynamic  equilibrium.  Heart  rate,  blood  pressure  and  cardiac 
output  are  intimately  related  and  must  be  able  to  respond  to  varying 
demands  that  the  body  may  place  on  them.  The  effect  of  these  regulatory 
mechanisms  on  the  heart  rate  has  been  an  area  of  significant  research  in  the 
study  of  mental  workload.  The  heart  rate  is  slowed  through  parasym¬ 
pathetic  (vagal)  stimulation  of  the  heart  which  primarily  slows  the  cardiac 
pacemaker  (sino-atrial  node).  Sympathetic  mechanisms  are  more  far  reach¬ 
ing  and  affect  the  systemic  system  as  well  as  the  heart.  Note  that  since  the 
heart  is  regulated  by  both  sympathetic  and  parasympathetic  mechanisms,  it 
is  possible  for  identical  HRs  to  arise  from  a  variety  of  circumstances.  The 
specifics  of  these  mechanisms  can  be  found  in  any  good  physiology  book  such 
as  Guyton  (1981). 


3.1.2.  Heart  Rate  Variability 

There  has  been  some  confusion  between  psychologist  and  physiologist 
over  the  term  sinus  arrhythmia  (SA).  For  the  purposes  of  this  report,  the 
term  SA  shall  be  used  to  refer  to  the  variation  of  heart  rate  due  to  respira¬ 
tory  effects  and  the  term  heart  rate  variability  (HRV)  to  refer  to  any  varia¬ 
tion  from  a  constant  heart  rate.  The  term  HRV  shall  not  refer  to  any 
specific  method  of  numerically  computing  the  variation  in  heart  rate  since 
numerous  methods  (Firth,  1973)  have  been  used  to  emphasize  different 
aspects  of  the  variability.  Interest  in  using  HRV  as  a  measure  of  mental 
load  arose  after  Kalsbeek  8c  Ettema  (1963)  reported  that  HRV  was  “gradu¬ 
ally  suppressed  when  increasing  the  difficulty  of  the  task.”  Since  that  time 
there  have  been  numerous  articles  (Kalsbeek,  1973;  Porges  et  al.,  1980; 
Akselrod  et  al.,  1981;  Sharit  8c  Salvendy,  1982)  and  a  symposium  leading  to 
an  entire  issue  of  Ergonomics  (Rolfe,  1973)  devoted  to  the  issue  of  HR  varia¬ 
bility.  Unfortunately  this  proliferation  of  research  has  led  to  much  confu¬ 
sion,  since  the  problem  is  not  as  straight  forward  as  it  might  have  appeared 
after  Kalsbeek’s  initial  findings. 

Many  researchers  simply  started  measuring  HR  under  real  life  working 
conditions  and  tried  to  draw  premature  conclusions  about  workload.  While 
changes  were  found  in  mean  heart  rate  (MHR),  instantaneous  heart  rate 
and  variance,  there  are  numerous  psychological  factors  which  could  account 
for  these  changes  under  normal  work  environments  (Firth,  1973).  Some 
researchers  (Sharit  8c  Salvendy,  1982;  Heslegrave  et  al.,  1979)  advocate  the 
use  of  the  mean  square  successive  differences  (MSSD)  as  a  measure  of  HRV 
and  the  MHR  instead  of  the  sample  variance  of  the  interbeat  interval. 
Mulder  (1973)  used  numerous  measures  of  HR  variability:  MHR,  second, 
third  8c  fourth  moments  of  R-R  interval,  sum  of  absolute  differences  of  R-R 
interval,  number  of  reversal  points,  etc.  in  an  attempt  to  measure  workload. 
It  was  found  that  there  were  many  side  effects  of  the  task  (e.g.,  attention 
motility,  motor  load  from  button  pushing,  respiration,  stimulus  frequency) 
which  contributed  to  overall  changes  in  HR  variability.  This  is  one  of  the 
major  stumbling  blocks  to  understanding  how  mental  load  affects  HR. 

Mulder  (1973)  concludes  his  work  with  “In  future  research  more  atten¬ 
tion  must  be  paid  to  the  types  of  tasks  and  the  behavioral  mechanism 
involved,  to  their  presentation  mode  and  to  a  more  sophisticated  analysis  of 
heart  rate  variability Sharit  8c  Salvendy  (1982)  used  two  contrasting  tasks, 
one  was  largely  based  on  visual  perception  (low  mental  load)  and  the  other 
required  the  mental  solution  (high  mental  load)  of  arithmetic  problems. 
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MHR  appeared  to  be  far  more  sensitive  than  HRV  in  distinguishing  between 
these  two  tasks.  This  is  in  contrast  to  Kalsbeek  &  Etterma  (1963)  who  con¬ 
tend  that  HRV  should  diminish  with  increased  load,  but  that  MHR  should 
remain  relatively  constant.  It  was  contended  that  the  differences  were  due 
to  attentional  loads  versus  informational  loads.  These  results  suffer  from  the 
lack  of  a  converging  criterion  (Kantowitz,  1985a,  1985b)  as  discussed  else¬ 
where  in  this  report. 

The  lack  of  a  precise  definition  of  mental  workload  still  clouds  much  of 
the  issue.  The  only  way  to  verify  that  HRV  depends  on  the  workload  is  to 
vary  the  workload,  but  that  is  exactly  what  we  do  not  know  how  to  meas¬ 
ure.  A  behavioral  measure  is  needed  to  provide  the  converging  criterion. 
The  level  of  vigilance  and  arousal  are  generally  considered  to  be  related  to 
mental  load.  It  has  also  been  shown  (Kalsbeek,  1963)  that  noise  and  open 
eyes  tended  to  suppress  some  of  the  HRV  found  with  the  eyes  closed.  It  may 
be  that  HRV  is  responding  more  to  arousal  levels  than  workload  levels. 
While  the  autonomic  controls  which  regulate  the  cardiovascular  system  are 
well  understood,  it  is  not  clear  how  workload  affects  these  regulatory 
mechanisms. 

3.1.3.  Spectral  Techniques 

The  most  promising  area  of  research  in  HR  variability  appears  to  lie 
with  spectral  processing  techniques  (Sayers,  1973,  1975;  Mulder  et  al.,  1973; 
Porges  et  al.,  1980;  Akselrod  et  al.,  1981;  DeBoer  et  al.,  1984;  Sharit  &  Sal- 
vendy,  1982;  Kantowitz,  1985b).  The  R-R  interval  is  the  variable  that  is  usu¬ 
ally  recorded  in  the  course  of  an  experiment.  The  data  can  be  interpolated 
to  form  an  equally  time  spaced  data  set  and  the  power  spectrum  estimate 
(PSE)  of  the  interbeat  interval  is  then  computed  (Kantowitz,  1985b;  Mulder 
et  al.,  1973).  Others  (DeBoer,  1984;  Sayers,  1975)  compute  the  PSE  directly 
from  the  R-R  intervals,  by  considering  the  series  to  be  equally  beat  sampled 
(rather  than  equally  time  sampled),  or  by  direct  transform  techniques,  which 
do  not  rely  on  samples  being  equally  spaced.  DeBoer  (1984)  has  shown  that 
these  two  direct  methods  are  intimately  related  and  should  lead  to 
equivalent  results  when  all  effects  are  accounted  for.  The  suppression  of 
HRV  activity  shows  up  as  a  reduced  PSE  with  the  remaining  energy  shifting 
closer  to  the  lower  frequencies.  Spectral  techniques  allow  us  to  use  several 
measures  of  HRV  (Kantowitz,  1985b)  from  different  bands  of  the  PSE  which 
are  associated  with  different  regulatory  mechanisms. 
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Figure  3.1.1  shows  a  PSE  of  the  interbeat  interval  computed  from 
workers  during  a  period  of  low  activity  and  from  a  period  of  industrial  work¬ 
load  of  an  unreported  nature.  The  components  in  Figure  3.1.1  related  to 
respiratory  activity  (~  0.35  Hz),  vasomoter  (blood  pressure)  activity  (~  0.1 
Hz)  and  thermal  activity  (~  0.025)  were  determined  by  independently  corre¬ 
lating  these  variables  previous  to  collecting  this  data.  The  respiratory 
activity  affects  the  blood  pressure  as  well  as  the  heart  rate.  The  activity  at 
0.1  Hz  related  to  vasomoter  activity  refers  to  the  internal  oscillatory  nature 
of  the  blood  pressure  regulatory  system.  The  blood  pressure  regulatory  sys¬ 
tem  is  complex  and  in  part  depends  on  the  property  of  smooth  muscle  tissue 
and  neural  control  from  the  brain  stem.  Due  to  the  nature  of  the  regulatory 
mechanism,  it  is  possible  for  an  external  periodic  signal  to  entrain  the  oscil¬ 
lation  (replace  the  normal  oscillation  with  the  external  driving  frequency). 
This  effect  is  similar  to  that  of  the  hierarchy  of  cardiac  pacing  or  that  of  a 
synchronizable  electronic  oscillator.  Respiration  has  been  known  to  entrain 
the  blood  pressure  regulation  under  certain  conditions  (Sayers,  1973).  Since 
the  brain  stem  is  an  integral  part  of  the  blood  pressure  regulatory  mechan¬ 
ism,  it  may  be  that  mental  loading  may  somehow  influence  this  regulation. 
Figure  3.1.1  shows  that  vasomoter  activity  is  greatly  reduced  with  the 
increased  loading,  but  one  should  not  jump  to  premature  conclusions. 

Spectral  techniques  not  only  allow  us  to  identify  some  of  the  external 
influences  on  HR  irregularity,  they  also  allow  us  to  correct  for  these  effects 
(Bendat  &  Piersol,  1971  chap.  5;  Porges  et  al.,  1981).  Caution  should  be 
exercised  before  we  blindly  correct  for  all  these  effects.  The  mechanisms  of 
mental  workload  are  not  clearly  understood  and  the  primary  effect  of  load¬ 
ing  may  indeed  be  on  the  regulatory  mechanism.  The  regulatory  mechanisms 
described  above  account  for  82%  (Sayers,  1973)  of  the  variance  in  the  HR. 
Clearly  there  are  some  effects  that  should  be  removed.  Mulder  et  al.  (1973) 
reports  a  peak  in  the  power  spectrum  at  the  stimulus  frequency.  This  effect 
should  be  removed  since  it  is  a  function  of  the  task  and  not  the  loading. 
The  effects  of  respiration  are  often  removed  before  computing  the  variance 
(Porges  et  al.,  1981;  Kantowitz,  1985).  Any  effect  of  workload  on  respiration 
and  its  effect  on  autonomic  control  is  then  lost.  If  one  deals  with  the  PSE 
rather  than  merely  the  variance  it  is  not  necessary  to  categorically  remove 
the  effects  of  the  autonomic  control,  since  they  can  be  examined  on  an  indi¬ 
vidual  basis. 

Clearly  there  are  workload  effects  on  the  HR  variability,  but  more  study 
needs  to  be  done  to  isolate  these  effects.  Experiments  should  be  designed  to 
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Figure  3.1.1  Components  of  Heart  Rate  Variability  (Sayers  1973) 


minimize  the  effects  of  parameters  unrelated  to  workload.  The  use  of  con¬ 
verging  operations  will  also  lend  more  weight  to  any  results  obtained.  It 
may  prove  that  HR  variability  has  limited  usefulness  in  real  world  workload 
situations  until  we  have  a  better  understanding  of  the  underlying  effects. 

3.2.  Respiration 

The  respiratory  rate  is  usually  measured  in  workload  studies  so  that  the 
effect  of  sinus  arrhythmia  (SA)  on  the  HR  can  be  removed  or  accounted  for. 
During  normal  respiration,  the  heart  rate  may  swing  by  5%  due  to  the 
respiratory  cycle,  but  this  may  increase  to  30%  during  deep  breathing  as 
shown  in  Figure  3.2.1  (Guyton,  1981).  Respiratory  rate  is  normally  quite 
variable  so  that  its  use  as  a  measure  of  mental  workload  may  lie  more  in  its 
relation  to  the  HRV  than  as  an  independent  measure  by  itself. 

Since  the  central  nervous  system  (CNS)  is  an  integral  part  of  the  auto¬ 
nomic  regulation,  loading  of  the  CNS  may  have  some  effect  on  the  perfor¬ 
mance  of  the  autonomic  regulation.  Inspiration  causes  an  occlusion  of  the 
veins  which  leads  to  a  temporary  decrease  in  blood  pressure.  The  carotid 
sinus  reflex  then  causes  the  heart  rate  to  increase  in  response  to  the  occlu¬ 
sion  of  the  veins.  Inspiration  also  stimulates  stretch  sensors  that  inhibit 
vagal  effects  (Porges  et  al.,  1981)  which  also  causes  the  heart  rate  to 
increase.  Respiration  is  also  inhibited  during  inspiration  by  these  stretch 
receptors  due  to  the  Hering-Breuer  reflex.  Valentinuzzi  and  Geddes  (1974) 
present  a  model  for  the  components  of  SA  (Figure  3.2.2).  They  claim  that 
the  dominant  mechanism  of  SA  during  normal  conditions  is  the  central  com¬ 
ponent  due  to  a  coupling  (K)  between  the  respiratory  center  and  cardiac 
center  in  the  medulla.  This  conclusion  is  based  on  sustained  heart  rate  vari- 
abilty  during  voluntary  apnea  (holding  breath)  in  humans  and  drug  induced 
apnea  in  dogs.  Another  possibility  is  that  of  independent  variations  of  fre¬ 
quency  within  the  cardiac  center  itself.  If  these  mechanisms  are  affected  by 
workload,  then  the  loading  may  show  up  as  a  change  in  the  heart  rate  itself 
or  in  the  relationship  between  heart  rate  and  respiration. 

3.3.  Event  Related  Potentials 
3.3.1.  Introduction 

There  has  been  such  a  cornucopia  of  research  in  the  area  of  signal  pro¬ 
cessing  of  electroencephlagram  (EEG)  and  event  related  potentials  (ERP) 
that  it  is  difficult  to  establish  a  starting  point  for  this  discussion.  The  most 
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arterial  pressure  anJ  blood  rcluin  in  the  heart. 


Figure  3.2.2  Model  for  Components  of  Sinus  Arrhythmia 
(Valentinussi  &  Geddes  1974) 
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recent  and  broadest  review  on  the  subject  of  “electromagnetic  signal  of  the 
human  brain”  was  done  by  Gevins  (1984),  though  a  more  specific  three  part 
review  on  “signal  processing  in  evoked  potential  research”  was  done  in  1981 
(Aunon  et  al. ;  Childers  et  al.;  McGiliem  et  al.).  In  order  to  limit  this  section 
to  a  reasonable  size,  we  will  limit  our  discussion  to  methods  that  are 
relevant  to  the  current  research  only.  This  will  primarily  include  discussion 
of  the  ERF  with  respect  to: 

1)  advances  in  measurement  technique 

2)  classification 

3)  cognitive  uses 


3.3.2.  Measurement  Techniques  for  ERP 

It  was  in  1875  that  Dr.  Richard  Canton  first  recorded  electrical  activity 
from  the  exposed  surface  of  rabbits’  and  monkeys’  brains  (Braizer,  1984). 
Fifty  four  years  later,  Hans  Berger  recorded  the  first  human  EEG  and  hoped 
that  it  would  provide  a  “window  on  the  mind”  (Gevins,  1984).  One  might 
argue  that  since  that  time  we  have  found  the  window,  but  that  it  is  so  dirty 
that  we  can’t  quite  make  out  what  it  is  that  we  see. 

ERPs  measured  from  the  scalp  typically  vary  from  a  few  tenths  of  a 
microvolt  to  several  microvolts  and  are  imbedded  in  the  ongoing  EEG 
waveform  whose  amplitude  is  typically  10  to  30  microvolts  (McGiliem  et  al., 
1981).  The  ERP  is  usually  considered  to  be  a  deterministic  signal  added  to 
a  random  noise  process  (on  going  EEG)  with  signal-to-noise  ratios  (SNR)  as 
low  as  -20  dB.  Improvements  in  signal  estimation  have  progressed  in  two 
directions,  averaging  techniques  and  filtering  of  single  trials.  Most  of  these 
techniques  assume  that  the  EEG  is  an  independent  noise  process,  which  may 
not  be  entirely  valid. 


The  earliest  “averaging”  technique  was  simple  photographic  superposi¬ 
tion  (Dawson,  1947).  The  most  common  and  simplest  average  is  the  conven¬ 
tional  straight  average,  which  predicts  an  increase  in  the  SNR  by  the  root  of 
the  number  of  trials  (Aunon  et  al.,  1981a).  Woody  (1967)  introduced  the 
cross-correlation  average  which  shifted  each  trial  by  the  amount  that  would 
cause  the  cross  correlation  between  that  trial  and  the  conventional  average 
would  have  a  maximum  at  zero  lag.  The  shifted  trials  were  then  averaged. 
This  allowed  one  to  compensate  for  random  latency  shifts  between  trials. 
Aunon  &  McGiliem  (1977,  1978)  expanded  upon  this  idea  and  developed  the 
latency  corrected  average  (LCA)  which  allowed  one  to  compensate  for 


individual  latency  shifts  for  each  peak.  Smoothing  was  added  to  make  the 
resulting  waveform  continuous  again.  These  techniques  have  been  compared 
(Aunon  et  al.,  1981;  McGillem  et  at.,  1985)  and  it  is  not  suprising  that  the 
continuous  LCA  yielded  the  more  enchanced  representation  of  the  ERP. 

Wiener  filtering  (Carlton  &  Katz,  1980)  is  probably  the  most  common 
filtering  that  has  been  applied  to  ERP  (McGillem  et  al.,  1981).  Minimum 
mean  square  error  (MMSE)  and  maximum  signal  to  noise  ratio  (MSNR) 
techniques  (Aunon  &  McGillem,  1975;  McGillem  &  Aunon,  1977;  Aunon, 
1978;  McGillem  et  al.,  1981)  have  also  been  used.  Recently  a  more  powerful 
class  of  time  varying  (Yu  &  McGillem,  1983)  and  multichannel  time  varying 
filters  (McGillem  &  Aunon,  1985;  Westerkamp,  1985)  have  been  utilized. 

3.3.3.  Classification  of  ERP 

Perhaps  the  simplest  way  to  classify  ERPs  is  by  visual  inspection. 
Unfortunately,  while  this  method  may  have  limited  success  when  applied  to 
the  average  ERP,  it  is  virtually  impossible  when  applied  to  a  single  trial 
ERP.  The  basic  theory  in  classification  is  to  develop  a  feature  set  which  can 
be  used  to  distinguish  between  the  classes.  There  are  numerous  heuristic 
and  statistical  manners  in  which  these  features  can  be  selected  (Fukunaga, 
1972;  Gevins,  1980;  McGillem  et  al.,  1981).  These  features  can  be  as  simple 
as  a  present  or  not  present  decision  or  they  can  be  projections  on  basis  func¬ 
tions  in  an  expansion  of  the  signal  space.  If  we  view  the  signal  as  being  the 
dependent  variable,  we  can  map  the  signal  in  a  space  that  is  spanned  by  the 
feature  set  (Franks,  1981).  If  each  class  represents  a  random  process  or  is 
the  sum  of  a  deterministic  signal  (e.g.  ERP)  plus  a  random  process  (e.g. 
EEG),  then  we  will  get  a  mapping  similar  to  the  one  shown  in  Figure  3.3.1 
(Nagy,  1968).  A  hyperplane  can  then  be  used  to  separate  the  classes  pro¬ 
vided  certain  requirements  are  met  (Lay,  1982).  The  hyperplane  selected 
and  the  manner  in  which  the  features  are  selected  define  the  classification 
algorithm.  This  is  a  very  general  overview  of  how  classifiers  are  used  and 
some  of  the  more  popular  ones  are  discussed  subsequently. 

A  popular  classification  technique  is  Linear  Step-wise  Discriminant 
Analysis  (LSDA)  where  features  are  selected  by  determining  the  data  vectors 
for  which  mean  values  of  each  class  are  most  different  as  measured  by  a 
one-way  analysis  of  variance  F  statistic  (Donchin  &  Herning,  1975). 
Discriminant  functions  can  also  be  selected  by  using  the  correlation  function 
(correlation  classifiers)  or  the  maximum  likelihood  classifier  (Bayes  rule) 
(Sencaj  et  al.,  1979).  The  Bayes  linear  classifier  will  provide  optimal  results 


h' 

3 

7.» 


I 

& 


i 
* « 


•s 

1 


»• 

.V 


i;:5 


•il 


\i\ 

IV 

pi 

b] 


/«• 


f/  »• 

/hj  /  M? 


CLASS  C1 


CLASS  C* 


Common  types  ol  linear  categorders  *  and  O  indicate  the  training  samples  in 
classes  C'  and  C’  respectively  The  ellipses  are  the  equiprobabilily  contours  on  the  postulated 
distributions  in  the  lest  data  Thesubscrtpls  associated  with  the  hyperplanes  and  weight  vectors 
pertain  to  the  following  caiegonicrs  (I)  distance  to  means.  (2)  correlation.  (J)  approsimale 
masimum  likelihood,  (4)  Anderson-Bahadur,  (5)  discrimim.nl  analysis.  (6) 
discriminant  analysis.  (7)  trainable  machine,  (8)  optimal  quadratic  boundary,  and  (9)  L5UA 


Figure  3.3.1  Linear  Classifiers  (Nagy  1868) 
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if  all  the  distribution  statistics  are  normal  (Fukunaga,  1972  p.  90).  If  the 
covariance  matrices  for  the  two  classes  are  not  equal  then  a  linear  classifier 
will  not  provide  optimum  classification  and  a  two  class  case  will  require  a 
quadratic  classifier  for  optimum  performance  (McGillem  et  al.,  1981).  Aunon 
et  al.  (1982)  gives  a  comparison  of  the  performance  of  linear  versus  qua¬ 
dratic  classifiers  for  ERPs. 

One  way  of  selecting  features  is  by  using  the  forward  sequential  feature 
selection  (FSFS)  method.  In  this  procedure  the  best  single  feature  is  selected 
first  using  a  "training  set"  of  data.  The  next  feature  selected  is  the  one  that 
works  best  in  combination  with  the  first  one,  and  so  on.  This  will  not  yield 
the  best  performance,  but  the  computational  burden  of  exhaustively  check¬ 
ing  the  performance  of  all  possible  combinations  of  features  can  be  enormous 
(Moser,  1984;  Halliday  et  al.,  1985). 

Other  methods  of  selecting  features  include  Principal  Component 
Analysis  (PCA),  where  features  (basis  functions)  are  selected  in  such  a  way 
as  to  maximize  the  rate  of  reduction  of  the  residual  variance  (Donchin,  1966; 
Van  Rotterdam,  1970;  John  et  al.,  1973).  This  is  accomplished  by  reordering 
the  eigenvectors  produced  by  the  covariance  matrix  in  order  of  decreasing 
eigenvalues.  Critics  of  this  method  claim  that  the  resulting  components  may 
have  little  bearing  to  the  original  vectors  and  hence  a  rotation  of  the  com¬ 
ponents  is  recommended  (varimax  rotation)  to  force  the  axes  formed  by  the 
components  to  align  with  a  signal  vector  (Rosier  &  Manzey,  1981;  Wastell, 
1981a;  Wood  &  McCarthy,  1984).  Childers  et  al.  (1982)  utilizes  a  feature 
extraction  algorithm  where  the  eigenvectors  are  ordered  by  decreasing  the 
value  of  Fisher’s  ratio  instead  of  simply  decreasing  order. 

3.3.4.  Cognitive  Uses  of  ERP 

The  overriding  majority  of  ERP  studies  done  for  cognitive  purposes 
have  focused  on  the  effects  on  the  positive  peak  located  approximately  300 
msec  (P300  or  P3)  after  stimulation  (Pritchard,  1981;  Ruchkin  et  al.,  1975; 
Naatanen,  1974;  Begleiter  et  al.;  1983,  Donchin  et  al.,  1973;  McCarthy  & 
Donchin,  1981;  Israel  et  al.,  1980).  At  one  point  it  was  believed  that  P300 
was  a  measure  of  the  amount  of  activity  of  a  general  purpose  cortical  proces¬ 
sor  (Donchin  et  al.,  1973),  however  at  this  point  there  appears  to  be  a  gen¬ 
eral  consensus  that  P300  is  predominantely,  if  not  solely,  a  measure  of  “per¬ 
ceptual  load”.  Israel’s  (1980)  classic  workload  experiment,  where  P300  was 
found  to  vary  with  loading,  utilized  a  visual  monitoring  task  (large  percep¬ 
tual  load)  to  vary  loading.  McCarthy  &  Donchin  (1981)  found  that  P300 


latency  was  affected  only  by  stimulus  discriminability.  Other  papers  detail 
the  effect  of  stimulus  incentive  value  (Begleiter  et  al.,  1983)  and  stimulus 
probability  (Ruchkin  et  al.,  1975).  Pritchard’s  paper  (1981)  provides  an 
excellent  review  of  P300,  while  Naatanen  (1973)  provides  a  older  review  of 
ERP  and  attention. 

A  noticeable  exception  to  the  P300  fixation  has  been  Gevins.  In  one 
experiment  (Gevins  et  al.,  1979a),  a  series  of  complex  cognitive  tasks  was 
performed  to  exercise  many  different  aspects  of  the  information  channel. 
The  results  were  analyzed  in  a  variety  of  ways  and  it  was  found  that  there 
was  significant  change  in  many  of  the  frequency  bands  (predominately  theta) 
which  correlated  with  the  different  tasks.  Unfortunately,  in  a  later  experi¬ 
ment  (Gevins  et  al.,  1979b),  where  the  noncognitive  aspects  of  the  task  (e.g., 
limb  movement,  eye  movements,  etc.)  were  more  tightly  controlled,  the 
results  could  not  be  replicated.  This  illustrates  one  of  the  major  difficulties 
with  cognitive  paradigms:  ensuring  that  the  noncognitive  aspects  of  an 
experiment  do  not  corrupt  the  data.  The  paradigm  in  the  proposed  research 
is  much  simpler  than  the  task  performed  in  these  studies  and  should  be  less 
subject  to  this  effect. 

3.4.  Biomagnetism 

3.4.1.  Introduction 

The  study  of  biomagnetic  phenomena  is  still  in  its  infancy  when  com¬ 
pared  to  its  older  brother,  bioelectric  phenomena.  Figure  3.4.1  (Katila, 
1981)  shows  a  comparison  between  major  bioelectric  and  biomagnetic 
phenomena,  along  with  their  landmark  dates.  Insight  into  the  current  state 
of  research  in  biomagnetism,  as  of  1984  (last  Biomagnetism  Conference),  is 
possible  by  looking  at  some  statistics  regarding  the  level  of  research  being 
conducted.  Cohen  (1985)  reports  that  there  are  ~  50  groups  worldwide  con¬ 
ducting  research  in  biomagnetisin,  with  a  total  of  200  to  250  researchers 
including  graduate  students.  Only  50  to  60  refereed  papers  are  published 
annually,  excluding  those  presented  at  the  conference  held  every  two  years. 
Approximately  one  third  to  one  half  of  the  papers  are  brain  related.  This 
does  not  represent  a  very  large  worldwide  effort  for  a  field  that  has  so  much 
potential.  There  is  still  a  lot  to  be  done  in  order  for  biomagnetism  to  reach 
a  mature  state. 
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Various  bioelectric  phenomena  and  their  blomagnetlc 
counterparts.  Corresponding  bioelectric  signals  do  not  exist 
for  pure  magnetization  phenomena.  References  to  the  pio¬ 
neering  works  of  various  biomagnetic  measurements  are  made  in 
the  reference  list  from  no  1  to  no  12.  All  known  biomagnetic 
fields  are  not  included  in  the  Table. 

Figure  3.4.1  Bioelectric  and  Biomagnetic  Fields  (Katila  1981) 


The  first  recording  of  the  magnetic  fields  produced  by  the  brain,  magne¬ 
toencephalogram  (MEG),  was  done  in  1968  and  took  eight  minutes  of 
averaging  to  produce  only  one  cycle  of  the  alpha  rhythm  (Cohen,  1985). 
MEG  study  was  not  practical,  due  to  the  extremely  low  intensity  of  the 
fields,  as  illustrated  in  Figure  3.4.2  (Williamson  &  Kaufman,  1981),  until  the 
development  of  the  Super  Conducting  Quantum  Interference  Device 
(SQUED).  MEGs  are  on  the  order  of  1  picotesla  (pT),  while  the  brain’s  mag¬ 
netic  evoked  fields  (EF)  are  on  the  order  of  0.1  pT  or  100  femtoteslas  (fT). 
The  earth’s  magnetic  field  is  on  the  order  of  70  microteslas  (almost  eight 
orders  of  magnitude  greater  than  EFs),  with  a  spatial  gradient  of  ~  3 
nanoteslas.  Current  state  of  the  art  measuring  systems,  DC  SQUID  in  a 
shielded  chamber,  have  a  noise  level  as  low  as  15  fT  /  \/Hz,  as  measured 
while  collecting  data  at  Wright  Patterson  Air  Force  Base  (WPAFB).  Figure 
3.4.3  shows  the  DC  SQUID  and  support  system  used  at  WPAFB,  while  Fig¬ 
ure  3.4.4  shows  the  shielded  room.  The  channel  was  bandpass  filtered  from 
0.1  to  25  Hz.  This  implied  a  channel  noise  of  75  fT  (25  Hz  @  15  fT/VHz), 
which  is  not  much  smaller  than  the  EF. 

3.4.2.  SQUIDS 

A  SQUID  system  is  basically  a  very  low  noise,  very  high  gain  amplifier 
connected  to  a  sensory  loop.  The  schematic  diagram  for  a  DC  SQUID  sys¬ 
tem  is  shown  in  Figure  3.4.5  from  Williamson  &  Kaufman  (1981)  and  a  com¬ 
plete  mathematical  explanation  is  available  in  that  reference  or  Sarwinski 
(1977).  The  SQUID  itself  is  a  superconducting  loop  (although  the  term  is 
often  used  to  imply  the  entire  measurement  system)  that  contains  either  one 
(RF  SQUID)  or  two  weak  links  (DC  SQUID).  A  superconducting  weak  link 
is  a  Josephson  junction  which  will  oscillate  when  biased.  This  oscillation 
must  have  an  integral  number  of  wavelengths  on  the  loop  since  it  is  a 
closed  loop.  It  can  be  shown  that  this  requirement  translates  to  an  integral 
number  of  flux  quanta  within  the  loop. 

The  weak  link,  Josephson  junction,  requires  a  certain  amount  of  bias 
current  to  cause  this  oscillatory  effect.  The  DC  or  RF  notation  refers  to  the 
manner  in  which  this  bias  current  is  supplied.  DC  SQUIDS  use  a  dc  current 
source  as  shown  in  Figure  3.4.5.  RF  SQUIDS  use  an  RF  source  which  is 
inductively  coupled  and  hence  a  larger  superconducting  loop  is  needed.  This 
larger  loop  size  is  responsible  for  the  greater  noise  in  the  RF  SQUID.  DC 
SQUIDS  became  available  in  the  early  1980s  and  yielded  almost  an  order  of 
magnitude  better  performance  than  the  RF  SQUIDS  due  to  their  lower  noise 
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Figure  3. 4. 3  DC  SQUID  and  Mounting  Systems 


Figure  3.4.4  Shielded  Room 


Figure  3.4.5  Schematic  for  a  DC  SQUID  System  (Williamson  &  Kaufman  1981) 


characteristic.  Large  numbers  of  trials  were  necessary  with  RF  SQUIDS  in 
order  to  record  an  EF.  This  number  has  dropped  considerably  with  the 
introduction  of  the  DC  SQUID  (Lewis  et  al.,  1985). 

An  input  signal  is  inductively  coupled  with  the  superconducting  loop. 
As  the  input  signal  increases  and  attempts  to  add  flux  to  the  loop,  the 
superconducting  quantum  nature  of  the  loop  forces  an  opposing  current  to 
ensure  that  the  total  flux  in  the  loop  remains  at  the  quantum  level.  Unregu¬ 
lated,  this  effect  would  continue  until  the  input  signal  increased  enough  to 
cause  the  flux  in  the  loop  to  jump  by  one  quantum  level,  at  which  point  the 
opposing  current  would  suddenly  drop  to  zero.  The  system  is  regulated 
though  current  feedback  which  is  used  to  keep  the  loop  from  jumping  quan¬ 
tum  levels.  The  amount  of  feedback  current  can  then  be  used  as  a  measure 
of  the  input  signal.  Note  that  while  the  input  signal  is  inductively  coupled 
into  the  squid  loop,  its  origin  does  not  have  to  be  magnetic  in  nature. 
SQUID  amplifiers  have  been  used  with  many  other  types  of  inputs,  but  for 
our  concerns,  i.e.  biomagnetic  signals,  the  input  is  provided  by  the  detector 
loop  shown  as  Ld  in  Figure  3.4.5. 

Numerous  types  of  detector  coil  configurations  have  been  used  as  shown 
in  Figure  3.4.6.  The  most  popular  of  these  at  the  moment  is  the  second 
order  gradiometer.  This  configuration  cancels  out  the  effect  of  constant 
fields,  e.g.  earth’s  magnetic  field,  and  slowly  changing  fields,  e.g.  nearby 
equipment  effects.  It  is  the  use  of  this  detector  loop  that  allows  SQUID 
measurements  to  be  made  in  a  normal  lab  environment  without  shielding. 
Shielded  rooms  (Figure  3.4.4)  are  not  necessary,  unless  you  are  located  in  an 
extremely  magnetically  noisy  environment,  such  as  at  WPAFB  where  the 
experimental  data  for  this  research  was  collected. 

3.4.3.  Models  for  Evoked  Fields 

In  order  to  understand  the  nature  of  the  evoked  field  it  it  necessary  to 
construct  a  model  as  to  how  it  is  produced.  The  most  common  model  used 
is  that  of  considering  the  head  to  be  a  series  of  concentric  homogeneous 
spheres  and  the  response  to  be  a  simple  current  dipole.  It  has  been  shown 
for  this  model  that  the  radial  component  of  the  dipole  will  produce  no  net 
external  fields  (Cuffin  &  Cohen,  1977)  due  to  the  model  symmetry.  The 
strength  of  the  external  field  decreases  rapidly  as  the  dipole  is  located 
deeper  in  the  brain  (Aunon  et  al.,  1984).  A  dipole  located  at  the  center  of 
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Detection  coils  with  the  following  configurations:  a) 
magnetometer,  b)  first-order  pradiometer,  c)  off-diagonal 
pradiometer,  d)  second-order  pradiometer,  c)  asymmetrical 
first-order  pradiometer,  and  f)  asymmetrical  second-order 
pradiometer. 

Figure  3.4.6  Detector  Coil  Arrangements  (Williamson  &  Kaufman  1981) 
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the  brain  would  produce  no  external  field  since  it  must  be  radially  oriented. 
The  sharpness  of  the  peaks  in  the  field  are  much  more  pronounced  for  shal¬ 
low  dipoles. 

This  model  also  predicts  that  the  radially  oriented  component  of  the 
field  will  not  be  affected  by  the  volume  conduction  associated  with  the 
dipole.  The  tangential  components  of  the  field  will  be  affected  by  the 
volume  flow  and  consequently  they  are  not  often  measured.  The  detector 
loops  shown  in  Figure  3.4.6  are  used  for  measuring  the  radial  component 
only  and  in  general,  when  we  refer  to  the  evoked  fields  we  are  referring  to 
the  radial  component  only. 

The  EEG  measured  on  the  surface  of  the  scalp  represents  a  spatial 
smearing  of  the  true  activity  beneath  the  skull  due  to  the  high  resistivity  of 
the  skull  and  the  volume  conduction  flow.  The  EF  is  unaffected  by  the  scalp 
and  its  radial  component  should  represent  a  true  measure  of  the  activity 
beneath  the  skull.  This  means  that  the  EF  responds  to  a  much  more  local¬ 
ized  effect  than  the  ERP  and  details  should  show  up  more  clearly.  Lewis  et 
al.  (1985)  was  able  to  report  measuring  a  single  trial  I  OF  (no  averaging)  for 
the  first  time  by  placing  the  detector  coil  over  the  localized  activity.  EFs 
from  shallow  sources  can  be  very  sharp  and  only  a  few  centimeters  may 
separate  the  positive  and  negative  peaks  in  a  field. 

The  focus  of  much  of  the  research  to  date  has  been  solving  the  inverse 
problem,  i.e.  locate  the  source  from  the  measured  field.  This  has  met  with 
considerable  success  (Aunon  et  al.,  1984;  Cohen  &  Cuffin,  1983).  In  general 
the  contour  mapping  of  the  EF  is  offset  from  the  mapping  produced  by  the 
ERP.  It  has  been  found  that  good  EFs  can  be  recorded  from  the  sides  of  the 
head  (Okada  et  al.,  1981;  Okada,  1983),  while  most  ERPs  measurements  are 
usually  recorded  from  the  top  of  the  head.  This  finding  was  confirmed  while 
at  Wright  Patterson  and  a  measurement  location  for  the  experiment  was 
chosen  to  be  7  cm  above  the  right  ear.  This  location  was  in  good  accor¬ 
dance  with  the  result  obtained  by  Okada  for  maximal  amplitudes  in  the  300 
to  450  ms  latency  range. 


CHAPTER  4 


EXPERIMENTAL  PROCEDURE 


4.1.  Introduction 

This  experiment  basically  builds  upon  the  groundwork  laid  out  by  Kan- 
towitz  (1985b  Experiment  1).  The  time  frame  for  his  original  work  only 
allowed  for  a  preliminary  (phase  one)  analysis  of  the  data.  In  the  current 
research,  more  sophisticated  data  analysis  is  performed.  Additional  elec- 
toencephalogram  electrodes  and  the  magnetoencephalogram  have  been 
added  to  the  data  set.  There  are  several  advantages  to  proceeding  in  this 
fashion.  The  experiment  is  well  designed  and  provides  for  converging  opera¬ 
tions.  The  task  differences  are  fairly  simple  so  that  the  differences  in  load¬ 
ing  are  clearly  defined.  Several  aspects  of  the  workload  are  varied.  The 
experiment  has  been  shown  to  give  reliable  results  and  results  can  be  com¬ 
pared  to  the  previous  research. 

4.2.  Paradigm 

The  PRP  paradigm  presents  two  stimuli,  preceded  by  a  warning  tone, 
in  a  sequence  to  form  a  single  trial.  First  a  (W)  warning  tone  (70  dB,  lKHz 
tone  of  500  msec  duration)  alerts  the  subject  to  the  start  of  a  new  trial.  A 
random  foreperiod  (W  offset  to  SI  onset)  with  a  mean  of  2.0  sec  was  created 
by  sampling  foreperiods  of  1.5,  2.0  and  2.5  secs  with  equal  probability.  Two 
stimuli  (LEDS)  were  presented  with  a  fixed  inter  stimulus  interval  (ISI)  of 
62.5  msec  and  250  msec.  A  baseline  condition  where  a  second  stimuli  did 
not  follow  the  first  was  also  recorded.  Figure  4.2.1  shows  a  diagram  of  the 
timing  for  a  single  trial. 

The  LEDS  were  drawn  with  equal  probability  from  either  a  four  choice 
or  two  choice  case.  The  background  illumination  was  dim,  _■  0.01  ft-L,  and 
each  LED  was  individually  balanced  to  have  an  illumination  of  0.3  ft-L.  A 
block  consisted  of  51  such  trials  and  the  number  of  choices  and  ISI  was  con¬ 
stant  for  each  block.  The  first  trial  from  each  block  was  discarded  as  a 
preventive  measure.  The  intertrial  interval  (W  onset  to  W  onset)  was  6 
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secs,  so  that  each  block  lasted  5  minutes.  The  subject  was  informed  of 
the  block  type  (ISI  and  number  of  choices)  before  the  start  of  each  block. 

The  experiment  was  performed  at  the  Aerospace  Medical  Research 
Laboratory  (AMRL)  at  WPAFB  due  to  the  availibility  of  a  DC  SQUID.  The 
subject  sat  in  a  shielded  room  (Figure  4.2.2),  viewing  a  set  of  4  horizontally 
mounted  LEDS,  slightly  below  eye  level,  presented  through  a  small  hole  in 
the  inner  shielded  wall.  A  head  support  with  a  strap  was  used  to  keep  the 
head  steady  for  the  SQUID.  The  subject  viewed  the  lights  from  a  distance 
of  122  cm.  This  was  the  maximum  possible  distance  due  to  the  size  of  the 
room  and  SQUID  support  system.  The  LEDS  were  placed  close  together,  4 
cm  span,  so  that  the  entire  display  was  considered  foveal. 

Responses  were  made  on  four  piano-like  keys,  operated  by  the  second 
and  third  fingers  of  each  hand.  The  switches  were  balanced  to  have  the 
same  approximate  characteristic  feel  and  throw.  Correct  responses 
extinguished  the  appropriate  LED  immediately,  while  incorrect  responses 
extinguished  the  LED  at  the  end  of  the  sampling  interval  for  that  trial. 
Each  switch  could  only  be  used  once  in  a  trial.  On  two  choice  trials  only  the 
inner  two  lights  were  used,  so  that  each  response  was  from  a  different  hand. 
This  was  not  necessarily  true  for  the  four  choice  trials. 

A  complete  experiment  (Figure  4.2.3)  consisted  of  a  initial  relaxation 
period,  a  series  of  4  blocks,  another  relaxation  period,  another  series  of  4 
blocks,  and  then  a  final  relaxation  period.  Relaxation  periods  were  5 
minutes  long,  during  which  time  the  doors  to  the  chamber  were  opened  and 
the  strap  holding  the  head  in  place  for  the  SQUID  was  released.  The 
number  of  choices  was  invariant  for  each  set  of  4  blocks,  and  half  of  the 
subjects  were  assigned  to  perform  the  four  choice  set  first,  while  the  other 
half  performed  the  two  choice  set  first.  The  4  blocks  within  a  set  consisted 
of  a  single  stimulation  practice  block  first  and  the  3!  orders  of  ISI  and  base¬ 
line  presented  in  a  random  order.  A  total  of  six  lab  personnel  were  used  as 
subjects  with  one  randomly  assigned  to  each  order.  It  took  ~  one  hour  to 
prepare  each  subject  (attach  electrodes,  etc)  and  one  and  a  half  hours  to 
perform  the  experiment.  Table  4.2.1  contains  a  summary  of  the  subjects 
that  were  used  in  this  experiment  and  Table  4.2.2  gives  the  order  of  the  trial 
blocks  for  each  subject. 


Relaxation  1 


Relaxation  3 


Sets  are  4  choice  or  2  choice 


SI  trial*  per  block 


ISIs  can  be  60  ms,  240  ms  or 


no  second  light 


Figure  4.2.3  Experiment  Sequence 


Table  4.2.1  Subject  Summary 


SUB  5 

SUB  6 

SUB  7 

SUB  8 
SUB  9 


WF  35  (lab  person)  Full  disk  caused  loss  of  several  blocks  of  data.  ! 

Lost  blocks  were  rerun  to  make  a  COMPLETE  DATA  SET  minus  data 
for  the  two  practice  blocks  which  were  not  rerun  in  interest  of  time. 


WF  28  (lab  person)  COMPLETE  DATA  SET  recorded.  Subject  was  less 

than  willing  at  times  and  data  is  highly  corrupted  by  eyeblinks. _ 

WM  27  (lab  person)  COMPLETE  DATA  SET  recorded  plus  1  additional 

block  with  51  trials  of  4  choice  125ms  ISI  stimulation. _ 

WM  25  (lab  person)  COMPLETE  DATA  SET  recorded.  Subject  was  well 

practiced  with  task  before  start  of  this  run. _ 

WM  28  (lab  person)  COMPLETE  DATA  SET  recorded _ 

WM  26  (lab  person)  COMPLETE  DATA  SET  recorded.  Subject  was  well 
practiced  with  task  before  start  of  this  run.  Two  additional 
blocks  with  25  trials  of  2  choice  single  stim  data  were  recorded  ! 

with  the  SQUID  in  alternate  locations  (Files  2,0x  &  2_0xx).  I 


Table  4.2.2  Test  Sequence  for  Subjects 


FILE 


testOO 


testO 


test4 

test5 

test6 _ 

test  7 _ 

testS 
test  9 
test  10 
testl 1 
test  12 
test  13 
testl4 


SUB  6 


SUB  7 


SUB  8 


SUB  9 


_2-0 

2_60  (45) 
2_240  (49) 
rest2 


4  240  II 


4  _80  ILL 

restl 
2  (36) 

2,0  (41) 

2,240  (34) 
2,60  (32) 
rest2 


2,240  49 


2,0  50 


2_60  (49 
restl 


2,60  (50 
48 


2_240  (51) 
restl 


4,240  48 


4,60  42 


4,0  (50) 
restl 


J _ j 

[49)  4  | 

[50) 

_ ? _ [47} 

4,240 

(50)  4,80 

(51) 

2,0  (30) 

4,0 

(49)  4,0 

(50) 

2,240  (43 

4,60 

(47) 

2,60  (36) 

4,240 

(50) 

rest2 

rest2 

4,120  (49) 


2  240  48 


2.60  (42 

2,0  (49' 

restl _ 

_! _ (47) 

4,240  (43 

4.60  (35 

4,0  (35 

rest 2 _ 

2_0x  (22 
2_0xx  (25 


Number  in_[  ^indicates  number  of  valid  trials  if  the  entire  2.5  secs  of  data  is  to  be  used 
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For  each  trial,  the  stimuli  (LED  numbers),  times  of  stimuli,  responses 
(switch  numbers),  and  times  of  each  response  was  recorded.  Times  were 
recorded  to  the  nearest  0.1  ms  relative  to  the  start  of  the  block.  It  was 
decided  to  record  all  of  this  information,  rather  than  just  the  result  (correct 
or  incorrect)  of  each  trial  as  a  precautionary  move.  If  unusual  results  occur 
for  a  particular  trial,  LED  or  switch  it  would  be  beneficial  to  have  this  infor¬ 
mation  to  see  if  a  pattern  was  developing. 

4.3.  Physiological  Measurements 

4.3.1.  Overview 

Event  related  potentials,  evoked  fields,  respiratory  and  cardiac  rate 
were  recorded  simultaneously  during  the  experiment.  To  the  best  of  the 
author’s  knowledge,  this  was  the  first  time  that  anyone  had  attempted  to 
record  the  magnetoencephalogram  while  simultaneously  recording  these 
other  physiological  measurements.  Four  channels  of  EEG  data  (C3,  C4,  P3, 
P4),  1  eyeblink  channel  (EOG),  and  1  MEG  channel  (EF)  were  digitized  for 
2.5  secs  (at  240  Hz  sampling  rate  —  600  data  pts),  commencing  0.5  secs 
before  the  first  stimulation  of  each  trial.  The  EF  was  recorded  from  a  single 
location  7cm  directly  above  the  right  ear  canal,  as  described  in  section  3.4.4. 
Only  the  time  of  each  R  wave  of  the  ECG  and  the  time  of  each  expiration 
from  the  respiration  were  recorded  as  these  waveforms  were  not  digitized.  A 
Grass  Model  79  Polygraph  was  used  to  produce  a  paper  recording  of  C3,  C4, 
P3,  P4,  EOG,  EF,  ECG,  and  respiration. 

4.3.2.  Respiration  &  Heart  Rate 

The  ECG  was  measured  from  a  lead  II  configuration  (Guyton,  1981). 
The  respiration  was  determined  by  using  a  electrolytic  strain  gauge  (Geddes 
8c  Baker,  1975  chap.  2)  as  part  of  a  voltage  divider  circuit.  Both  signals 
were  then  amplified  using  Grass  7P511J  amplifiers  and  fed  to  the  data 
acquisition  interface  box. 

The  signals  were  threshold  detected  at  this  point  and  used  to  generate 
system  interrupts  to  record  the  time  of  each  R  wave  and  the  start  of  expira¬ 
tion.  Times  were  recorded  relative  to  the  start  of  the  experiment  and  to  the 
nearest  0.1  msec.  The  levels  for  threshold  detection  were  known  to 
correspond  to  a  specific  pen  deflection.  This  facilitated  setting  the  amplifier 
gains  and  correlating  the  paper  record  to  the  recorded  times  and  digitized 
EEG  when  necessary. 
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It  was  important  that  the  plane  of  the  loop  formed  by  the  electrolytic 
strain  gauge  was  kept  perpendicular  to  the  plane  of  the  SQUID  detector 
loop.  This  prevented  the  magnetic  field  produced  by  the  current  in  this 
loop  from  being  coupled  into  the  SQUID  channel.  The  current  was  limited 
to  several  microamperes  and  it  was  found  that  there  was  no  measurable  cou¬ 
pling  unless  the  plane  of  the  strain  gauge  was  parallel  to  the  detector  loop. 

Respiratory  and  cardiac  rate  were  also  recorded  during  the  rest  period 
between  each  set  of  blocks.  This  was  originally  done  in  order  to  check  for 
changes  in  baseline.  Unfortunately,  the  subjects  were  very  restless  during 
the  rest  periods  due  to  d.scomfort  of  the  head  restraint  during  the  stimulus 
blocks,  so  these  measurements  may  be  biased. 

4.3.S  Electroencephalogram  and  EOG 

ERP  data  was  collected,  using  Beckman  minature  electrodes,  from  the 
International  10-20  electrode  locations  (Jasper,  1958)  C3,  C4,  P3,  and  P4 
with  the  linked  mastoids  used  as  the  reference  electrode.  The  electroculo- 
gram  (EOG)  was  recorded  to  allow  for  off  line  rejection  of  eyeblink  artifact. 
Grass  amplifiers  7P511J  were  used  to  amplify  four  channels  of  EEG  data 
(P3,  P4,  C3,  C4)  and  one  eyeblink  channel  (EOG)  were  amplified  using  a 
Grass  7P511J  amplifiers  and  then  digitized  using  the  Tecmar  labmaster 
board.  These  amplifiers  as  well  as  the  Grass  amplifier  attached  to  the 
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SQUID  output  were  calibrated  in  terms  of  (A/D  units  out  of  the  digitizer)  / 
(micro  volts  in).  Offline  analysis  was  done  to  reject  trials  due  to  eyeblink 
and  amplifier  saturation. 

A  waveform  analyzer  was  used  to  determine  the  time  delay  associated 
with  the  Grass  amplifiers.  A  DC  -  20  Hz  linearly  swept  signal  was  input  into 
the  amplifier  and  the  resulting  output  was  then  correlated  with  the  input. 
The  peak  of  the  correlation  waveform  was  used  as  a  measure  of  the  delay 
time.  It  was  found  that  the  amplifiers  had  a  2.7  msec  delay  when  set  for  a 
.1  to  100  Hz  passband  and  a  4.7  msec  delay  when  the  60  Hz  notch  filter  was 
added  in.  Delay  times  are  important  since  the  electroencephlographic  data 
and  magnetoencephlographic  data  are  processed  by  different  channels  and 
their  digitized  waveforms  have  to  be  shifted  by  the  difference  of  their  delay 
times  in  order  to  preserve  the  time  synchronization  of  the  measurements. 
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4.3.4  Magnetoencephalogram  f 

Due  to  the  high  level  of  electrical  activity  within  the  building,  AMRL 
had  a  magnetically  shielded  room  built  by  EG&G.  The  room  (Figure  3.4.4) 
was  basically  a  double  layer  shielded  box  with  the  outer  layer  being  a  120 
inch  cube  and  the  inner  layer  a  90  inch  cube.  A  wood  floor  was  supported 
by  wooden  rods  passed  through  small  holes  in  the  bottom  of  the  chamber 
which  provided  mechanical  isolation  from  the  perma-alloy  walls  and  floor.  A 
double  set  of  sliding  doors  was  used  to  provide  access  to  the  shielded  room. 
When  closed,  doors  were  pneumatically  held  against  the  frame  to  ensure 
electrical  contact.  A  series  of  electrical  feedthroughs  were  used  to  send  sig¬ 
nals  in  and  out  of  the  chamber.  Several  removable  two  inch  portholes  were 
placed  in  the  walls  which  allowed  ventilation  and  the  use  of  non-standard 
cabling.  An  intercom  system  was  constructed  by  using  a  piezoelectric 
speaker  to  amplify  the  sound  passed  through  the  chamber  walls  by  a  hollow 
tube. 

A  simple  program,  "MUTUAL",  written  at  the  Naval  Personnel  and 
Development  Research  Center  (NPRDC)  in  San  Diego,  California,  (Appendix 
B)  was  used  to  compute  the  mutual  inductance  between  a  test  loop  and  the 
SQUID  detector  loops.  This  allowed  the  computation  of  the  magnetic  flux 
input  into  the  SQUID  based  on  the  current  in  the  test  loop  and  was  used  to 
verify  the  SQUID  calibration.  If  the  SQUID  is  functioning  properly,  the  cali¬ 
bration  from  the  test  loop  should  match  the  manufacturer’s  specification 
(Figure  B1  in  Appendix  B).  A  small  loop  of  wire  wrap  wire  was  placed 
around  the  tail  of  the  dewar  containing  the  SQUID.  The  remaining  wire  was 
twisted  to  prevent  stray  mutual  inductance  between  the  test  loop  and  the 
detector  loops. 

There  are  several  advantages  to  using  a  calibration  loop  instead  of  sim¬ 
ply  relying  on  the  manufacturers  calibration.  The  initial  calibration  may 
drift  if  the  alignment  of  the  probe  within  the  dewar  changed  due  to 
remounting  the  probe  or  allowing  the  dewar  to  warm.  SQUID  performance 
can  also  deteriorate  if  ice  should  buildup  under  the  coil  form  at  the  dewar 
tail.  This  has  happened  at  NPRDC.  Use  of  the  test  loop  allows  a  complete 
operational  performance  check  which  would  turn  up  any  unsuspected  prob¬ 
lems. 

Another  important  advantage  is  that  it  allows  one  to  compute  the 
entire  system  transfer  function  rather  than  just  a  part  of  it.  It  was  found 
that  the  delay  time  from  the  detector  loop  to  the  Grass  Amplifier  output 


was  22.5  msec  (SQUID  control  unit  =  3.3  msec,  Khronhite  25  Hz  Low  Pass 
Filter  =  14.5  msec,  Grass  Amp  (with  60  Hz  Notch)  =  4.7  msec).  The 
difference  in  time  between  the  SQUID  channel  and  the  electrode  channels 
was  19.8  msec  (22.5  -  2.7)  which  is  ~  5  sample  points  (at  4.167  ms  per 
point).  This  time  shift  has  largely  been  ignored  by  other  researchers,  though 
it  can  be  quite  significant  depending  on  the  amount  of  low  pass  filtering  that 
is  added  to  the  SQUID  channel.  The  filtering  used  at  NPRDC  amounted  to 
^  30  msec  of  delay  which  would  make  peak  alignment  between  the  EEG 
and  MEG  quite  difficult  if  this  went  uncorrected. 

Any  metal  that  the  subjects  were  wearing  or  carrying  was  first  removed 
before  they  entered  the  shielded  chamber.  The  subjects  wore  a  nylon  cap 
over  the  electrodes  and  a  circle  was  drawn  on  this  cap  to  aid  in  SQUID 
placement.  The  subjects  sat  with  their  head  strapped  into  a  padded  head 
rest  at  a  ~  45  degree  angle.  The  SQUID  was  lowered  into  place  until  it 
just  made  contact  with  the  cap  and  then  was  backed  off  slightly  (see  Figure 
4.2.2).  Since  the  shielded  room  and  SQUID  support  at  WPAFB  had  not 
been  used  for  experimentation  before,  it  was  necessary  to  try  several  types 
of  subject  placement  before  deciding  on  the  final  arrangement. 

Due  to  the  distraction  and  time  it  took  to  open  and  close  the  pneumati¬ 
cally  sealed  doors,  the  doors  were  only  opened  during  the  rest  periods 
between  the  sets  of  4  blocks  and  not  between  each  block.  This  meant  the 
subject  was  sealed  in  the  room  for  periods  of  ~  30  minutes  until  the  doors 
were  opened  during  the  rest  periods.  The  ventilation  system  was  inadequate 
to  keep  the  chamber  well  ventilated  over  such  a  long  period  of  time.  The 
seating  and  head  support  were  also  not  adequately  comfortable  for  the  long 
durations  inside  the  chamber.  One  must  consider  that  there  is  some 
refinement  period  after  the  construction  of  a  new  system  before  it  can  reach 
it’s  optimal  performance. 

4.4.  Data  Acquisition  System 

An  extremely  powerful,  yet  flexible  data  acquisition  system  (Figure 
4.4.1)  was  built  for  this  experiment  and  to  serve  the  future  needs  of  the  EEG 
laboratory  at  Purdue  University.  The  system  was  built  around  an  IBM  AT, 
equipped  with  internal  and  external  hard  disk  drives  and  a  Tecmar  Labmas- 
ter  data  acquisition  board.  A  Grass  Model  79,  8  channel  polygraph  and 
EEG  amplifier  was  used  to  amplify  and  make  paper  recordings  of  the  sig¬ 
nals.  A  complete  description  of  the  system  along  with  all  the  applicable 
software  and  documentation  is  given  in  Appendix  C. 


Figure  4.4.1  Equipment  Configuration 


An  interface  box  was  made  connecting  the  IBM  and  Tecmar  Labmaster 
to  the  outside  world.  This  interface  brought  all  of  the  Tecmar  data  and 
control  lines  to  a  standard  printed  circuit  (PC)  board  edge  connector.  The 
interface  also  supplied  power  and  made  22  output  lines  available  via  a 
choice  of  three  standard  types  of  connectors  mounted  on  the  outside  of  the 
box.  All  of  the  necessary  hardware  for  this  experiment  was  then  built  on  a 
single  PC  board  and  inserted  in  the  box.  The  advantage  of  this  type  of  sys¬ 
tem  is  that  its  flexible  nature  allows  many  different  experiments  to  be  built 
without  worrying  about  all  the  connections.  Each  new  experiment  is  simply 
constructed  on  a  PC  board  and  all  the  necessary  connections  are  already 
provided. 

The  data  acquisition  program  was  written  primarily  with  the  intent  of 
being  used  for  real  time  processing,  however  this  experiment  did  not  demand 
real  time  processing.  The  start  of  each  trial  is  initiated  by  the  controlling 
program  and  the  entire  data  collection,  stimulus  presentation  and  response 
recording  is  interrupt  driven.  A  flag  is  raised  at  the  end  of  the  trial  so  that 
the  controlling  program  knows  to  start  the  next  task.  The  controlling  pro¬ 
gram  is  free  to  handle  other  task  while  data  is  being  collected.  The  advan¬ 
tage  of  this  method  is  that  only  5%  of  the  systems  resources  are  used  in  the 
data  acquisition  and  95%  of  the  resources  are  available  for  real  time  pro¬ 
cessing.  The  data  acquisition  code  is  extremely  flexible  and  is  able  to  handle 
a  variety  of  tasks  (e.g.,  turning  oil  the  light  if  the  correct  switch  is  pushed) 
in  a  background  mode  that  can  respond  to  both  hardware  and  software 
interrupts. 

All  data  collection  is  double  buffered.  This  means  that  while  the  inter¬ 
rupt  driven  code  is  collecting  the  data  for  the  current  trial,  all  results  from 
the  previous  trial  are  available  to  the  controlling  program.  The  controlling 
program  can  process  that  data  and  a  decision  can  be  made  as  to  whether 
that  data  should  be  stored.  The  interrupt  driven  code  checks  to  see  if  a 
reject  flag  has  been  raised  by  the  controlling  program  before  storing  the 
data  on  the  next  trial,  hence  the  name  double  buffered.  The  time  frame  for 
this  research  did  not  allow  the  implementation  of  real  time  error  checking, 
however  all  of  the  mechanisms  for  adding  this  checking  to  the  package  are 
already  in  place.  All  of  the  design  problems  were  worked  out  at  Purdue  and 
then  the  entire  system  was  taken  to  WPAFB,  to  incorporate  their  SQUID 
into  the  design. 


CHAPTER  6 


PRELIMINARY  DATA  ANALYSIS 

6.1.  Introduction 

This  chapter  discusses  the  procedures  and  results  from  the  preliminary 
data  analysis  phase.  The  EEG  and  MEG  channels  are  screened  for  artifact 
rejection  and  then  filtered.  Average  responses  and  power  spectral  estimates 
from  the  prestimulus  interval  are  computed.  Mean  inter-beat-intervals  for 
the  heart  and  two  measures  of  heart  rate  variability,  sample  variance  and 
mean  square  successive  difference,  are  computed.  Statistics  based  on  subject 
performance  are  also  computed.  Data  from  six  subjects  are  used  in  this 
phase  of  the  analysis. 

6.2.  Behavioral  Responses 

The  response  performance  is  used  to  provide  the  converging  criteria 
called  for  by  Kantowitz  (1985a).  The  "RESPOND"  program  was  used  to 
compute  the  response  performance.  For  each  trial,  the  response  fell  into  one 
of  six  categories: 

1)  The  subject  did  not  respond  to  the  presented  stimuli. 

2)  The  subject  responded  before  the  stimulus  was  presented. 

3)  Both  responses  were  wrong. 

4)  First  response  was  wrong,  second  response  was  correct. 

5)  First  response  was  correct,  second  response  was  wrong. 

6)  Both  responses  were  correct. 

The  program  counted  the  number  of  trials  that  fell  into  each  of  these 
categories  and  computed  the  average  response  times  and  standard  deviations 
for  the  last  four  categories. 

In  this  report,  the  term  "correct  response"  shall  be  taken  to  mean  that 
both  responses  were  correct  unless  specifically  stated  otherwise.  Since  the 
firs*  trial  was  always  rejected  as  practice,  there  was  a  maximum  of  50 
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possible  correct  responses.  We  shall  only  consider  the  response  time  statis¬ 
tics  for  the  "both  correct"  case  since  the  majority  of  the  responses  fall  in  this 
category  and  incorrect  responses  do  not  appear  to  have  the  same  mean  as 
correct  responses.  RTl,  RT2  represent  the  average  response  times  and 
<JrTi  >  0RT2  represent  the  standard  deviations  of  the  response  times  to  the 
first  and  second  stimulus  respectively.  In  comparing  average  response  times 
between  loading  conditions  for  a  given  subject,  it  is  important  to  consider 
the  standard  deviations  of  the  response  times.  In  Figure  5.2.1,  bar  graphs 
are  used  to  compare  response  times  to  the  first  stimulus.  The  upper  part  of 
each  bar  is  shaded  to  indicate  the  range  corresponding  to  the  average  plus 
and  minus  the  standard  deviation.  The  midpoint  of  the  shading  is  the  aver¬ 
age  response  time  for  that  loading.  Figure  5.2.2  is  a  similar  graph  for  the 
response  to  the  second  stimulus.  Figure  5.2.3  is  a  multiple  bar  graph  which 
shows  the  percentage  of  correct  responses  to  both  stimuli  for  all  subjects  at 
all  loading  levels.  In  most  cases  the  accuracy  was  in  excess  of  90%,  i.e. 
more  than  45  out  of  50  correct  responses  to  both  stimuli. 

One  of  the  ways  that  loading  can  be  determined  is  from  the  task 
description.  The  goal  of  this  research  is  to  find  behavioral  and  physiological 
correlates  to  the  predicted  increase  in  loading.  We  can  evaluate  the  task 
descriptions  (i.e.  number  of  choices  and  ISI)  in  light  of  the  hybrid  capacity 
model  to  assign  relative  loadings.  There  are  15  possible  binary  comparisons 
for  the  six  loading  conditions,  however  not  all  of  the  outcomes  can  be 
predicted  from  the  model  and  task  description.  The  model  cannot  provide 
any  insight  as  to  the  relative  loading  between  a  2  choice  62.5  msec  ISI  load¬ 
ing  and  a  4  Choice  250  ms  ISI  for  example.  Table  5.2.1  summarizes  the 
predicted  outcomes  for  the  behavioral  variables  for  six  of  these  binary  com¬ 
parisons  based  on  the  model  and  task  description.  Three  of  these  comparis¬ 
ons  measure  the  effect  of  choice  by  holding  the  ISI  (or  single  stimulation) 
constant.  Two  of  these  comparisons  measure  the  effect  of  ISI  on  double 
stimulation  by  holding  choice  constant.  The  last  binary  comparison  is 
between  the  highest  (4  Choice  62.5  msec  ISI)  and  lowest  (2  Choice  250  msec 
ISI)  loading  for  the  double  stimulation  cases.  Note  that  the  hybrid  capacity 
model  does  not  specifically  require  a  slower  response  time  to  the  first 
stimulus  when  the  ISI  is  decreased  from  250  msec  to  62.5  msec,  however  it 
does  allow  for  this  to  occur  while  some  of  the  previously  discussed  models  do 
not. 

Table  5.2.2  summarizes  the  measured  outcomes  of  these  comparisons 
across  all  six  subjects.  Mean  response  times  to  both  stimuli  for  each  subject 
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were  compared  for  the  different  loading  cases  by  using  a  modified  student  t- 
test  which  allowed  for  different  variances  and  sample  sizes  (Sachs  p.  271, 
1984).  Only  the  trials  with  correct  responses  were  used  in  computing  the 
mean  and  variance  for  the  t-test.  The  response  time  entries  in  Table  5.2.2 
indicate  those  subjects  whose  mean  response  times  were  in  agreement  (at  the 
0.025  level  of  significance)  with  the  outcome  predicted  in  Table  5.2.1.  A 
means  test  could  not  be  performed  on  the  accuracy  since  scoring  for  each 
trial  was  not  based  on  a  numerical  value.  The  accuracy  entries  in  Table 
5.2.2  are  a  simple  comparison  of  the  mean  accuracies  for  each  subject.  One 
would  expect  that  response  time  2  and  accuracy  would  be  the  best  indicator 
of  loading  for  the  double  stimulation  cases  and  that  response  time  1  and 
accuracy  would  be  the  best  indicators  for  the  single  stimulation  control 
cases.  It  appears  that  at  the  higher  levels  of  loading  (columns  3,5,6),  the 
measured  performance  is  consistent  with  the  predicted  performance.  At  the 
lowest  loading  level  (column  l),  the  measured  performance  is  not  consistent 
across  the  subjects  for  either  accuracy  or  response  time. 

Our  task  description  and  behavioral  measures  have  resulted  in  the  same 
effective  loadings  for  the  higher  loading  comparisons.  What  remains  is  for 
the  physiological  measures  to  also  converge  on  the  same  relative  loadings. 
This  would  provide  the  converging  criteria  that  are  called  for  by  Kantowitz 
(1985a). 

5.3.  Heart  Rate 

The  mean  IBI  and  two  measures  of  its  variance  were  computed  in  this 
phase  of  the  analysis.  The  data  acquisition  program,  "DAC",  recorded  the 
time  of  each  R  wave  as  determined  by  threshold  detecting  the  ECG.  It  was 
necessary  to  first  prescreen  this  data  to  check  for  ectopic  beats  and  shifts  in 
baseline  which  might  cause  a  P  wave  to  cross  the  threshold.  The 
"INTBEAT"  program  was  used  to  produce  two  output  files,  one  containing 
the  actual  time  of  each  beat  and  the  other  containing  the  inter-beat- 
intervals  (IBI).  The  program  would  also  flag  beats  whenever  the  IBI  was  not 
within  an  acceptable  range  (500  -  1500  ms).  Plots  of  the  EBI  were  then  made 
for  each  loading  case  and  examined  visually  for  unusual  deviations  that  the 
computer  program  may  have  missed.  The  statistics  for  the  IBIs  were  then 
computed  from  the  corrected  data. 

Mean  heart  IBI  ranged  from  a  low  of  757  msec  (subject  5)  to  a  high  of 
1218  msec  (subject  9)  as  computed  over  the  5  minute  interval  that  it  took  to 
complete  one  trial  block.  The  threshold  detector  was  set  at  ~  30%  of  the  R 
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Figure  5.3.1  Mean  Heart  Inter-Beat-Interval  (FBI) 


Table  5.2.1  Behavioral  Performance  Predicted  by  Task  Description 
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Table  5.2.2  Measured  Behavioral  Performance 
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wave  amplitude,  to  favor  an  extraneous  threshold  crossing  (c.g.  P  wave 
crossing)  over  a  missed  R  wave.  Only  one  recorded  IBI  was  found  to  exceed 
1500  msec  and  the  paper  recording  of  the  ECG  showed  it  to  be  due  to  a  long 
beat.  It  is  extremely  unlikely  that  there  were  any  missed  beats  due  to  the 
setting  of  the  threshold  level  and  that  the  INTBEAT  program  and  visual 
screening  of  the  data  did  not  indicate  any. 

There  were  several  instances  were  a  recorded  IBI  was  measured  below 
500  ms.  These  questionable  beats  were  evaluated  on  an  individual  basis  by 
examing  the  paper  record  of  the  ECG.  Ectopic  beats  which  inhibited  the 
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next  regular  beat  were  simply  moved  to  the  mid  point  between  the  two  adja¬ 
cent  beats.  This  can  be  justified  since  ectopic  beats  are  of  local  origin  to  the 
cardiac  tissue  and  are  not  regulated  by  the  autonomic  mechanisms.  P  wave 
crossings  were  removed  entirely.  Only  one  P  wave  crossing  (across  all  six 
subjects)  was  found  during  the  testing  blocks,  though  several  others  were 
found  during  the  rest  periods  when  the  subjects  tended  to  move  around 
more.  The  only  ectopic  beats  found  were  for  subject  six,  who  had  5  ectopic 
beats  out  of  a  total  of  ~  2000  beats. 


The  "HR_AVAR"  program  was  used  to  compute  the  mean  inter-beat- 
interval,  FBI,  and  two  measure  of  the  heart  rate  variability,  s2  and  c>2/2  from 
the  corrected  IBIs.  s2  is  the  sample  variance  which  provides  an  unbiased 
estimate  of  the  population  variance,  /2  is  Vz  mean  square  successive 

difference  (MSSD)  and  also  provides  an  unbiased  estimate  of  o' jgj.  These 
values  were  computed  in  the  following  fashion: 


EBI  =  —  Vffil,  s2  =  — 


VIBI;2  -  (IBI)2/n  Vj  (EBIi+1  -  IBI,)2 

T -  <?  j2  =  '=1  „ - 


n  —  1 


2(n  -  1) 


It  has  been  claimed  that  ^ /2  is  a  better  estimate  of  heart  rate  variability 
(HRV)  than  s2  (Heslegrave  et  al,  1979;  Sharit  and  Salvendy,  1982)  since  the 
former  removes  the  effects  of  nonrandom  gradual  trends. 


Figure  5.3.1  shows  the  IBI  for  all  subjects  and  loading  cases.  There 
does  not  appear  to  be  any  consistent  pattern  of  IBI  with  loading.  The  IBI 
remains  constant  for  some  subjects,  increases  with  loading  for  other  subjects 
and  decreases  with  loading  for  others.  Figure  5.3.2  shows  the  sample  stan¬ 
dard  deviation,  \/s2  and  Figure  5.3.3  shows  the  standard  deviation  com¬ 
puted  from  the  MSSD  method,  v^/2.  Again,  no  consistent  pattern 
emerges  to  distinguish  the  workload.  Note  that  /2  is  consistently  smaller 
than  s2,  with  the  exception  of  subject  9  four  choice  loading  conditions.  It 
appears  that  there  may  be  trends  present  which  inflate  the  value  of  s2. 

s2  should  match  /  2  exactly  for  independent  samples  from  a  normally 
distributed  population.  The  ratio  t)2  /2  /  s2  can  be  used  as  a  test  of  indepen¬ 
dence  for  normal  data  (Sachs  p.  373,  1984).  If  we  assumed  50  sample  points 
for  each  estimate,  then  this  ratio  would  have  to  be  less  than  .772  to  reject 
the  hypothesis  of  independence  at  the  .05  level  of  significance.  This  would 
correspond  to  a  threshold  of  .878  (VY772)  for  the  standard  deviation  esti¬ 
mates  which  are  shown  in  Figures  5.3.2  and  5.3.3.  A  quick  visual  inspection 
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Figure  5.3.2  Standard  Deviation  of  Heart  IBI 
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Figure  5.3.3  Standard  Deviation  of  Heart  IBI  (MSSI)  Method) 


of  these  two  graphs  show  that  the  LBIs  are  not  independent  with  the  possible 
exception  of  subject  9.  A  standard  run  test  (Sachs  p.  375,  1984)  was  per¬ 
formed  on  the  IBIs  and  the  hypothesis  of  independence  was  rejected  at  the 
.0001  level  of  significance  for  all  cases.  This  is  to  be  expected  since  the 
heart  rate  is  a  controlled  process. 

One  of  the  largest  factors  in  HRV  is  due  to  respiratory  sinus  arrhyth¬ 
mia  (SA).  The  IBI  should  stretch  with  each  expiration  and  shrink  with  each 
inspiration.  Figure  5.3.4  shows  a  plot  of  IBI  over  a  100  sec  time  interval 
during  the  middle  of  the  4  choice  -  62.5  msec  loading  block  for  Subject  7. 
The  start  of  each  expiration  is  indicated  on  this  graph  by  a  There 

appears  to  be  a  trend  of  increasing  IBI  following  the  start  of  expiration. 
Note  that  not  every  expiration  will  cause  the  EBI  to  increase  since  there  are 
many  other  factors  that  also  increase  IBI.  SA  will  be  dealt  with  more  in  the 
advanced  analysis  section. 

Figure  5.3.5  shows  IBI  over  the  course  of  an  entire  block  (same  subject 
and  block  as  above)  along  with  response  time  1  (RTl)  plotted  and  aligned 
on  the  same  time  scale.  It  does  not  appear  that  RTl  is  correlated  with  IBI. 
Also  the  IBI  starts  at  750  msec  and  quickly  rises  to  its  mean  value  1000  msec 
about  which  it  varies  for  the  remainder  of  the  block.  This  pattern  of  a  sud¬ 
den  increase  in  IBI  was  not  consistent  across  other  blocks  though  it  did 
appear  in  some. 

5.4.  Respiration 

The  respiration  channel  was  processed  in  a  similar  fashion  to  the  heart 
rate  channel.  The  data  was  first  prescreened  by  the  computer  to  flag  all 
breaths  that  did  not  fall  within  between  1.8  secs  and  8  secs.  These  breaths 
were  corrected  by  examining  the  paper  record.  The  resulting  data  was  then 
graphed  for  all  of  the  cases  and  was  visually  inspected  for  abnormal  varia¬ 
tions  that  might  need  correction. 

The  chest  circumference  was  the  variable  that  was  measured  as  an  indi¬ 
cator  of  respiratory  activity.  Unfortunately  this  measurement  does  not 
necessarily  make  monotonic  transistions  from  inspiration  to  expiration  so 
simple  threshold  detection  does  not  produce  the  most  reliable  estimate  of 
respiratory  rate.  A  significant  problem  was  partial  breaths  which  did  not 
always  cross  the  threshold  level.  Careful  selection  of  the  amplifier  gain  and 
threshold  level  for  each  subject  reduced  the  number  of  beats  that  had  to  be 
resolved  by  visual  inspection  to  approximately  \%  for  subjects  5,6,7  and  8, 
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Figure  5.3.4  Respiratory  Sinus  Arrhythmia  as  Detected  using  IBI 
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but  was  at  a  10%  level  for  subjects  4  and  9.  Subject  4  liked  to  hold  her 
breath  and  on  several  occasions  this  exceeded  14  seconds  in  length.  Subject 
9  liked  to  take  partial  breaths  which  did  not  cross  the  threshold  level.  In 
hand  correcting  the  data,  the  rule  that  was  used  was  that  a  partial  breath 
had  to  at  least  half  the  size  of  the  surrounding  breaths  in  order  to  be 
counted. 

Figure  5.4.1  shows  the  mean  respiratory  intervals  for  the  corrected 
data.  Figure  5.4.2  shows  the  standard  deviation  of  the  respiratory  intervals. 
Note  the  large  standard  deviation  for  Subject  4  due  long  periods  of  holding 
her  breath.  At  the  start  of  this  research  it  was  believed  that  it  was  only 
necessary  to  record  the  respiratory  rate  so  that  its  contribution  to  the  power 
spectrum  of  the  heart  rate  could  be  accounted  for.  In  retrospect  it  would 
have  been  advatageous  to  digitize  the  respiratory  waveform  so  that  the 
correlation  function  and  depth  of  breathing  could  be  considered. 

5.5.  Artifact  Rejection  and  Prefiltering 

The  "PERFORM"  program  is  used  to  determine  records  to  be  rejected. 
It  produces  two  output  files,  a  reject  file  and  a  performance  file.  The  reject 
file  is  a  list  of  records  to  be  rejected  along  with  a  code  number  which  indi¬ 
cates  the  reason  for  the  rejection.  The  performance  file  contains  the  subject 
response  times  along  with  a  performance  code  which  indicates  the  subjects 
accuracy  for  each  of  the  non-rejected  trials.  The  "ANALYZE. H"  file  con¬ 
tains  a  list  of  these  codes. 

A  trial  was  rejected  due  to  any  of  the  following  reasons: 

1)  If  it  was  the  first  trial. 

2)  If  an  eyeblink  occurred  during  the  response. 

3)  If  channel  saturation  occurred  during  the  response. 

4)  If  the  subject  did  not  respond  to  the  stimulus. 

5)  If  the  subject  responded  before  the  stimulus  was  presented. 

An  eyeblink  is  defined  as  a  change  exceeding  75  microvolts  on  the  electrocu- 
logram  (EOG)  channel  during  a  100  msec  interval.  Saturation  is  defined  by 
any  of  the  channels  (EEG,  MEG  or  EOG)  exceeding  the  input  range  of  the 
digitizer  for  more  than  2  sample  pts  (~  8  msec).  Table  4.2.2  shows  the 
number  of  valid  trials  for  each  block  if  the  entire  2.5  secs  of  data  was 
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Figure  5.4.1  Mean  Respiratory  Intervals 
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Figure  5.4.2  Standard  Deviation  Respiratory  Intervals 


searched  for  artifacts.  In  each  stage  of  the  analysis,  the  PERFORM  pro¬ 
gram  was  rerun  to  reject  only  those  artifacts  that  occurred  within  the  time 
intervals  that  would  be  used  in  that  analysis. 

The  "SEP_CHAN"  program  is  used  to  separate  and  filter  the  individual 
channels  from  the  original  data  file  which  was  recorded  by  the  "DAC"  pro¬ 
gram.  The  program  reads  the  reject  file  so  that  only  valid  records  are 
separated  out  and  used  in  computing  the  average.  The  program  has  the  abil¬ 
ity  to  undersample  the  data  and  also  produces  an  average  of  the  filtered 
data  for  each  channel. 

The  data  for  each  channel  is  first  calibrated  and  the  linear  trend  that  is 
formed  by  joining  the  first  and  last  point  is  then  removed.  This  trend  is 
removed  to  force  the  first  and  last  point  to  zero  and  remove  step  effects 
when  the  filter  is  applied.  Note  that  this  is  different  from  the  "best  fit  fine" 
formed  by  a  linear  regression  of  all  the  data  points.  A  lowpass  filter  designed 
using  the  Remez  Exchange  Algorithm  (McClellan  et  al  1979)  is  then  applied. 
The  filter  used  was  a  39  point  finite  impulse  response  filter  with  a  3  dB 
cutoff  at  25  Hz.  The  passband  ended  at  21  Hz  (.075  dB  ripple)  and  the  stop- 
band  started  at  30  Hz  (36  dB  attenuation).  The  data  was  then  undersam¬ 
pled  by  selecting  every  fourth  data  point.  This  reduced  the  effective  sam¬ 
pling  rate  to  60  samples  per  second  and  the  corresponding  nyquist  rate  to  30 
Hz. 

6.6.  Power  Spectral  Estimate  of  the  Prestimulus  Interval 

For  each  trial,  there  was  a  l/z  sec  prestimulus  interval  which  consisted  of 
30  data  points  after  undersampling.  The  "BT"  program  was  used  to  com¬ 
pute  the  power  spectral  estimate  (PSE)  from  this  interval.  The  mean  of 
each  trial  was  removed  and  the  autocorrelation  function  (equivalent  to  the 
autocovariance  for  a  zero  mean  process)  was  then  computed.  The  indivi¬ 
dual  autocorrelations,  Rxx,  were  then  averaged  to  produce  an  average  auto¬ 
correlation,  Rxx.  Rxx  was  then  normalized  by  dividing  by  Rxx(0).  A  30  point 
Hamming  window  was  applied  and  a  discrete  symmetric  cosine  transform 
was  taken  to  produce  a  normalized  (o2 =1)  estimate  of  the  power  spectrum 
(PSE)  with  a  frequency  spacing  of  2  Hz. 

If  we  consider  the  ongoing  EEG  (i.e.  prestimulus  interval)  to  be  a  ran- 

-  2 

dom  process,  then  Rxx(0)  is  an  estimate  of  the  power,  O  ,  in  this  process. 

This  estimate  has  a  mean,  <7  =  Rxx(0),  and  a  standard  deviation, 
0\2  =  <7r  (g).  Rxx(0)  is  the  scaling  factor  for  the  PSB1  and  Or  (q)  is  the 


standard  deviation  of  the  scaling  factor.  Appendix  A  contains  graphs  of  all 
of  the  normalized  power  spectra  obtained  in  this  fashion.  There  are  30 
graphs  (6  subjects  x  5  channels)  with  6  levels  of  loading  per  page  to  allow 
easy  comparisons  of  the  PSEs.  1^(0)  is  indicated  on  these  graphs  as  the 
"avg  power"  and  Or  ^  is  indicated  as  "std  dev". 

In  general,  two  cases  can  be  separated  by  a  classifier  if  the  difference 
between  the  means  of  the  estimator,  mj—  m2,  is  significantly  greater  than  the 
sum  of  the  standard  deviations  <7i+^2  (Figure  3.3.1).  Since  Or  jqj  is  typi¬ 
cally  on  the  order  of  Rxx(0),  the  power  in  the  prestimulus  interval  would  not 
provide  a  good  measure  of  the  mental  loading.  If  any  particular  frequency 
(or  frequency  band)  of  the  PSE  was  to  be  used  as  a  classifier,  then  the 
respective  PSEs  would  have  to  be  markedly  different  due  to  the  large  vari¬ 
ance  of  the  scaling  factors.  Looking  at  the  graphs  in  Appendix  A,  this  does 
not  appear  to  be  the  case.  It  may  still  be  possible  however  to  use  the  PSE 
as  an  indicator  of  workload  by  comparing  the  ratio  of  power  in  different 
bands,  i.e.  PSEJ2HZ)  /  PSEj(8Hz)  .vs.  PSE2(2Hz)  /  PSE2(8Hz).  Due  to 
the  large  amounts  of  data,  this  will  require  significant  analysis  before  trends 
may  appear  and  is  not  covered  in  this  phase. 

6.7.  Average  of  the  ERP  and  the  EF 

Appendix  A  also  contains  30  graphs  (six  subjects  x  five  channels)  of  the 
average  low  pass  filtered  ERPs  and  EF  with  six  loading  levels  per  page.  The 
time  scale  for  these  graphs  was  chosen  to  show  the  entire  data  record  ;  from 
500  msec  before  the  first  stimulus  was  presented  to  2000  msec  after  it.  In  all 
discussions  and  graphs  of  ERPs  and  EFs  in  this  report,  the  time  of  the  first 
stimulus  is  indicated  as  time  =  0,  with  the  prestimulus  interval  indicated  by 
negative  time.  Recall  that  the  first  step  in  the  data  analysis  was  to  check 
for  artifacts.  In  appendix  A,  the  entire  data  records  (-500  msec  to  2000 
msec)  were  checked  for  artifacts,  so  that  only  artifact  free  records  were  used 
when  computing  the  average. 

There  are  advantages  and  disadvantages  to  checking  the  entire  record 
for  artifacts.  The  purpose  of  the  graphs  in  appendix  A  was  to  show  a  base¬ 
line  condition  (prestimulus),  a  response  (ERP)  and  then  a  return  to  baseline. 
A  single  unremoved  eyeblink  anywhere  in  the  region  of  interest  can  have  a 
detrimental  effect  on  the  average.  The  danger  in  rejecting  trials  from  an 
artifact  anywhere  in  the  record  is  that  you  may  reject  more  of  the  data  than 
you  want.  Subject  5  was  someone  who  might  be  described  as  an 


uncooperative  subject.  In  the  high  level  loading  case  of  4  Choice  -  62.5  msec 
ISI,  there  were  only  6  unrejected  trials  out  of  a  possible  50.  If  the  interval 
for  checking  artifacts  was  reduced  to  a  1  second  interval  (0  to  1000  msec), 
the  number  of  valid  trials  increases  to  36.  If  the  region  of  interest  was  only 
the  response  (0  to  1000  msec),  36  valid  trials  would  yield  much  better  results 
than  6.  Throughout  the  data  analysis,  trials  are  rejected  only  if  the  artifact 
falls  within  the  region  of  interest  for  the  current  analysis.  Trial  rejections 
due  to  artifacts  (-500  to  2000  msec  range)  ranged  from  an  average  high  of  23 
rejections  per  block  (50  trials)  for  Subject  5,  to  less  than  1  per  block  for  Sub¬ 
ject  7,  with  Subject  6’s  data  being  artifact  free. 

Average  responses  are  normally  computed  by  aligning  the  stimulus  of 
each  trial  and  computing  an  average.  The  peaks  in  the  ERP  and  EF  which 
are  synchronized  to  the  stimulus  become  enchanced  by  this  average.  What 
about  peaks  that  may  be  synchronized  with  the  response  instead  of  the 
stimulus?  The  "AVER"  program  was  used  to  compute  unfiltered  averages 
which  were  aligned  with  the  stimulus,  first  button  push  (response)  and  the 
second  button  push.  Figure  5.7.1  shows  a  conventional  average  (stimulus 
aligned)  overlaid  on  an  average  computed  by  aligning  the  first  response. 
The  response  aligned  average  was  computed  by  shifting  each  trial  so  that  it 
had  an  effective  response  time  of  0  msec,  i.e.  the  first  response  was  shifted 
to  the  first  stimulus.  In  plotting  the  graph,  the  conventional  average  was 
shifted  forward  by  the  average  response  time  to  the  first  stimulus  and  the 
x-axis  was  scaled  accordingly  so  that  0  msec  still  represented  the  time  of  the 
first  stimulus.  Figure  5.7.1  is  taken  from  the  P3  channel  of  Subject  6  with  a 
2  Choice  -  240  msec  ISI  loading.  The  average  response  time  for  the  first 
stimulus  (RTl)  was  314  msec  with  a  standard  deviation  (^rtj)  of  58  msec. 
The  response  aligned  average  appears  to  have  a  slightly  higher  peak  in  the 
575  msec  range  and  shows  a  negative  and  positive  peak  in  the  800  msec 
range  which  do  not  appear  in  the  conventional  average.  It  is  too  soon  to 
draw  any  general  conclusion  from  this  during  the  preliminary  analysis  stage. 

Figure  5.7.2  is  identical  to  Figure  5.7.1  except  that  the  EF  (MEG)  chan¬ 
nel  was  used  instead  of  P3.  It  was  hypothesized  that  the  oscillations  in  the 
conventional  average,  centered  about  900  msec,  were  due  to  head  movement. 
The  oscillations  in  the  response  aligned  average  have  the  same  form  but  are 
50%  larger  in  amplitude.  The  MEG  has  been  shown  to  be  relatively 
impervious  to  eyeblink  contamination  due  to  it’s  localized  measurement 
area.  It  is  therefore  unlikely  that  this  oscillation  is  a  muscle  artifact  from 
the  button  push,  since  the  arm  was  resting  on  a  table  and  the  finger  only 
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had  to  move  a  small  amount.  It  appears  likely  that  in  a  rush  to  respond 
(speed  was  stressed  over  accuracy)  a  subject  might  cause  their  upper  body 
to  oscillate  slightly  when  hitting  the  buttons.  Recall  that  the  buttons  were 
large  piano  like  keys  which  only  required  a  small  light  throw  since  they  were 
attached  to  microswitches.  The  subject  rested  his  fingers  on  the  keys,  but 
the  sheer  size  of  the  keys  invited  the  subject  to  use  all  his  force  when  push¬ 
ing  down  on  the  keys. 

In  addition  to  possible  head  movement  when  a  subject  "pounced"  on  the 
response  keys,  the  head  restraint  could  not  keep  the  head  properly  posi¬ 
tioned  under  the  SQUID  with  a  high  degree  of  precision.  The  SQUID  was 
placed  at  the  beginning  of  each  set  and  checked  at  the  end  of  each  set  (4 
blocks  and  ~  30  min.  later).  The  worst  case  movements  were  measured  to 
be  a  2  cm  lateral  shift  and  a  1  cm  in-line  shift  (i.e.  dewar  tail  was  1  cm  from 
the  scalp).  Since  it  took  several  minutes  to  open  and  close  the  doors  to  the 
shielded  room,  it  was  not  practical  to  check  the  SQUID  at  the  end  of  each 
block.  There  appear  to  be  oscillations  in  many  of  the  EFs  and  it  does  not 
appear  at  this  point  that  any  useful  information  can  be  obtained  from  these 
measurements.  Time  constraints  on  this  research  did  not  allow  the  oppor¬ 
tunity  to  wait  for  a  more  optimal  support  structure  which  would  have 
reduced  this  effect. 


CHAPTER  6 


ADVANCED  ANALYSIS 


6.1.  Introduction 

This  chapter  discusses  the  procedures  and  results  from  the  advanced 
data  analysis  stage.  Spectral  techniques  are  used  in  analyzing  the  heart  IBI. 
The  latency  corrected  average  is  used  to  further  analyze  the  ERP  and  iden¬ 
tify  significant  peaks  in  the  response.  These  peaks  are  then  checked  for 
correlation  with  the  behavioral  measures  of  response.  All  six  subjects  are 
used  in  the  heart  rate  analysis,  but  only  three  subjects  are  used  for  the 
latency  corrected  averaging  due  to  the  computational  effort. 

6.2.  Spectral  Estimates  of  Heart  Rate  Variability 

In  the  previous  chapter,  heart  rate  variability  was  measured  in  terms  of 
the  mean  inter-beat  interval,  IBI  and  two  estimates  of  the  variance  of  the 
beat  interval,  <7  jgj.  In  this  chapter  we  break  down  the  heart  rate  variabil¬ 
ity  into  its  spectral  components  to  get  a  clearer  picture  of  the  effect  of  load¬ 
ing.  Several  methods  of  computing  the  Power  Spectral  Estimate  (PSE)  were 
tried  and  compared  before  deciding  to  use  the  PSE  computed  from  the  auto¬ 
correlation  of  the  Low  Pass  Filtered  Event  Series  (LPFES). 

There  are  two  main  ways  of  computing  the  the  spectral  content  for  a 
series  of  point  events.  The  spectrum  of  intervals  is  based  on  computing  the 
power  spectrum  of  the  inter-beat  intervals  using  the  beat  number  as  a 
transform  axis  while  the  spectrum  of  counts  is  based  on  using  time  as  a 
transform  axis  (Dekwaadsteniet,  1982;  Sayers,  1973,  1975).  In  computing  the 
spectrum  of  counts,  one  can  either  work  with  a  series  of  unevenly  spaced 
delta  functions  (R-waves)  or  interpolate  to  an  evenly  spaced  time  series  of 
inter-beat  intervals.  Mohn  (1975)  discusses  the  computational  complexity  of 
the  various  methods.  Rompelman  et  at.  (1982)  advocates  interpolating  by 
rounding  off  the  location  of  the  R-wave  to  the  nearest  msec  to  preserve  the 
binary  nature  of  the  signal  and  then  computing  the  PSE  by  use  of  a  sparse 
discrete  fourier  transform.  Kantowitz  (1985b)  simply  performed  a  linear 
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interpolation  from  IBI  to  instantaneous  heart  rate  at  1  sec  spacing.  The  sig¬ 
nal  can  also  be  interpolated  by  passing  the  delta  functions  corresponding  to 
the  R-waves  through  an  ideal  low  pass  filter,  which  amounts  to  a  Sine  func¬ 
tion  (Sin(x)/x)  interpolation.  The  French-Holden  algorithm  (Peterka  et  ah, 
1978).  is  a  computationally  efficient  means  of  applying  an  ideal  low  pass 
filter  with  a  frequency  cutoff  at  V2  the  desired  sampling  rate  by  utilizing  the 
binary  nature  of  the  series  of  point  events  The  output  of  this  filtering  has 
been  called  the  low-pass  filtered  event  series  (LPFES).  DeBoer  et  al.  (1984) 
has  shown  that  for  small  variations  from  the  mean,  the  spectrum  of  intervals 
and  the  spectrum  of  counts  will  yield  very  similar  results  if  all  factors  are 
taken  into  account. 

One  problem  when  computing  the  PSE  from  the  lBIs  is  that  the  data 
from  each  block  does  not  have  the  same  number  of  points.  Each  data  block 
consist  of  ~  311  seconds  of  data  which  may  contain  anywhere  from  a  low  of 
254  heart  beats  to  a  high  of  404  heart  beats  depending  on  the  subject  and 
loading.  The  first  approach  investigated  was  to  compute  the  spectrum  of 
intervals  for  the  first  and  last  128  beats  in  each  trial  in  an  attempt  to  see  if 
the  spectrum  was  changing  with  time.  This  was  computed  by  taking  the 
magnitude  of  the  Fast  Fourier  Transform  (FFT)  of  these  two  data  segments 
after  multiplying  by  a  Hamming  window  and  removing  the  mean.  While  the 
two  PSEs  both  identified  peaks  in  the  same  location,  the  variance  in  the 
estimates  were  too  large  to  draw  any  specific  conclusions.  For  normally  dis¬ 
tributed  data,  the  variance  of  the  estimate  of  the  PSE  computed  in  this 
fashion  can  be  as  large  as  the  estimate  itself.  To  reduce  this  variance,  the 
data  is  usually  segmented  and  the  PSEs  from  each  segment  are  averaged  to 
compute  the  periodogram  (Kay  &  Marple,  1981).  This  reduction  in  variance 
is  obtained  at  the  cost  of  reducing  the  frequency  resolution. 

There  are  still  several  limitations  with  proceeding  in  this  fashion.  The 
number  of  data  points  for  each  case  are  still  not  equal  so  that  the  processing 
and  variance  for  each  estimate  will  be  slightly  different.  If  it  is  desired  to 
view  the  PSE  on  a  frequency  (Hz)  axis  by  scaling  the  spectrum  of  intervals 
axis  by  1  /IBI,  then  each  case  will  have  a  different  delta  frequency  (df),  the 
measurement  between  frequency  components,  and  a  point  by  point  com¬ 
parison  between  the  cases  may  not  be  accurate.  We  can  not  break  the  data 
into  many  segments  since  we  need  a  resolution  fine  enough  to  distinguish 
between  the  energy  in  the  thermal  band  (0.025  Hz)  and  the  blood  pressure 
band  (0.07  -  0.1  Hz)  as  shown  in  Figure  3.1.1.  A  frequency  increment  of  0.01 
Hz  in  the  PSE  would  correspond  to  a  100  point  data  segment  if  the  IBI  was 


equal  to  1  sec.  If  the  processes  which  mediate  the  heart  rate  are  functions 
of  time  rather  than  beat  number,  their  energy  in  the  transform  domain  will 
be  spread  out  due  to  the  nonuniform  sampling  of  a  periodic  signal. 

The  spectrum  of  counts  approach  was  deemed  to  have  too  many 
shortcomings  and  was  abandoned  after  a  cursory  look  did  not  indicate  any 
consistent  differences  between  the  loading  cases  or  the  PSE  from  the  first 
and  last  128  beats  in  a  trial.  It  was  decided  to  pursue  a  spectrum  of  counts 
approach  based  on  the  French-Holden  algorithm  to  provide  the  interpola¬ 
tion.  The  output  of  the  French-Holden  algorithm  was  inverted  to  convert  it 
from  a  rate  to  a  set  of  interpolated  equal  spaced  IBls  which  was  then  used 
as  the  LPFES.  The  first  10  and  last  10  points  were  discarded  to  remove 
interpolation  effects  and  any  initial  change  in  heart  rate  with  the  start  of  a 
block.  The  mean  was  then  removed  from  the  data.  The  autocorrelation 
was  then  computed  and  windowed  by  a  Hamming  window  to  64  positive  lag 
values.  The  Fourier  transform  of  the  windowed  autocorrelation  function  was 
then  used  as  the  PSE  for  analyzing  the  heart  rate  variability.  This  is  the 
classical  Blackman  Tukey  PSE  (Kay  &  Marple,  1981). 

In  order  to  fully  understand  the  effects  of  the  interpolation,  the  PSE 
was  computed  using  the  Blackman  Tukey  method  from  both  the  original  1BI 
data  and  the  equally  time  spaced  LPFES.  Figure  6.2.1  shows  the  original 
IBI  series  superimposed  on  the  LPFES  series  at  two  different  time  scales  for 
Subject  6  at  a  2  Choice  -  250  ms  ISI  loading.  The  original  data  is  plotted  as 
the  IBI  versus  the  time  of  each  beat,  while  the  LPFES  series  is  simply  plot¬ 
ted  at  1  sec  increments.  The  60  sec  time  scale  in  Figure  6.2.1a  shows  that  it 
takes  several  data  points  for  the  interpolation  to  settle  down  which  is  why 
the  first  and  last  10  points  were  discarded  when  computing  the  PSE.  The 
300  sec  time  scale  in  Figure  6.2.1b  shows  that  the  interpolation  algorithm  is 
capable  of  following  the  data  for  the  duration  of  a  trial. 

Figure  6.2.2  is  a  graph  of  the  Blackman  Tukey  PSE  of  the  original  IBI 
series  overlayed  with  the  Blackman  Tukey  PSE  of  the  LPFES  when  the  time 
increment  (dt)  was  set  to  the  mean  IBI  for  the  same  data  as  shown  in  Figure 
6.2.1.  No  normalization  was  performed  on  either  the  x  or  y  axis  since, 
according  to  DeBoer  et  al.  (1984),  these  two  graphs  should  be  equivalent  for 
small  coefficients  of  variation  (er/mean)  of  the  inter-beat  interval.  The  mean 
IBI  for  this  data  was  906  ms  while  the  standard  deviation  was  53  ms.  A 
time  increment  of  906  ms  in  the  French-Holden  algorithm  implies  a  cut  off 
frequency  of  0.55  Hz  (l/(2  x  906  m^)).  Scaling  the  frequency  axis  of  the 
spectrum  of  intervals  PSE  according  to  the  method  outlined  by  DeBoer 
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would  also  place  ttie  maximum  frequency  at  0.55  Hz.  The  frequency  axis  for 
the  spectrum  of  counts  method  (LPFES)  is  a  true  linear  frequency  scale. 
The  frequency  axis  for  the  spectrum  of  intervals  method  is  a  linear  approxi¬ 
mation  with  the  degree  of  linearity  increasing  as  one  gets  closer  to  the  origin 
and  the  coefficient  of  variation  decreases. 

Note  the  emergence  of  a  sharp  peak  in  the  PSE  from  the  LPFES  at 
0.167  Hz  and  0.333  Hz.  While  these  peaks  are  also  present  in  the  spectrum 
of  counts,  they  are  not  as  pronounced.  A  peak  was  consistently  located  at 
0.167  Hz  in  the  PSE  for  all  subject  and  loading  conditions,  though  its  ampli¬ 
tude  varied  greatly  according  to  subject.  It  is  believed  that  this  peak  is  due 
to  entrainment  of  the  heart  rate  by  the  warning  tone  which  was  presented 
every  six  seconds  to  signify  the  start  of  a  new  trial.  Note  that  the  stimuli 
were  not  periodic  since  their  location  was  varied  relative  to  the  warning 
tone.  Mulder  et  al.  (1973)  also  reports  finding  a  peak  at  the  stimulus  fre¬ 
quency  in  the  power  spectrum.  The  peak  at  0.333  Hz  is  believed  to  be  a 
harmonic  of  the  peak  at  0.167  Hz  though  it  could  not  always  be  identified 
due  to  its  overlap  with  the  respiratory  band  in  the  power  spectrum.  The 
mean  inter  respiratory  rate  for  this  subject  and  loading  condition  was  2808 
ms  with  a  standard  deviation  of  394  ms.  This  would  center  the  respiratory 
peak  at  0.356  Hz.  The  degree  of  sinus  arrhythmia,  as  indicated  by  the  per¬ 
cent  of  power  in  the  respiratory  band  of  the  PSE,  varied  greatly  from  sub¬ 
ject  to  subject.  The  current  subject,  Subject  6,  showed  one  of  the  lowest 
degrees  of  sinus  arrhythmia. 

The  superiority  of  the  spectrum  of  counts  (LPFES)  PSE  over  the  spec¬ 
trum  of  intervals  (original  IBI)  PSE  in  identifying  time  synchronized  (as 
opposed  to  beat  synchronized)  events  is  the  primary  reason  that  it  was  used 
in  the  analysis.  One  possible  indicator  of  loading  that  was  considered  was 
the  relationship  between  these  two  PSEs.  If  the  time  increment  of  the 
LPFES  is  set  to  the  mean  IBI,  the  points  on  the  PSE  frequency  axis  should 
roughly  align  with  those  on  the  PSE  from  the  original  IBI  data.  The 
waveform  obtained  by  a  point  by  point  subtraction  of  these  two  PSE  did  not 
appear  to  yield  any  measures  consistent  with  the  loading.  For  the 
remainder  of  the  analysis  of  heart  rate  variability,  the  time  increment  for 
the  French-Holden  algorithm  was  set  to  1  see  (Fmax  0.5  Hz).  Figure 
6.2.3  is  an  overlay  of  the  LPFES  (dt  906  ms)  PSE  shown  in  f  igure  6.2.2 
with  the  PSE  computed  for  a  LPFES  (dt  1.0  sec)  to  dispel  any  doubts 

that  the  peaks  may  somehow  be  dependent  on  t  he  low  pass  filter  in  t  ho 
interpolation  algorithm. 


The  mean  of  a  series  is  usually  removed  before  computing  the  power 
spectrum  in  order  to  remove  the  large  DC  component  from  the  spectrum.  If 
it  is  believed  that  the  series  has  a  trend  which  is  not  relevant  to  the  process 
being  measured,  then  the  trend  can  also  be  removed  before  computing  the 
power  spectrum.  This  will  remove  the  power  at  DC  and  the  extreme  low  fre¬ 
quency  end  of  the  spectrum.  Some  of  the  data  records  did  indicate  a  trend, 
though  it  did  appeared  to  vary  by  subject  and  loading  condition.  A  linear 
regression  was  performed  on  all  of  the  cases,  but  the  slope  did  not  appear  to 
be  consistent  with  the  workload.  Working  witli  the  differences  of  the  points 
is  another  technique  that  is  used  when  it  is  believed  that  extraneous  trends 
are  present.  This  is  the  dual  of  taking  the  derivative  of  a  continuous 
waveform,  which  would  multiply  its  resultant  PSE  by  a  .J  window.  The  net 
effect  is  to  reduce  the  power  in  the  low  frequency  bands  while  emphasizing 
the  power  in  the  high  frequency  bands. 

In  order  to  be  complete,  PSEs  were  computed  using  each  of  the  above 
prefiltering  operations:  demeaning  -  removal  of  the  mean,  detrended  -  remo¬ 
val  of  the  best  fit  line,  and  differenced  -  using  the  difference  between  suces- 
sive  1BI.  A  PSE  was  also  computed  using  the  heart  rate,  which  is  the  output 
of  the  French-Holden  algorithm,  rather  than  the  IBI.  The  mean  rate  was 
subtracted  out  before  the  PSE  was  computed.  The  four  types  of  PSE  were 
all  computed  using  the  LPFES  with  a  time  increment  of  1  sec.  Note  that 
this  time  increment  was  greater  than  the  mean  IBI  for  Subjects  4,  5  &  6,  less 
than  the  mean  IBI  for  Subjects  8  &  9  and  about  the  same  as  the  mean  IBI 
for  Subject  7.  Figure  6.2.4  shows  four  PSEs  computed  in  the  above  fashion 
for  Subject  8  under  both  the  4  Choice  -  62.5  ms  ISI  (high)  loading  and  the  2 
Choice  -  250  ms  ISI  (low)  loading.  The  large  scale  on  the  power  axis  is  due 
to  working  in  a  millisecond  time  scale  for  the  inter-beat  intervals.  Note 
than  when  using  the  French-IIolden  algorithm  it  is  possible  to  get  aliasing  in 
the  PSE.  The  PSE  from  the  differenced  beats  emphasize  the  aliasing  due  to 
its  nonlinear  weighting  of  the  PSE.  It  was  judged  that  the  aliasing  was  not 
significant  in  the  bands  of  interest  and  the  aliasing  for  the  low  load  condi¬ 
tion  shown  in  Figure  6.2.4  was  typical.  The  degree  of  aliasing  indicated  by 
the  high  load  condition  in  Figure  6.2.4  was  fairly  rare. 

Note  the  large  reduction  in  power  in  the  high  loading  case  compared  to 
the  lower  loading  case.  Six  loading  conditions  for  six  subjects  and  four  types 
of  prefiltcring  yields  144  graphs  of  the  type  shown  in  Figure  6.2.4.  Then'  are 
too  many  conditions  to  compare  manually.  It  was  decided  to  analyze  this 
data  by  comparing  the  power  in  four  frequency  bands.  The  bands  were 
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chosen  to  correspond  with  the  stimulus  frequency  as  well  as  the  three  main 
physiological  factors  known  to  elfect  heart  rate  (Sayers  1973);  thermal  regu¬ 
lation,  blood  pressure  control,  and  respiration. 

The  location  of  the  stimulus  band  was  easy  to  determine  since  it  is  cen¬ 
tered  at  l/(stimulus  period)  or  0.167  Hz.  It  was  decided  to  sum  up  the 
power  in  the  band  from  0.1484  to  0.1796  Hz  (5  points  wide)  to  represent  the 
contribution  of  the  stimulus  to  the  PSE.  It  was  not  always  clear  from  the 
PSE  where  the  respiration  peak  was  located  so  it  was  decided  to  use  the 
respiration  measure  itself  to  define  the  respiration  band.  For  the  low  load 
condition  shown  in  Figure  6.2.4  the  mean  respiratory  interval  was  3553  ms 
(~  0.28  Hz)  with  a  standard  deviation  of  537  ms.  For  the  high  load  condi¬ 
tion  it  was  3210  ms  (~  0.31  Hz)  with  a  standard  deviation  of  298  ms.  The 
respiratory  band  was  defined  to  start  at  l/(mean+oj  and  end  at 
l/(mean— rr).  This  corresponds  to  bands  of  0.245  -  0.332  Hz  and  0.285  - 
0.343  Hz  for  the  low  and  high  loading  cases  respectively.  Note  that  this 
aligns  with  the  peak  evident  in  the  high  load  case  and  allows  us  to  make  an 
objective  decision  in  the  low  load  case  where  a  single  peak  cannot  be  iso¬ 
lated. 

The  choices  for  the  thermal  and  blood  pressure  band  are  not  as  straight 
forward  as  for  the  stimulus  and  respiration  band.  We  do  not  have  indepen¬ 
dent  measures  of  body  temperature  or  blood  pressure.  Sayers  (1973)  has 
correlated  the  heart  rate  with  independent  measures  of  body  temperature 
and  blood  pressure  He  reports  the  thermal  band  to  be  centered  at  0.025  Hz 
and  the  blood  pressure  band  to  vary  from  0.07  to  0.10  Hz  (see  Figure  3.1.1). 
The  frequency  increment,  when  using  the  aforementioned  procedure  for 
computing  the  PSE,  is  0.00781  Hz.  The  power  at  a  particular  frequency  is 
spread  over  several  adjacent  points  in  the  frequency  domain,  commonly 
called  leakage  effect,  when  the  PSE  is  computed  from  a  finite  data  length. 
Based  on  the  above  facts  and  from  viewing  PSE  across  the  subject  and  load¬ 
ing  cases  for  peak  locations,  it  was  decided  to  define  the  thermal  band  as 
DC  to  0.0391  Hz  and  the  blood  pressure  band  as  0.0547  to  0.1250  Hz.  The 
blood  pressure  activity  is  not  clearly  indentifiable  in  Figure  6.2.4,  but  in 
other  subjects  this  band  was  found  to  align  with  significant  activity.  The 
thermal  band  chosen  was  in  good  correlation  with  the  activity  shown  in  Fig¬ 
ure  6.2.4. 

The  bar  graph  shown  in  Figure  6.2.5  shows  the  total  power  (DC  to  0.5 
Hz)  for  all  subjects,  loading  and  prefiltering  conditions.  Figures  6.2.6,  6.2.7 
and  6.2.8,  6.2.9  are  bar  graphs  which  show  the  power  in  the  thermal,  blood 
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pressure,  stimulus,  and  respiration  bands  respectively.  The  axis  for  the 
power  (ordinate)  in  Figures  6.2.5  to  6.2.9  is  an  offset  dB  scale  that  ranges 
from  0  to  30  dB.  The  offsets  used  between  bar  graphs  were  not  constant  in 
order  to  be  able  to  show  the  bars  at  a  reasonable  height.  The  offset  within 
an  individual  bar  graph  was  constant  and  is  indicated  on  the  axis. 

Before  discussing  these  graphs  which  summarize  the  spectral  nature  of 
the  heart  beat,  the  behavioral  measures  of  response  (Table  5.2.2)  should  be 
reviewed.  The  most  significant  changes  in  behavioral  performance  occured 
for  three  of  the  binary  comparisons  between  loading  (columns  3,  5  &  6  in 
Table  5.2.2).  The  three  binary  comparisons  were: 

4  Choice  -  62.5  ms  vs.  4  Choice  -  250  ms(ISl  Loading) 

4  Choice  -  62.5  ms  vs.  2  Choice  -  62.5  ms(Choice  Loading) 

4  Choice  -  62.5  ms  vs.  2  Choice  -  250  ms(High-Low  Loading) 

Note  that  all  of  the  comparisons  involved  the  highest  loading  level  of  4 
Choice  -  62.5  ms.  It  may  be  that  comparisons  between  the  other  levels  of 
loading  are  not  as  significant  because  the  subject  was  not  sufficiently  loaded 
to  produce  a  marked  behavioral  or  physiological  difference.  We  will  refer  to 
these  three  companions  of  loading  by  the  above  names  which  are  given  in 
parentheses.  In  looking  at  the  bar  graphs  (Figures  6.2.5  to  6.2.9)  this 
corresponds  to  the  sixth  bar  in  each  group  being  lower  than  the  third, 
fourth,  and  fifth  bar  if  the  power  in  each  band  decreased  with  loading.  ISI 
loading  is  indicated  by  comparing  the  fourth  bar  to  the  sixth  bar.  Choice 
loading  is  indicated  by  comparing  the  fifth  bar  to  the  sixth  bar.  High-low 
loading  is  indicated  by  comparing  the  third  bar  to  the  sixth  bar. 

In  examining  total  power  shown  in  Figure  6.2.5,  the  demeaned  and 
detrended  method  yield  the  correct  relationships  for  all  three  comparisons 
for  all  subjects  with  the  exception  of  Subject  6.  The  heart  rate  yields  simi¬ 
lar  results  except  that  it  also  fails  for  the  high-low  comparison  for  Subject  5. 
The  differences  method  fails  all  three  comparisons  for  three  of  the  subjects 
and  passes  all  three  comparisons  for  the  other  three  subjects.  This  can  be 
attributed  to  its  extremely  nonlinear  windowing  of  the  PSE.  Note  that  all 
four  prefiltering  methods  passed  all  three  comparisons  for  the  mean  of  the 
six  subjects.  Extreme  caution  should  be  exercised  when  simply  working  with 
the  mean  since  the  mean  of  the  difference  method  passed  all  three  tests 
while  half  of  the  individual  subjects  failed. 

In  comparing  the  four  filtering  methods  when  looking  at  the  individual 
frequency  bands,  all  methods  gave  roughly  similar  performance.  The 
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problem  with  the  nonlinear  weighting  of  the  frequency  axis  by  the  differences 
method  is  not  very  severe  when  the  power  is  computed  from  only  a  few  adja¬ 
cent  points.  There  was  no  single  band  for  any  of  the  prefiltering  methods 
that  could  correctly  classify  any  of  the  loading  comparisons  for  all  six  sub¬ 
jects.  If  the  acceptance  criterion  is  dropped  to  5  out  of  6  subjects,  then  the 
thermal  band  is  the  only  band  that  can  correctly  classify  some  of  the  loading 
comparisons.  The  demeaned  thermal  band  can  distinguish  the  choice  com¬ 
parison,  while  both  the  detrended  and  differenced  thermal  band  can  distin¬ 
guish  the  high-low  and  choice  comparison. 

The  relative  heights  of  the  bars  in  the  graph  of  the  PSE  for  the 
demeaned  prefiltering  condition,  shown  in  Figure  6.2.5,  should  correspond  to 
the  relative  heights  of  the  bars  of  the  graph  of  the  standard  deviation, 
shown  in  Figure  5.3.2.  The  only  notable  exception  to  this  is  Subject  7  under 
the  \  choice  -  62.5  msec  loading  condition.  The  IBI  plotted  in  Figure  5.3.5  is 
from  this  subject  and  loading  case.  The  initial  few  seconds  at  the  start  of 
the  block  had  a  IBI  which  was  significantly  lower  than  mean  (by  2-  300  ms). 
In  computing  the  PSE,  these  points  were  truncated.  They  were  not  trun¬ 
cated  in  the  computation  of  the  standard  deviation  which  explains  the 
inflated  estimate.  If  this  error  is  compensated  for,  then  the  standard  devia¬ 
tion  will  have  same  classification  performance  as  the  total  power.  The  stan¬ 
dard  deviation  estimate  computed  by  the  mean  square  successive  difference 
(MSSD),  as  shown  in  Figure  5.3.3,  should  not  be  subject  to  this  bias.  The 
MSSD  estimate  of  the  standard  deviation  could  only  classify  four  out  of  the 
six  subjects  for  the  choice  comparison.  Its  performance  for  the  other  com¬ 
parisons  was  lower  so  that  its  use  as  a  classifier  is  dubious  at  best. 


A  relative  measure  of  the  power  in  a  band  may  indicate  effects  even 
though  the  absolute  measure  of  the  power  in  the  band  might  not.  Figure 

6.2.10  is  a  bar  graph  which  indicates  the  percent  of  the  total  power  lying  in 
each  of  the  four  bands  using  the  demeaned  method  of  prefiltcring.  Six 
binary  comparisons  can  be  made  when  using  four  distinct  bands.  Figure 

6.2.11  shows  the  relative  power  between  the  bands  for  the  six  binary  com¬ 
parisons.  Figure  6.2.11  uses  a  dB  scale  so  that  a  negative  number  indicates 
that  the  first  band  had  less  power  than  the  second  band.  Unfortunately, 
none  of  the  graphs  in  Figures  6.2.10  &  6.2.11  can  distinguish  between  any  of 
the  loading  comparisons  for  all  six  subjects.  At  the  five  out  of  six  subject 
acceptance  level,  there  is  a  slightly  greater  percentage  of  power  in  the 
stimulus  band  for  the  high-low  comparison  except  for  Subject  7.  This  also 
shows  up  as  a  decrease  in  the  ratio  of  the  thermal  band  power  to  the 
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stimulus  band  power  for  the  high-low  comparison  except  for  Subject  5.  The 
ratio  of  blood  pressure  power  to  stimulus  band  power  also  decreases  for  the 
choice  comparison  except  for  Subject  7.  The  best  classifier  appears  to  be  the 
total  power  in  the  demeaned  or  detrended  PSE  of  the  LPFKS. 

6.3.  P300  as  a  Measure  of  Loading 

The  P300  or  positive  peak  occurring  in  the  event  related  potential 
(ERP)  waveforms  at  ~  300  ms  after  stimulation  has  been  extensively  studied 
as  an  indicator  of  cognitive  activity.  Unfortunately,  the  procedures  to  iden¬ 
tify  and  measure  P300  have  been  far  from  consistent  and  even  vary  within 
the  same  research  group.  The  use  of  P300  as  a  measure  involves  three  steps; 
prefiltering,  identification  of  the  peak  and  computation  of  the  measure  from 
the  identified  peak. 

The  ERP  is  usually  low  pass  filtered  to  reduce  any  alpha  activity  that 
might  hinder  the  identification  of  P300.  Magliero  and  Donchin  (1984)  used 
3.5  Hz  (-3  dB)  low  pass  prefiltering,  while  Kutas  and  Donchin  et  al.  (1977) 
used  6.29  Hz  low  pass  prefiltering.  As  part  of  the  signal  processing  for  this 
research,  a  sharp  cutoff  low  pass  filter  was  designed.  The  Remez  Exchange 
Algorithm  (McClellan  et  al.  1979)  provides  the  design  of  a  low  pass  filter  by 
specifying  the  passband  and  stopband.  This  is  a  safer  approach  than  simply 
specifying  the  3  dB  point  of  a  filter  since  it  eliminates  the  possibility  that 
the  transition  region  of  the  filter  may  extend  into  the  desired  pass  and  stop- 
bands.  It  was  decided  that  the  prefiltcring  should  consist  of  a  low  pass  filter 
with  a  passband  ending  at  3.5  Hz  and  a  stopband  starting  at  8  Hz  to  block 
the  alpha  activity.  The  resultant  filter  provided  20  dB  of  attenuation  in  the 
stop  band  with  a  0.15  dB  ripple  in  the  passband  and  a  -3  dB  point  at  5.6 
Hz.  The  result  of  this  filtering  on  the  average  response  (P3  electrode)  for 
several  subjects  is  shown  in  Figures  6.3.1,  6.3.2  and  6.3.3  for  the  single 
stimulation,  250  rns  IS  I ,  and  62.5  ms  ISI  loading  conditions,  respectively. 
Both  2  choice  and  4  choice  conditions  are  shown  in  the  figures. 

Peak  identification  has  traditionally  consisted  of  finding  the  largest 
amplitude  within  a  specified  interval.  The  interval  used,  however,  varies 
with  the  researcher.  Kavis  and  Donchin  (1983)  used  a  300  to  750  ms  inter¬ 
val,  while  McCarthy  and  Donchin  (1981)  used  a  200  to  1500  ms  interval. 
Virtually  every  paper  that  deals  with  P300  uses  a  different  interval.  The 
current  research  allows  the  recorded  data  to  determine  the  interval  based  on 
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a  statistical  test  rather  than  experimenter  bias.  The  mechanisms  for  this 
procedure  is  described  in  the  subsequent  section  entitled  "Latency  Corrected 
Average". 

The  purpose  in  locating  the  peaks  is  so  that  some  metric  can  be  calcu¬ 
lated  which  may  prove  useful  in  determining  levels  of  cognitive  activity. 
Some  researchers  avoid  the  issue  of  trying  to  locate  a  specific  peak  in  each 
trial  by  using  a  discriminant  function  (Squires  &  Donchin,  1976;  Duncan- 
Johnson  &  Donchin;  1981).  This  discriminant  function  uses  a  sum  of 
weighted  values  applied  to  the  amplitude  at  several  latencies  and  across 
several  electrodes.  These  weightings  were  determined  by  comparing 
classification  accuracies  using  a  Step-wise  Discriminant  Analysis  (SWDA) 
procedure  for  cognitive  tasks  across  a  large  subject  base.  Other  measures 
that  have  been  computed  include  the  peak  latency  and  amplitude  or  area 
computations.  Latency  is  usually  computed  from  the  peak  in  the  average 
(McCarthy  &  Donchin,  1981;  Magliero  &  Donchin  1984)  and  has  been 
reported  to  vary  from  210  ms  to  550  ms  after  stimulation  (John,  1977). 
Amplitude  measures  can  be  taken  from  the  amplitude  of  the  largest  peak  in 
the  conventional  average,  the  amplitude  of  the  largest  peak  in  the  Woody 
Correlation  Average  (Kutas  &  Donchin  et  ah,  1977)  or  the  average  from 
aligning  all  of  the  largest  peaks  found  in  the  individual  trials  (Duncan- 
Johnson  &  Donchin,  1981).  The  area  of  an  interval  around  300  ms  has  also 
been  used  as  a  measure  since  it  minimizes  the  effect  of  latency  jitter  on  the 
peaks  (Israel  et  ah,  1980).  Israel  also  computed  the  latency  jitter  as  a  meas¬ 
ure  derived  from  P300.  The  present  research  uses  the  latency  corrected 
average  to  locate  and  compute  measures  for  P300. 

6.4.  The  Latency  Corrected  Average 

The  (LCA)  latency  corrected  average  (McGillem  &  Aunon,  1977;  Aunon 
&  McGillem,  1979;  Begleiter  et  ah,  1983)  was  chosen  as  the  method  of  identi¬ 
fying  cognitive  measures  from  the  ERP.  For  the  purposes  of  discussing  the 
LCA,  the  term  "LCA  peak"  shall  refer  to  the  average  peak  computed  by  the 
program  and  the  term  "trial  peak"  shall  refer  to  a  peak  from  an  individual 
trial  that  was  used  in  computing  the  average  peak. 

The  LCA  procedure  was  selected  for  several  reasons.  Trial  peaks  are 
selected  only  from  statistically  significant  regions  which  are  determined  indi¬ 
vidually  for  each  block  (50  trials).  An  LCA  peak  is  computed  for  each  valid 
region,  rather  than  just  using  a  single  region  to  represent  the  entire 
response.  The  LCA  allows  for  the  possibility  of  a  trial  to  have  a  null 


response;  i.e.  it  is  not  required  that  every  trial  have  a  valid  peak  in  each 
region.  Trial  peak  statistics  (percent  occurrence,  mean  latency,  latency  vari¬ 
ance  and  mean  amplitude)  are  computed  only  from  valid  peaks  within  the 
region.  The  latency  of  an  LCA  peak  is  determined  by  averaging  the  trial 
peak  latencies,  rather  than  averaging  the  peaks  themselves,  and  then 
searching  for  the  peak  in  the  average. 

The  LCA  procedure  is  implemented  in  the  "MLCAY"  program  and 
starts  by  computing  and  applying  a  minimum  mean  square  error  filter 
(MMSE)  filter  to  the  normalized  and  detrended  data.  All  peaks  within  the 
individual  trials  are  detected  by  crosscorrelating  the  trials  with  a  template 
having  the  generalized  shape  of  a  peak  in  the  ERP.  Histograms  of  both  the 
positive  and  negative  peaks  are  computed  at  each  latency  from  the  ensemble 
of  trials.  A  sliding  window  is  then  applied  to  the  histograms  and  a  non- 
parametric  sign  test  (Sachs  sec  4.2.4,  1984)  is  performed  to  determine  if  the 
data  in  the  window  is  from  a  positive,  negative  or  zero  median  process  at  a 
specified  confidence  level  (95%).  If  the  sign  test  indicates  a  positive  median 
for  that  window  then  the  region  is  considered  to  be  one  of  significant  posi¬ 
tive  peaks.  The  record  is  then  divided  up  into  regions  of  positive,  negative 
and  no  peaks.  For  each  positive  region,  the  largest  positive  peaks  from  each 
trial  are  aligned  at  the  mean  trial  peak  latency  and  averaged  over  an  eleven 
sample  point  width.  Negative  peaks  are  processed  in  a  similar  fashion. 

The  output  of  the  "MLCAY"  program  is  a  listing  of  the  locations  of  the 
positive  and  negative  peaks  (up  to  a  maximum  of  10  each)  and  the  average 
(11  sample  points  wide)  of  the  trial  peaks  used  to  compute  that  LCA  peak. 
Statistics  for  the  LCA  peak  are  also  computed  (percent  occurrence  of  the 
trial  peaks  in  the  region,  mean  amplitude,  latency  mean  and  variance).  A 
listing  of  the  location  and  amplitude  of  the  trial  peaks  used  for  each  LCA 
peak  is  also  output.  Intermediate  results,  such  as  histograms,  MMSE  filter, 
etc  may  also  be  output  to  check  program  operation. 

Figure  6.4.1  shows  a  plot  of  the  LCA  peaks  computed  from  the  same 
data  shown  as  an  average  in  Figure  6.3.3.  Note  that  the  LCA  identifies 
some  peaks  that  do  not  show  up  in  the  average  due  to  smearing  caused  by 
latency  jitter.  The  P300  peak  has  been  clearly  identified  for  all  three  sub¬ 
jects  in  the  400  to  500  ms  range.  The  4  choice  -  62.5  ms  ISI  P300  response 
clearly  occurs  later  than  the  2  choice  -  62.5  ms  ISI  P300  which  is  in  agree¬ 
ment  with  previous .  findings  that  P300  latency  increases  with  loading 
(Review  in  Pritchard,  1981).  The  negative  peak  in  the  100-200  msec  range 
(commonly  called  N200)  also  appeared  to  be  slightly  delayed  in  latency  for 
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the  higher  loading.  The  comparisons  indicated  on  this  graph  are  for  the 
choice  loading  comparison  that  was  discussed  in  section  6.2. 

Figures  6.4.2,  6.4.3  &  6.4.5  show  the  results  of  the  LCA  for  the  P4  elec¬ 
trode  for  Subjects  6,  7,  &  8  respectively.  The  highest  loading  case  in  these 
figures  are  the  dashed  line  in  the  bottom  graph.  The  choice  comparison 
compares  the  dashed  line  in  the  bottom  graph  to  the  solid  line  in  the  bottom 
graph.  The  ISI  comparison  compares  the  dashed  line  in  the  bottom  graph  to 
the  dashed  line  in  the  middle  graph.  The  high-low  comparison  compares  the 
dashed  line  in  the  bottom  graph  to  the  solid  line  in  the  middle  graph. 

The  shift  in  the  P300  for  the  choice  comparison,  as  shown  in  Figure 
6.4.1  for  the  P3  electrode,  is  also  evident  for  all  of  the  subjects  for  the  P4 
electrode.  It  is  hard  to  judge  the  ISI  and  high-low  comparisons  because  the 
output  of  the  LCA  for  the  250  ms  ISI  is  somewhat  erratic.  It  does  not 
appear  to  find  the  N200  or  the  P300  for  Subject  7  under  the  4  Choice  -  250 
ms  ISI  loading  condition  as  shown  in  Figure  6.4.3.  For  Subject  6,  the  LCA 
appears  to  find  a  N200  and  P300  for  each  of  the  stimuli  that  were  presented 
250  ms  apart  (Figure  6.4.2  middle  graph).  It  may  be  that  in  some  of  the 
subjects  the  stimuli  presented  250  ms  apart  are  be  fused  and  treated  as  a 
single  stimulation.  The  P300  peak  does  appear  to  shift  for  both  the  ISI  and 
high-low  loading  comparisons  for  Subjects  7  &  8.  If  we  only  consider  the 
first  "P300"  peak  of  Subject  6,  then  these  comparisons  also  hold  true. 

Figure  6.4.5  shows  a  scatter  plot  of  reaction  time  to  the  first  stimulus 
(RTl)  and  the  P300  latency.  Reaction  time  has  been  shown  to  correlate 
with  P300  latency  under  certain  conditions  (McCarthy  &  Donchin,  1981) 
Kutas  and  Donchin  (et  al.,  1977)  have  shown  that  there  is  a  higher  degree  of 
correlation  in  tasks  that  stress  accuracy  considerations  over  speed  considera¬ 
tions.  In  the  present  research,  the  subject  was  instructed  that  speed  was 
more  important  than  accuracy  which  may  account  for  the  low  correlation 
between  RTl  and  P300  latency. 

6.5.  Conclusions 

The  task  description  based  on  the  hybrid  capacity  model  was  initially 
designed  to  produce  6  distinct  levels  of  loading.  The  behavioral  variables 
were  not  able  to  distinguish  between  the  6  distinct  levels.  This  means  that 
we  did  not  have  a  converging  criterion  to  justify  that  we  had  indeed 
changed  the  loading  significantly  in  all  6  cases.  The  heart  rate  information 
showed  a  decrease  in  the  IBI  variabilty  when  comparing  the  highest  level  of 
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loading  with  the  three  other  levels  that  were  found  to  have  a  significant 
difference  as  measured  by  the  behavioral  response.  The  IB1  variability  was 
not  able  to  distinguish  between  the  levels  in  cases  where  the  behavioral 
response  could  not  distinguish  between  the  levels.  The  EEG  also  indicated 
an  increase  in  P300  latency  which  correlated  with  the  increase  in  load  and 
decrease  in  variability.  Unfortunately  no  useful  information  was  obtained 
from  the  MEG  data  due  to  head  movement.  These  results  indicate  that 
there  are  significant  changes  in  the  physiological  variables  and  the  impor¬ 
tance  of  providing  a  converging  criterion. 
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APPENDIX  B 

"MUTUAL"  -  SQUID  Calibration  Program 


"MUTUAL"  -  SQUID  Calibration  Program 


The  mutual  program  is  a  short  program  that  computes  the  mutual 
inductance  and  coupling  between  the  SQUID  detector  coils  and  a  test  loop 
that  is  coaxial  to  them.  It  was  originally  written  for  and  tested  on  a  BTI 
(formarlly  S.H.E.  Corp.)  single  channel  DC  SQUID  with  a  second  order  gra- 
diometer  probe.  Tables  from  Grover  (1946,  1973)  are  used  by  the  program 
to  compute  the  mutual  inductances.  Once  the  mutual  inductance  has  been 
computed,  the  coupling,  or  flux  induced  in  the  detector  for  a  given  current  in 
the  test  loop  is  computed. 

The  program  was  written  in  the  C  programming  language  by  Michael 
Plonski  of  Purdue  University  while  on  temporary  assignment  to  the  Naval 
Personnel  Research  Center  (NPRDC)  in  San  Diego.  All  measurements  are 
relative  to  the  bottom  of  the  coil  form  on  which  the  detector  coils  are 
wound.  The  bottom  of  the  dewar  tail  is  typically  0.394  inches  from  the  bot¬ 
tom  of  the  coil  form.  The  detector  coil  dimensions  are  stored  in  define  state¬ 
ments  at  the  beginning  of  the  program.  These  measurements  should  be 
verified  with  those  supplied  by  the  manufacturer.  This  program  has  also 
been  used  at  the  Aerospace  Medical  Research  Lab  at  Wright  Patterson  Air 
Force  Base  where  these  values  had  to  be  changed  to  meet  the  specifications 
of  their  particular  SQUID. 

The  source  code  for  the  program  used  at  NPRDC  is  included  here.  A 
1.2  inch  radius  test  loop  was  made  out  of  wire  wrap  wire  and  was  used  to 
test  the  program.  This  size  loop  fit  snuggly  around  the  dewar  tail  and  the 
remainder  of  the  wire  was  twisted  tightly  with  a  drill  to  prevent  any  stray 
coupling  from  the  wire  to  the  detector  loop.  The  test  loop  was  moved  from 
the  bottom  of  the  dewar  tail  (0.394  inches  below  the  coil  form  bottom)  to 
two  inches  above  the  dewar  tail  (1.6  inches  above  the  coil  form)  and  meas¬ 
urements  where  taken  at  0.1  inch  increments. 

Figure  Bl  shows  the  comparision  between  the  coupling  predicted  by  the 
program  and  coupling  measured  by  the  SQUID.  The  measured  coupling  was 
simply  the  number  of  femtoteslas  (fT)  measured  by  the  SQUID  using  the 
manufacturer  supplied  calibration,  divided  by  the  current  in  the  loop  as 
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measured  with  a  oscilloscope  and  series  resistor.  All  measurements  were 
made  using  a  10  Hz  sine  wave  and  peak  to  peak  values  were  compared.  The 
predicted  coupling  was  computed  over  a  larger  range  from  0.5  inches  below 
the  coil  form  bottom  to  2.5  inches  above  to  demonstrate  the  symmetry  as 
the  loop  passes  over  the  various  detector  coils.  The  predicted  and  measured 
coupling  closely  overlay  in  the  -0.4  to  1.6  inch  range  over  which  the  meas¬ 
ured  values  were  taken.  The  peak  coupling  occurs  when  the  test  loop  is  cen¬ 
tered  over  the  four  loops  that  comprise  the  middle  set  of  the  detector  loops. 
Two  nulls  are  present  when  the  test  loop  is  placed  between  the  midlle  and 
lower  set  of  coils  and  again  when  it  is  placed  between  the  middle  and  upper 
set  of  coils.  This  is  because  the  middle  set  of  loops  is  wound  in  the  opposite 
direction  of  the  other  loops. 

A  test  fixture  can  be  made  by  wrapping  a  coil  around  a  wooden  form 
that  could  be  placed  over  the  dewar  tad.  Care  should  be  taken  when  desig- 
ing  a  test  fixture  so  that  the  test  loop  aligns  over  an  portion  of  the  curve 
with  a  low  slope  (preferably  centered  over  the  lower  set  of  coils)  in  order  to 
minimize  errors  due  to  placement.  A  second  loop  could  also  be  placed  over 
one  of  the  nulls  to  check  that  the  probe  aligment  has  not  shifted  over  time. 
The  test  fixture  would  also  allow  quick  operational  checks  and  ther  ability 
to  compute  the  transfer  function  for  the  entire  measurement  system  rather 
than  just  the  amplifier  and  filter  stage. 


"Inductance  Calulations"  1946,  1973,  F.  Grover,  published  by  the  Instrument 
Society  of  America  uner  an  arrangement  with  Dover  Publications,  Inc. 


Figure  Bl  Predicted  vs  Measured  Coupling 
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PROGRAM  MUTUAL 

t^tttt************ ********************** ********************************* 


THIS  PROGRAM  WILL  COMPUTE  THE  EFFECTIVE  MUTUAL  INDUCTANCES 
BETWEEN  A  SINGLE  LOOP  AND  THE  SECOND  ORDER  G RADIOMETER 
DETECTION  COILS  USING  EQS  77,78  &  TABLE  13  OF  "INDUCTANCE 
CALCULATIONS"  BY  GROVER,  DOVER  PUBLICATIONS.  ALL  VALUES  ARE 
IN  INCHES  FROM  BOTTOM  OF  QUARTZ  COIL  FORM.  APPROXIMATE 
DISTANCE  FROM  BOTTOM  OF  COIL  FORM  TO  BOTTOM  OF  DEWAR  TAIL 
IS  .394  In.  PROGRAM  WAS  CHECKED  USING  DBX  (C-debugger)  TO 
ALTER  THE  VALUES  TO  THOSE  USED  IN  EXAMPLE  25  OF  GROVER. 

Written  by.  Mike  Plonski  NPRDC,  6/85 

*************************************************************************/ 


#include  <stdio.h> 

#include  <math.h> 

/*  FUNCTIONS  7 

extern  float  fvai(); 

#define  sq(x)  ((x)*(x))  /*  square  function  */ 

# define  abs(x)  ((x)  >  0?  (x):  (-(x)))  /*  absolute  value  func  */ 

#define  intocm(x)  ((x)*2.54)  /*  inch  to  cm  convert  */ 

I*  DETECTOR  COIL  PARAMETERS  &  TABLE  13  */ 

float  xlQ  =  {.027,  .057},  /*  location  of  lower  coils  */ 

xm[]  =  {1.247,  1.277,  1.307,  1.337},  /*  middle  coils  */ 

xu[]  =  {2.527,  2.557},  /*  upper  coils  */ 

rd  =  .4885,  /*  radius  of  det  coils  */ 

kinc  =  .01,  /*  increments  in  k  for  table  13  input  */ 

f[]  =  {  /*  table  13  from  grover  */ 

1.0,  .021474,  .017315,  .014937,  .013284, 

.012028,  .011017,  .010179,  .009464,  .008843, 

.008297,  .007810,  .007371,  .006974,  .006611, 

.008278,  .005970,  .005685,  .005420,  .005173, 

.004941,  .004723,  .004518,  .004325,  .004142, 

.003969,  .003805,  .003649,  .003500,  .003359, 

.003224,  .003095,  .002971,  .002853,  .002740, 

.0026317,  .0025276,  .0024276,  .0023315,  .0022391, 

.0021502,  .0020646,  .0019821,  .0019026,  .0018259, 

.0017519,  .0016805,  .0016116,  .0015451,  .0014808, 

.0014188,  .0013585,  .0013004,  .0012443,  .0011900, 

.0011374,  .0010865,  .0010373,  .0009897,  .0009436, 

.0008990,  .0008558,  .0008141,  .0007736,  .0007345, 

.0006968,  .0006600,  .0006246,  .0005903,  .0005571, 

.0005251,  .0004941,  .0004642,  .0004353,  .0004074, 

.0003805,  .0003545,  .0003295,  .0003054,  .0002823, 

.00025998,  .00023859,  .00021806,  .00019840,  .00017959, 


main() 


.00016162,  .00014450,  .00012821,  .00011276,  .00009815, 
.00008438,  .00007146,  .00005940,  .00004824,  .00003798, 
.00002866,  .00002035,  .00001312,  .00000708,  .00000249, 
0.0}; 


/*  INTERNAL  VARIABLE  DECLARATIONS 
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float 


rc, 

sc, 

ml, 

mm, 
mu, 
mt, 
coup 
pi= 


/*  location  of  cal  coil  relative  to  form  bottom  * / 
/*  radius  of  calibration  coil  */ 
/*  scale  factor(cm)=sqrt(rc*rd)  */ 
/*  mutual  inductane  with  lower  coil  *  / 
/*  middle  coil  * / 

/*  upper  coil  *  / 

/*  effective  mutual  ind  =  mu  -  mm  +  ml 
,  /*  coupling  coeff.  in  ft/micro  amp  */ 

3.14159265;  /*  standard  definition  */ 


/*  GET  CALIBRATION  COIL  PARAMETERS 


V 


printf("Note  Coil  Form  Bottom  is  .394  in  from  Bottom  of  Dewar  Tail.0); 
printf("Distance  (inches)  from  Calibration  Coil  to  Coil  Form  Bottom  ?  "); 
scanf("%f",&xc); 

printf("Radius  (inches)  of  calibration  coil  ?  “); 

scanf("%f",&rc); 

sc  =  intocm(sqrt(rc*rd)); 

/*  COMPUTE  MUTUAL  INDUCTANCES  •/ 

ml  =  sc*(fval(rc,abs(xc-xl[0]))  +  fval(rc,abs(xc-xl[l])));  /*  eq77  * / 

mu  =  sc*(fval(rc,abs(xc-xu[0]))  +  fval(rc,abs(xc-xu[l]))); 
mm  =  sc*(fval(rc,abs(xc-xm[0j))  +  fval(rc,abs(xc-xm[l]))); 
mm  +=  sc*(fval(rc,abs(xc-xm(2]))  +  fval(rc,abs(xc-xm[3]))); 
mt  =  abs(mu  -  mm  +  ml); 

/*  PRINT  OUT  RESULTS  */ 

printf("  MUTUAL  INDUCTANCES  (micro  henrys)  0); 
printfj"  Lower  Coil  %(  0,ml); 

printfj"  Middle  Coil  %f  0,mm); 

printf("  Upper  Coil  %(  0,mu); 

printfj"  Effective  Total  %{  0,,nt); 

/*  COMPUTE  AND  OUTPUT  COUPLING  COEFFICIENT 

/*  coupling  -■  mutual  inductance  /  effective  area,  for  compatibilty  with  * / 
/*  S.H.E.  field  sensitivity  value  use  combined  area  of  bottom  two  turns  */ 


coup  =  10000000.0  *  mt  /  (2*pi*sq(intocm(rd))); 

printff'  Coupling  Coefficent  —  %8.1f  femto-tesla  /  micro-amp  0,roup); 


exit(O); 


& 

fv 

tv 
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float  fval(rc,d)  /*  return  interpolated  value  from  table  13  */ 

float  rc,  /*  radius  of  calibration  coil  */ 

d;  /*  distance  between  cal  coil  &  current  det  coil  */ 

< 

float  k,  /*  input  parameter  for  table  13  */ 

kdelta,  /*  norm  diff  between  true  k  &  table  entry  used  */ 

dl,  /*  first  difference  in  table 

d2,  /*  second  difference  in  table 

value;  /*  final  value  to  return  */ 

int  ik;  /*  floor  index  for  k  in  table 
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k  =  (sq(rd-rc)  +  sq(d))/(sq(rd+rc)  +  sq(d));  /*  eq.  78  */ 

ik  =  k/kinc; 

kdelta  =  (k  -  ik*kinc)/kinc; 

if(k  <  .1  ;;  k  >  .9)  /**  WARNING 
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{ 


} 


printf(”warning:  k’2  =  %f  is  outside  desired  range  for  this  tableO, k); 
printf("  see  Grover  (ch.  11)  for  alternate  tables  0); 


if(ik  ==  0  |]  ik  >  97) 


r 


ERROR 
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printff'error:  k’2  =  %t  is  out  of  range  for  this  table  0,k); 
exit(l); 

} 

dl  =  f[ik  +  l]  -  f [ik]; 

d2  =  f[ik+2]  -  2.0*f[ik+l]  +  f [ik]; 

value  =  f[ik]  +  kdelta*(dl  -  (1.0  -  kdelta)*d2/2.0); 

return(value); 


APPENDIX  C 


lr 


DATA  ACQUISITION  CODE  MAINTAINACE  GUIDE 


OVERVIEW: 

The  data  acquisition  package  was  designed  to  be  very  flexible  in  allow¬ 
ing  the  Tecmar  Labmaster  board  and  IBM  PC/AT  to  collect  data  and  con¬ 
trol  experiments.  The  basic  system  is  shown  in  Figure  Cl.  The  package  is 
capable  of  presenting  and  recording  8  bits  of  stimulus  and  response  data 
while  timing  two  sets  of  external  events  and  digitizing  eight  channels  of 
data.  Warning  tones  can  also  be  presented  to  alert  the  user  to  prepare  for  a 
stimulus.  All  of  the  data  collection  is  interrupt  driven  which  allows  the  user 
to  perform  real  time  analysis  of  the  data  while  it  is  being  collected.  The 
A/D  channels  are  double  buffered  so  that  the  previous  trial  can  be  analyzed 
while  the  data  for  the  current  trial  is  being  collected.  The  program  is  con¬ 
trolled  by  the  input  parameter  file  which  may  be  prepared  in  advance  and 
modified  at  run  time.  The  makepar  program  is  a  program  that  will  generate 
a  parameter  file  based  on  prompting  the  user  for  inputs.  The  data  collection 
of  the  D/A  channels  is  broken  down  into  trials  and  the  user  can  have  as 
many  trials  as  disk  space  will  allow.  The  external  interrupt  lines  are  moni¬ 
tored  continuoulsy  during  the  course  of  the  experiment  and  are  independent 
of  the  trials.  Each  trial  can  present  two  independent  8  bit  stimuli  and 
record  two  8  bit  responses.  The  responses  interrupt  handler  can  be  used  to 
turn  the  stimuli  off.  Up  to  64K  bytes  of  digitized  data  can  be  recorded  per 
trial. 

The  timing  for  an  individual  trial  is  shown  in  Figure  C2.  The  delay 
time  (in  msec),  foreperiod  (in  msec),  ISI  (in  sample  points)  and  warning  tone 
frequency  can  be  varied  on  a  trial  by  trial  basis.  Both  stimuli  can  also  be 
varied  on  a  trial  by  trial  basis.  The  warning  tone  duration,  sampling  length 
and  prestimulus  interval  (for  the  first  stimulus)  are  fixed  at  the  start  of  each 
experiment.  All  of  these  variables,  as  well  as  the  sampling  rate,  number  of 
responses  and  number  of  external  interrupt  lines  are  controlled  by  the 
parameter  file. 


TECMAR  CONTROL  USING  MNEMONIC  CODES: 


Much  of  the  mystery  of  programming  the  teemar  lab-master  card  has 
been  alleviated  through  the  use  of  mnemonic  definitions  contained  in  the 
teemar. h  and  param.asm  files.  The  teemar. h  file  is  not  an  exhaustive  list  of 
command  codes,  however  it’s  scope  is  more  than  adequate  to  handle  all  of 
the  command  codes  neccessary  for  the  current  software.  The  param.asm  file 
contains  a  subset  of  the  definitions  that  are  in  the  teemar. h  file  along  with 
general  definitions  for  the  assembly  language  routines.  General  definitions 
for  the  C  language  routines  are  in  the  param.h  file. 

The  first  section  of  teemar. h  contains  the  address  of  all  the  internal  tec- 
mar  registers  along  with  their  mnemonic  equivalent.  The  remainder  of  the 
file  consist  of  general  command  codes  as  well  as  application  specific  codes  for 
these  registers.  An  example  of  a  general  mnemonic  definition  from  teemar. h 
is: 

#define  LOADTX  0x08  /*  T_CONT  =  LOADTX  ;  TX;  T_DATA  ->  x  load  reg  */ 

The  convention  of  using  capital  letters  for  C  constants  is  followed  and  also 
used  for  assembly  language  constants.  The  Ox  prefix  is  the  standard  prefix 
to  indicate  a  hexadecimal  number.  Functionally,  this  statement  only  sets 
the  constant  LOADTX  equal  to  8.  Organizationally,  this  allows  us  to  write 
mnemonic  rather  than  cryptic  code.  The  comment  to  the  right  of  the 
definition  explains  its  use.  When  the  timer  control  register  (T_CONT)  is  set 
to  the  command  code  for  load  register  selection  (LOADTX)  is  or’ed  with  the 
register  number  (Tl,  T2,  T3,  T4  or  T5),  the  timer  data  register  will  address 
the  specified  load  register  (1,  2,  3,  4,  5  respectively).  An  example  of  the  use 
of  this  definition  is: 

outb(T  CONT, LOADTX  j  T5);  /*  point  timer  data  rcg.  to  timer  5  load  reg  */ 

out2b(T  DATA,  INIT  DELAY);  /*  set  timer  5  load  reg  to  the  initial  delay  */ 

T5  is  just  a  constant  equal  to  5.  This  may  seem  like  an  excesivc  use  of 
mnemonics,  however  some  command  codes  identify  timers  by  bit  number. 
The  command  code  to  arm  and  load  timer  5  (A_LSX  j  S5)  is  one  such  exam¬ 
ple.  Sf>  is  equal  to  0x10  (bit  5  set),  while  T5  is  equal  to  0x05  (the  number  5) 
yet  both  are  used  to  identify  timer  5,  but  to  different  command  codes. 


Mnemonic  codes  are  used  for  all  references  to  the  tecmar  board,  with 
the  execption  of  the  number  0,  to  enhance  the  readability  of  the  code. 
While  it  may  be  possible  to  program  the  tecmar  board  soley  through  thre 
use  of  these  mnemonics,  nothing  can  replace  an  understanding  of  how  the 
board  actually  functions.  The  command  "outb(P_CONT,  PA_2M);"  which 
configures  port  A  for  mode  2  operation  does  not  tell  you  very  much  if  you  do 
not  know  what  mode  2  is.  The  mnemonic  codes  can  also  be  a  valuable  asset 
while  reading  the  tecmar  manual  as  they  can  help  clarify  some  of  com¬ 
mands.  NOTE  that  approximately  20%  of  the  mnemonics  have  not  been 
used  in  the  current  software  and  are  therefore  untested.  One  may  wish  to 
verify  the  formulation  of  specific  mnemonics  before  using  them  for  other 
applications. 

MEMORY  ALLOCATION: 

The  1B.M-PC  operating  under  the  DOS  3.x  operating  system  is  a  16  bit 
machine  and  consequentially  a  data  segment  contains  64K  bytes  of  storage. 
The  small  model  C  compilers  assume  a  single  segment  for  data  storage  and 
therefore  use  16  bit  pointers.  This  is  not  adequate  for  our  needs,  so  that  the 
large  model  compilers  must  be  used  which  allow  for  multiple  data  segments 
through  the  use  of  32  bit  pointers.  It  is  impcritive  to  remember  that  all 
pointers  contain  both  a  segment  and  offset  value  when  interfacing  to  assem¬ 
bly  language  routines. 

There  are  four  types  of  data  structures  which  are  used  in  the  program 
and  defined  in  the  param.h.  The  param.asm  file  contains  an  equivalent 
definition  of  the  IC_STRUCT  data,  since  this  is  the  only  structure  used  by 
the  assembly  language  code. 

1)  INT_STRUCT  interrupt  structure  which  is  used  to  store  the  original 

ibm  interrupt  vectors  before  revectoring  by  the  dac  program 

2)  RESP_STRUCT  response  structure  used  to  store  the  stimuli, 

responses  and  responses  time 

3)  PAR_STRUCT  parameter  structure  which  is  used  to  hold  all  of 

the  parameters  which  specify  the  program  operation 
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4)  IC_STRIJCT  interrupt  control  structure  is  used  to  pass  information 
between  the  high  level  C  routines  and  the  low  level 
assembly  interrupt  handlers.  It  is  crucial  that 
the  C  structure  definition  in  param.h  match  the 
assembly  definition  in  param.asm  EXACTLY. 

The  dac.c  file  allocates  storage  for  the  program  parameters,  responses 
and  interrupt  revectoring  using  the  first  three  data  structures.  The 
storage. asm  file  allocates  several  segments  for  storage,  including  the  dataseg- 
ment  which  contains  the  interrupt  control  structure.  Two  64K  segments 
(buffferO  and  bufferl)  are  allocated  to  store  the  data  from  the  A/D  conver¬ 
sions  in  2  byte  interger  form.  On  odd  trial  numbers,  data  is  written  to 
bufferl,  while  on  even  trials  it  is  written  to  bufferO.  This  double  buffering 
allows  us  to  write  out  the  results  from  trial  n-1  while  data  is  collected  for 
trial  n.  Two  smaller  segments  (externalO  and  externall)  are  allocated  to 
provide  storage  for  the  timing  of  external  events  in  4  byte  long  interger 
form.  Storage  is  allocated  in  the  storage. asm  file. 


INTERRUPT  HANDLEING: 


There  are  three  types  of  interrupts  (timer,  response,  and  external)  that 
may  occur  and  each  type  has  its  own  interrupt  handler.  The  interrupt 
request  lines  are  defined  in  the  param.h  file  and  the  comments  in  the 
tec_init  routine  show  how  the  board  should  be  jumpered  for  the  current 
configuration.  The  tec_init  routine  initilizes  all  of  the  tecmar  control  regis¬ 
ters  and  uses  the  set_int  routine  to  alter  the  normal  IBM  interrupt  mask 
and  vectors.  Any  premature  termination  of  the  program  which  takes  place 
after  the  tec_init  routine,  MUST  call  the  error  routine  to  reset  the  IBM-PC 
to  it’s  original  state.  Failure  to  restore  the  interrupt  masking  &  vectoring 
will  result  in  ERRATIC  MACHINE  BEHAVIOR.  This  can  also  occur  if  one 
"breaks"  out  of  the  program.  The  IBM-PC  can  be  rebooted  to  restore  the 
IBM-PC  to  it’s  correct  state.  A  real  time  32  bit  clock  (timers  3&4)  running 
at  10  KHz  is  started  at  the  beginning  of  the  experiment  and  is  used  by  all  of 
the  interrupt  routines  to  record  the  time  of  events  (e.g.  stimulation, 
response,  external  interrupts).  The  clock  has  an  accuracy  of  0.1  milliseconds 
over  a  period  that  can  be  as  long  as  119  hours. 


TIMER  INTERRUPT: 


The  tm_int  routine  is  the  interrupt  handler  for  timer  5.  Timer  5  is  used 
for  both  delay  and  sample  timing  and  is  operated  in  a  count  down  mode 
which  generates  an  interrupt  when  the  count  reaches  zero.  If  the  sampling 
bit  (F_AD)  in  the  flag  byte  of  the  interrupt  control  structure  (ic. flags)  is  set, 
then  a  set  of  samples  from  the  A/D  converters  is  taken  (200  microseconds  to 
sample  6  channels).  The  routine  then  examines  the  sample  number  to  deter¬ 
mine  if  a  stimulus  should  be  turned  on  or  if  this  was  the  last  set  of  samples 
to  take  for  this  trial.  Note  that  the  second  stimulus  sent  to  the  output  port 
(Port  B)  is  equal  to  the  current  stimulus  in  the  port  register  exclusively  or’ed 
with  the  second  stimulus  stored  in  the  parameter  structure.  This  allows  the 
stimuli  to  be  set  and  reset  (by  the  response  routine)  independently  of  each 
other;  or  to  use  the  second  stimulus  to  reset  the  first  stimulus  by  using  the 
same  bits  for  both.  The  time  of  the  first  stimulus  is  read  from  the  real  time 
clock  and  stored  to  faciliate  determination  of  response  times  and  data  align¬ 
ment  with  externally  timed  events.  The  time  of  the  second  stimulus  can  be 
deduced  from  the  number  of  sample  points  between  the  stimuli.  If  the  sam¬ 
pling  bit  is  not  set,  then  the  routine  assumes  that  only  a  delay  was 
requested  and  sets  the  delay  finished  flag  (FJDELAY).  Since  no  sampling  is 
involved,  this  type  of  interrupt  can  be  handled  in  approximately  30  micro 
seconds.  The  interrupt  line  for  a  delay  is  generated  from  a  1  KHz  clock,  so 
that  this  line  is  forced  low  within  the  interrupt  routine  to  prevent  the  same 
clock  cycle  from  causing  multiple  interrupts. 


RESPONSE  INTERRUPT: 

The  res_int  routine  is  used  to  record  both  the  response  time  and 
response  byte.  It  will  record  up  to  two  responses  per  trial  and  the  interrupt 
from  the  second  response  is  ignored  if  the  second  response  is  the  same  as  the 
first  response.  This  means  that  if  you  press  the  first  response  five  times  and 
then  press  the  second  (different)  response,  only  the  first  time  that  the  first 
response  was  pressed  and  the  time  of  the  second  response  will  be  accepted 
and  stored.  The  current  response  byte  is  compared  to  the  appropriate 
stimulus  byte,  if  they  match  the  current  stimulus  is  and’ed  with  the  inverse 
of  the  current  response.  Note  that  this  is  different  from  simply  clearing  the 
stimulus  byte  and  allows  independent  operation  of  the  two  stimuli  and 
responses.  Flags  (ic. flags  byte,  bits  F_RESP0  &.  F_RFSPl)  are  used  to  indi¬ 
cate  which  response  byte  is  active  and  interrupts  rccievcd  with  no  active 


response  flag  byte  are  ignored.  Both  response  flags  are  set  at  the  start  of 
sampling  so  that  responses  given  during  the  prestimulus  interval  are 
recorded  with  a  negative  response  time.  This  allows  us  to  check  for  antici¬ 
pation  on  the  part  of  the  subject.  When  we  write  the  modified  stimulus  to 
Port  B  (output  port),  we  also  clear  the  interrupt  line  associated  with  B. 
This  can  pose  a  serious  problem,  since  the  Port  B  interrupt  line  is  used  to 
indicate  an  external  event  on  one  of  the  two  data  bits  assigned  to  Port  C. 
This  situation  is  discussed  further  in  the  external  interrupt  handler  section 
and  the  hardware  section  which  describes  the  hardware  solution  to  this 
problem.  The  hardware  is  configured  so  that  any  button  push  will  latch  the 
switch  configuration  into  Port  A  (input  port)  and  generate  an  interrupt. 
Switches  were  also  made  bounceless  by  using  an  RC  charge  and  discharge 
network  in  combination  with  a  single  pole  single  throw  switch.  The  response 
is  read  from  the  Port  A  register  and  the  time  is  read  from  the  real  time 
clock. 

EXTERNAL  INTERRUPTS: 

The  ext_int  routine  is  the  interrupt  handler  that  is  used  to  respond  to 
external  interrupts  occuring  on  the  two  data  bits  assigned  to  Port  C  (bits  C6 
&  C9)  and  generated  using  the  Port  B  interrupt  request  line  (INTRb).  Port 
B  is  currently  configured  as  an  output  port  and  drives  a  series  of  LEDs. 
Port  C  is  currently  used  to  threshold  detect  the  EKG  and  respiratory 
waveform.  The  interrupt  request  line  on  Port  B  in  this  configuration  was 
designed  to  function  as  a  ready  to  recievc  signal.  Since  we  are  only  driving 
LEDs  we  do  not  need  this  function  and  instead  use  the  Port  B  interrupt  to 
indicate  that  an  external  event  has  occurred  on  one  of  the  two  data  bits  of 
Port  C.  Unfortunately,  every  time  that  we  write  to  Port  B  we  clear  the 
interrupt  request  line.  A  hardware  solution  (see  hardware  section)  was  con¬ 
structed  based  on  generating  a  second  interrupt  on  Port  B  if  the  first  inter¬ 
rupt  line  was  cleared  in  less  time  than  it  took  for  the  ext_int  routine 
(approximately  30  microseconds)  to  handle  the  interrupt.  This  solution  is 
satisfactory  proved  that  the  Port  B  interruprt  request  line  priority  is  at  a 
higher  level  than  the  other  interrupt  priorities.  Port  C  only  has  two  data 
bits  available  since  the  other  six  bits  are  used  for  handshaking  on  Ports  A 
and  B.  The  interrupt  handler  simply  reads  in  Port  (■  and  determines  which 
external  event  occured  (bit  C6  or  C7  set)  to  determine  which  data  segment 
to  store  the  time  of  the  event  in. 


HARDWARE: 


An  interface  box,  shown  in  Figure  Cl,  was  contructed  to  connect  the 
Tecmar  data  acquisition  to  the  outside  world.  The  input  port  (Port  A)  and 
the  output  port  (Port  B)  were  both  buffered  in  this  box  to  prevent  damage 
to  the  Tecmar  board  if  an  improper  connection  was  made.  The  main  circui¬ 
try  in  the  box  was  the  interrupt  arbitrator  circuitry  shown  in  Figure  C3. 
This  circuit  served  two  purposes.  The  circuit  had  to  be  able  to  correctly 
handle  simultaneous  interrupts  on  both  external  interrupt  input  lines.  It 
also  had  to  address  the  problem  of  using  the  Port  B  interrupt  line  for  data 
on  Port  C.  This  problem  was  addressed  in  was  discussed  in  the  external 
interrupt  section.  The  term  simultaneous  can  be  defined  as  two  events 
occuring  within  the  time  required  to  service  a  single  interrupt  (approxi¬ 
mately  30  microseconds).  This  is  a  more  than  adequate  definition  of  simul¬ 
taneity  since  the  real  time  clock  only  measured  time  in  increments  of  100 
microseconds. 

The  resistor  capacitor  combination  at  the  input  of  the  circuit  is  used  to 
set  the  threshold  level  and  provide  some  low  pass  filtering  of  the  input  lines. 
Note  that  in  this  configuration  the  respiration  channel  was  triggered  by  the 
start  of  an  expiration  and  the  EKG  channel  was  triggered  by  the  rising  edge 
of  the  R  wave  if  the  input  polarities  were  sert  properly.  The  Schmitt  trigger 
inverter  (U9)  is  used  to  prevent  multiple  threshold  crossings  from  signals 
with  a  slow  slew  rate.  The  multivibrator  (137)  and  flip  flop  (U4)  were  used  to 
handle  the  problem  of  pseudo  simultaneous  threshold  crossings.  The  line 
which  crosses  the  threshold  second  will  generate  a  second  interrupt  as  soon 
as  the  system  finishes  with  the  first  interrupt.  The  operation  of  this  circuit 
was  tested  and  verified  with  pulses  that  were  seperated  by  1  microsecond. 
The  next  set  of  multivibrators  (U8)  are  used  to  reinterrupt  the  system  if  the 
interrupt  line  was  cleared  prematurely  by  one  of  the  other  routines  writing 
to  Port  B.  This  circuit  was  successfully  used  on  10  different  subjects  and 
the  only  problems  encountered  were  with  partial  breaths  which  did  not  cross 
the  threshold  level. 
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SUBROUTINE  -  CHECKPAR(fp) 


PURPOSE: 

check  the  parameter  structure  &.  allow  user  to  review 


INPUT: 

external  -  par:  parameter  structure  defined  in  param.h 
OUTPUT: 

external  -  par  (modified) 


DEPENDENCIES: 

get)ine()  -  function  used  to  read  in  1  line  at  a  time 

sldio.h 

inath.h 

param.h 


/ 


rf  include  ''stdio.h  - 
^include  <math.h  • 
if  include  "param.h" 

#  define  gctp(form.v)  if(gctfine(s)  O)sscai.f(s,form,v) 
void  rheckpar() 


extern  PAR.  STRUCT  par; 
int  rO,  c  1 ; 
char  sjMAXLINE); 
int  i, 

min  delay,  max^delay, 
min.wtper,  max_wtper, 
min.  wtfore,  max_wtfore, 
min  isi,  max  isi; 


/*  parameter  structure  */ 

/*  temp  char  to  make  stimulus  word  */ 

/*  string  for  inputting  parameters  */ 

/*  temporary  index  */ 

/*  min  &'  max  IS]  delays  */ 

/*  min  &  max  of  warning  tone  periods  */ 
/*  min  &  max  of  warn  tone  foreperiod  */ 

/*  min  max  of  H  samps  between  slims  */ 


printf("Do  you  wish  to  review  the  parameters  selected  7  ); 
get  line  (s); 

whilefs'O]  ’Y’  1!  s|0)  —  ’y’){ 

printff'Entcr  new  value  or  CR  to  keep  old  0); 

printff "number  of  channels  Ct  d  7  .par.nclian); 
getp("%d",4rpar.nchan); 

printff'number  of  responses  expected  (0,1  or  2)  ^od  7  , par. response); 


getp("<7id",A  par.  response); 


printff'number  of  external  interrupt  lines  (0,|or  2)  %<i  ?'‘,par.cxtcrnal); 

getp(  par. ex  ter  nal); 

printff'number  of  trials  (max  ?<’d)  ‘"<1  T,MAX  TRIAL, par. ntrial); 

getpf' ffd'.Apar.n trial); 

printff'sampling  period  in  usee  (max  -  G-1000)  ‘7d  ?",par.s.  per); 

getpf'^rd  ",&par.s  per); 

printff  warning  tone  duration  in  msec  '7d  par.wt  dur); 
get  pf’fcd  ' par.w  t  dur); 

print f(  "number  of  samples  per  trial  rfd  7  .par  s  total); 
getpf'ftid  '.At  par.s  total); 

printff'sample  #  after  which  to  turn  stimulus  0  on  (7d  7", par.s  stimO); 
getp(  Tud’  .&par.s  stimO); 

min  .delay  -  max  delay  par. delays  ();;  /  *  initialize  min  A  max  */ 

min.  wtper  —  max  wtper  par.wt  persO  ; 

min.wtfore  max  wtfore  par.wt  fore  0  ; 

min  isi  max  isi  par.s  isi!0  ; 

for  (r  0;  i  par. n t rial;  i  •  t  ) [ 

min  delay  minfmin  delay,  par. del, ays  i;); 
max  delay  maxfmax  delay,  par.delay.7i  ); 
min  wtper  minjtnin  wtper,  par.wt  persji  ). 
max  wtper  max(max  wiper,  par.wt  pers  i  ); 

min  wtfore  min(min  wtfore,  par.wt  fore  i  ); 

max  wtfore  inaxjmax  wtfore,  par.wt  f..re  i  ); 
min  isi  minjmin  isi,  par  s  isi  i  ); 
max  isi  max(max  isi,  par.s  isi(i ; ); 

l 

printff  1st  trial  stimulus  '7x  A  last  trial  stimulus  '7x  %x0, 

par.stimslOj  A  OxOOfT, (par.stims[0j  •  •  H)  A  OxOOIf.par.stims  par.ntrial 
-1;  A  OxOOfT,  (par.stimsjpar.ntrial  -1,  •  H)  A  OxOOIf); 

printff  "ISI  (in  ft  samples)  range  from  rfjd  to  (7d  O.tnin  isi, max  isi); 
printff  inter  trial  delay  range  from  'Yd  to  ‘  fid  insert), min  delay, max  delay); 
printf( "warning  tone  periods  go  from  'Yid  to  r<’d  insert), min  wtper, max_wtper); 
printff  VVT  fore  periods  range  from  ‘7d  to  '7d  nisei  0, min  w  tTore, max.  wtfore); 
printfj  'Do  you  wish  to  review  all  f7d  trials  par.ntrial); 
getline(s); 

if(s;0)  -  y  ;  so;  ’Y  ){ 

printf(  ”0im(hex  pair),  I S I ( i/  samps),  trial  delay(ms),  \VT  period(us),  WT  foreperiod(ms))0); 
forfi  0;  i  par.ntrial;  i  *  *  )j 

print f( ’(‘'ex  rrX,  ‘Yd,  *7d,  ‘  <.d,  *7 d)  7 ".par.stims  i  A 

OxOOfT, (par. stims1ij  •  ■  R)  A  OxOOfT, par.s  isijij, par. delays  ij, 
par.wt  prrs[ii, par.wt  forejii); 
if(getline(s)  1  0)( 

sscanffs,  "'7x  ro»  fed  *7  d  ‘7d  '7  d",  .DO,  Arl,  Apar.s  isii’, 

Apar. delays  i  ,  Apar.wt  prrsjil); 


e.  e.  s 


par.stimslij  —  cO  [  ((  cl  ■'  •  8)  &  ~OxFF); 


printf)"  Oeview  Parameters  Again  7"); 
getline(s); 


/*  return  to  top  of  loop  */ 


CHECK  RANGE  OF  PARAMETERS  */ 


if( par . response  >  2)  error)  "tnkpar:  par. response  2"); 

if ( p a r . external  2)  error)  rlieckpar:  par. external  >  2” ); 

iffpar.ntrial  >  MAX.  TRIAL  (error)  checkpar:  p.ir.r.trial  •  MAX  TRIAL  ) 

for ( i  —  0;  is  par.ntrial;  i  +  -t ){ 

if(par.delavs!ij  par.wt.dur  4  par.wt.  fore [i i ) 

error)  checkpar:  par. delays  •  par.wt.dur  4  par.wt  foreperiod"); 
if) pa r .s  i s i ! i .  ;  0) 

error)'  checkpar:  par.s_isi  =-  0"); 

iflpar.s  is i  i  ;  par.s  stimO  par.s  total) 

error)  checkpar:  par.s  stOl  -t  par.s. stimO  par.s  total  ); 


PROGRAM  - DAC 


i 

i 


PURPOSE: 

data  acquisition  program 


. * . . . * . / 

<f  include  <  stdio.h> 

(/include  "termar.h" 

(/include  "param.h" 

PAP  STRUCT  par;  /*  allocate  space  for  parameter  structure  */ 

RESP  STRUCT  resps[MAX„TRLAL|;  /*  allocate  space  for  responses  £  times  */ 

INT  STRUCT  tm,  /*  allocate  space  for  timer  int  vector  */ 

res,  /*  allocate  space  for  response  int  vectors  */ 

ext;  /*  allocate  space  for  external  int  vectors  */ 

main)) 


setup!); 
tec  inil(); 
while  (samp())i 

/*  enter  task  to  be  done  while  sampling  here  */ 

delay!); 

/*  enter  task  to  be  done  during  delay  here  */ 


tec  donef); 

storef); 

exit(O); 


/*  reset  £  disable  tecmar  £  int  */ 
/*  store  £  close  all  Dies  */ 

/*  normal  program  exit  */ 


FUNCTION  SUBROUTINE  -  DELAY() 


. . . . . . . .  1 

1 

l 

PURPOSE:  ! 

Perform  all  neccessary  operations  after  sampling  is  • 

completed  for  a  trial.  j 

1)  sampling  done  flag  set  ->  call  error  &  terminate 
sampling  finished  before  high  level  proccesing  —  fatal  error 

2)  wait  for  sampling  done  flag  to  be  set 

3)  reset  sampling  done  flag 

5)  check  if  this  was  last  trial  , 

7)  start  delay  in  background  j 

INPUT:  ! 

none  -  external  ic,  par  ! 

OUTPUT:  I 

f 

none  ( 

i 

RETURN  (int):  j 

I 

1  =  normal  sampling  -  inter  trial  delay  A-  warning  tones  initiated  j 

0  —  normal  finish  -  no  delay  initiated  i 

-  use  startsamp  to  store  last  trial  j 

DEPENDENCIES:  | 

tin  intj)  -  performs  actual  delay  in  background 
errorj)  -  error  &  terminate  routine 
ic  -  interrupt  delay  structure  (extern) 

par  -  experiment  parameter  structure  (extern) 
teemar.h  -  parameters  for  teemar  board 
param  h  -  general  parameters 
stdio.h  -  C  standard  input/output 


#  include  ■  stdio.h  ■ 
# include  "param. h" 
H  include  "teemar.h" 

int  delay() 


extern  PAR  STRUCT  par, 
extern  IC  STRUCT  ic; 
int  temp; 


/*  experiment  parameter  structure  */ 
/*  interrupt  control  structure 


/*  CHECK  SAMI’!. INC.  DONE  FL  .O  */ 


if(ic. flags  ft-  F.  SAMI')  error)  "cndsamp:  sampling  finished  early  "); 


u 

A 


r  WAIT  FOR  SAMPLING  FLAG  ft  THFN  RESPT  IT  */ 


whilc(!(ic .flags  ft;  F_SAMP)); 

ic .flags  &=  ~(F  SAMP  !  F  AD  |  F  RESPO  ;  F_RESP1); 

/*  clear  samp  done,  a/d,  waiting  f< < r  response  flags  */ 


/*  CHECK  IF  FINISHED  WITH  TRIALS  -  RETURN  (0)  -  last  trial  not  stored  yet  */ 


if  (+  t  par.itria)  =  —  par.ntrial)return(O); 


/*  START  DELAY  ft  WARNING  TONES  */ 


print  f(  "starting  delay  for  trial  ji  /<  d  0. par. it  ria  1); 

out  l.(T_  CONT,T5); 

<>ul2l)(T_ DATA,TM5ms);  /*  set  timer  a  for  msec  delays  */ 

out2h(T.  DATA, par. delays'par.itrialj);  /*  set  de  lay  in  load  reg*/ 
out(T_  CONT.A  _LSX  Sn);  /*  start  delay  */ 


if  (par.wt.dur  !—  0)]  /*  produce  warning  tones  ‘/ 

outl>(T  .CONT.LOADTX  !  Tl);  /*  T  DATA  •  timer  1  load  reg  */ 
temp  -  par.delay.sjpar.itrial]  -  pa  r.wt..  f«re ‘par. it  rial!  -  par.wt.dur  + 

(par.s  per  *  par.s  stirnO)  /  1000  ; 
out2b(T  DATA,  temp); 

outh(F  CONT.LOADTX  '  T2);  /"  T.DATA  -  timer  2  load  reg  */ 
out2b(T  DATA.  par. wt.  porsjpar.it  rial  ); 

out  b(T_CONT,  A_LSX  ;  SI  S2);  /*  start  warning  tone  generator  */ 

iff!  I IGH.FREQ  •  -  Fl  /  (long)par.wt  pers  par. it  rial  ){ 

outb(HU  ODA.OxOh);  /*  output  0  volts  (0  -  10  range)  */ 


outbfl.H  ODA.OxOO); 

■  else | 

outb(HB  ODA.OxOO); 
outbjlJL  0DA.  0x00); 


/ ‘  output  avohs  (0  -10  range)  */ 


/*  RETURN  NORMALLY  */ 


return)  1  ); 


mmmss&samiftBp 


N  % 


SUBROUTINE  -  ERROR(msg) 

4444444444444444444444444444444444444444444444444444444444444444444 

PURPOSE: 

print  error  message  Sr.  terminate  program 

INPUT: 

msg  -  pointer  to  error  message  to  be  printed  out 
OUTPUT: 

none  -  execution  terminates 

DEPENDENCIES: 

exi t ( )  —  standard  C  exit  routine 

tecmar.h  —  tecmar  constant  definitions 

param.h  —  general  parameter  &  structure  definitions 

stdio.h  —  C  standard  input/output 

******************************************************************* 


it  include  <stdio.h  '  ■ 

^/include  "param.h" 

%  include  "tecmar.h" 

void  error(msg) 

char  *msg;  /*  error  message  to  be  printed  out  */ 

{ 

extern  PAR_STRUCT  par; 

outb(T_CONT,  DISAUTO);  /*  disable  timer  interrupts  * / 

outb(P_CONT,  PC_RSET  ;  1NTE_A1N);  /*  diable  PORTA  interrupts  */ 

outb(P_CONT,  PCLRSET  ;  INTE..B);  /*  disable  PORTB  interrupts  */ 

if(inb(INTA01 )  k  BITO)  outb(INTAO]  ,par.int^mask);  /*  restore  8259  mask  * / 
fprintf(stderr,  "%s  0,  msg); 

exit(-l);  /*  abnormal  program  termination  */ 


t  W VI  "wV'.'A.T'S.NNiAA Wl-MS V  .Wl 


*  JL*Wi 

title  EXTERNAL  INTERRUPT  HANDLER  -  extjnt 


» 


I 

J 

include  pa.ram.asm  j 

codeseg  segment  para  public  ’code’  I 

assume  cs:codeseg  j 

public  ext_int_,extO_p_,extl_p_,ext_f_ 


extO_p_ 

label  byte 

;acces  offset  &  segment  using  C  pointer 

extO_off  dw 

0 

joffset  of  extO  buffer 

extO_seg 

dw  0 

;segment  of  exlO  buffer 

extl_p_ 

label  byte 

;acces  offset  &  segment  using  C  pointer 

extl_offdw 

0 

joffset  of  ext]  buffer 

extl_seg 

dw  0 

jsegment  of  ext]  buffer 

cxt_f_ 

label  byte 

;cxt  int  Dag  =  #  external  int  lines 

ext_f 

db  0 

ext_int_  proc 

far 

;entry  point  for  interrupt  (ext.int  in  C) 

pushreg 

;save  ax,  di,  ds,  dx,  es 

cmp 

cs:ext_f,F .JCXTO  ;Don’t  bother  to  read  PORTC  if  only 

je 

get_exO 

;1  external  interrupt  line 

mov 

dx,  PORTC 

in 

al,dx 

;get  port  C  input  lines 

test 

al,F_EXT0 

;check  if  PORTC  bit  F_EX0  is  high 

j* 

ck_exl 

get_exO: 

getime 

cs:extO_off 

;store  time  of  interrupt  in  extO  buffer 

mov 

cs:extO_o£T,di 

;inc  offset  to  point  to  next  free  byte 

cmp 

di,  BUF_EXT0 

;check  if  next  write  will  overflow  buf 

jb 

ext_int_exit 

mask_off: 

in 

al,  1NTA01 

;  get  current  interrupt  mask 

or 

al,  P0RT13_IRQ_BIT  ;  mask  off  further  interrupts 

out 

INTA01,  al 

jmp 

short  ext_int. 

_exit 

ck._cxl : 

test 

al.FJiXTl 

;check  if  PORTC  bit  F_EXl  is  high 

ext_int_exit 

i 

i 

i 

i 

i 


getime 

cs:extl_off 

;store  time  of  interrupt  in  extl  buffer 

mov 

cs:extl_off,di 

;inc  offset  to  point  to  next  free  byte 

cmp 

di,  BUF_EXTl 

;check  if  next  write  will  overflow  buf 

jnb 

mask_off 

ext_int_exit: 

mov 

dx,PORTB 

in 

al,dx 

;get  current  stimulus 

out 

dx,al 

;reset  PORTB  interrupt 

mov 

al,  EOI 

out 

INTAOO,  al 

jacknowledge  interrupt  to  8259 

popreg 

jrestore  ax,  di,ds,dx,es 

iret 

ext_int_  endp 

codeseg  ends 
end 


******** 


FUNCTION  SUBROUTINE  GETLINE(s) 


PURPOSE: 

getline  will  get  a  line  of  input  and  return  the  number  of 
characters  in  the  input  line.  0  =  blank  line 

INPUT: 

s  =  character  array  of  length  MAXLINE 
OUTPUT: 

s  =  iput  line  from  standard  input 

RETURNS: 

length  of  input  line 

0  =  carriage  return  only  -  no  blanks  or  characters 
DEPENDENCIES: 


^include  <stdio.h> 

^include  "param.h" 

int  gctlinc(s) 
char  s[); 

{ 

int  c, 

i  =  0, 

lim  =  MAXLINE; 

while(—  lim  >  0  &&  (c  =•■  gctchar())  !=  EOF  &.&: 


'y  V'V.*;'-  ...  V  V'  ‘  /  '  ■  V 


I 


V* 

$ 

$ 

.•2 


,y<i 


a 


I 

L 

si 


Si 


J* 

i 

$ 


ji 


| 

j! 


1 


SUBROUTINE  -  GETPAR(fp) 


PURPOSE: 

get  the  parameter  struct  from  the  default  file 


INPUT: 


fp  -  filename  containing  parameters  for  input 
external  -  par:  parameter  structure  defined  in  param.h 


OUTPUT: 

external  -  par  (modified) 


DEPENDENCIES: 

gctline()  -  function  used  to  read  in  1  line  at  a  time 

stdio.h 

math.h 

param.h 


^include  <stdio.h> 
^include  <math.h> 
^include  "param.h" 


void  gctpar(fp) 


FILE  ‘fp; 


/*  parameter  file  pointer  */ 


extern  PAR_STRUCT  par;  /*  parameter  structure  */ 

int  cO,cl;  /*  temp  byte  to  make  stimulus  word  */ 

int  i;  /*  temporary  index  */ 


if(fscanf(fp,"%d",&par.nchan)  =  —  0)  /*  number  of  channels  */ 

errorf'getpar:  unable  to  read  par.nchan"); 


if(fscanf(fp,"%d",&par. response)  =====  0)  /*  #  responses  expected*/ 

errorf'getpar:  unable  to  read  par.responsc"); 


if(fscanf(fp,"%d",&par.extcrnal)  =====  0)  /*  #  ext.  interrupts  PORTB_lRQ  */ 

error("gctpar:  unable  to  read  par. external'  ); 


if(fscanf(fp,"%d",&par.ntrial)  =====  0)  /*  number  of  trials  */ 

errorf’getpar:  unable  to  read  par.ntrial"); 


if(fscanf(fp,"%d",&par.s_per)  =====  0) /*  sample  period  in  usee  */ 
errorf'getpar:  unable  to  read  par.s.per"); 


if(fscanf(fp,"%d",&par.wt_dur)  =  =  0)  /*  warning  tone  duration  msec  */ 

error("getpar:  unable  to  read  par.wt_dur '); 


if(fscanf(fp,"?cd",&par.s_total)  =  =  0)  /*  number  of  sample  pts  */ 

errorf'gctpar:  unable  to  read  par.s_total"); 


if(fscanf(fp,"9c/d",&par.s_stim0)  ==  0)  /*  stim  0  on  after  stim  #  */ 

crror( "getpar:  unable  to  read  par.s. stiinO’  ); 


/♦  GET  TRIAL  DEPENDENT  DATA  &  CHECK  RANGE  ON  DEJ  AYS  */ 


for  (i=0;  i  <  par.ntrial;  in — h){ 

if  (fscanf(fp,"%x  %x  %d  %d  %d  %d", 

&c0,  &cl,  /*  stimulus  is  hex  pair  (two  each)  */ 

/*  1st  two  =  stimO,  2nd  two  —  stiml  */ 
&par.s_isi[i],  /*  #  samps  after  stimO  until  stiml  * / 

&par.deiays|i],  /*  delay  is  in  msec  * / 

&par.wt_pers[i],/*  wt_pers  —  period  in  usee  */ 

&par.wt_fore jij)  <  6)  /*  wt_forc  =  foreperiod  in  msec  */ 

crror("getpar:  unable  to  read  trial  dependent  data"); 
par.stims(i)  =  cO  ]  ((cl  <<  8)  &  ~0xff); 


fclose(fp); 


/*  CLOSE  KILE  V 


r.VtVi,*’  1  '  V  V  '  ’•  v,‘  »'  V  v  v  v  i  't  Li  i 


PROGRAM  -  MAKEPAR 


*4**44*t***t**l*4M**t****4»***4***4t4«»**4***44»4****l*4t*»*******«* 

PURPOSE: 

make  a  parameter  file  to  be  used  by  the  dac  program 
user  is  prompted  for  output  file  name  &  various  parameters 

DEPENDENCIES: 

gctline()  -  used  to  read  in  one  line  at  a  time 
writepar()  -  used  to  write  out  parameter  file 
error()  -  standard  error  writing  routine 
stdio.h 
parain.h 

44*44444444444444444444444444444444444444444444444444444*444444444444 


^include  <stdio.h> 
it  include  "parani.h" 

H  define  getp(form,v)  if(getline(s)  !=  O)sscanf(s,form,v) 

PAR. STRUC  T  par  =  {5,2,2,0,0,0,0,101,  /*  parameter  structure  */ 

0,4167,500,600,120}; 


char 

pc  ’n’; 

/*  pc  —  y  practice  case,  else  pc  —  n  */ 

char 

sjMAXLINEj; 

/*  string  for  inputting  parameters  */ 

int  i, 

seed; 

/*  temporary  index  k  seed  for  random  * / 

int 

rnd,temp; 

/*  temp  variables  for  randomising  * / 

int 

nchoice  =  4; 

/*  number  of  choices  in  this  block  */ 

int 

isi  =  15; 

/*  IS!  in  #  of  sample  pts  */ 

int 

wtper2  =  100; 

/*  2  choice  case  has  low  warn  tone  */ 

int 

wtper4  =  1000; 

/*  4  choice  case  has  high  warn  tone  * / 

int 

delay  [MAXJTR1AL]; 

/*  temp  for  randomising  delay  / 

int 

wtfore[MAX_TRLAL]; 

/*  temp  for  randomising  fore  periods  */ 

int 

stim_wt[MAX_TRLAL][2j; 

/*  temp  for  randomising  stimuli  &  wt  */ 

int 

smask2[2]  =  0x204,0x402; 

/*  stimulus  mask  for  2  choice  case  */ 

int 

smask4[l2j  =  0x102,0x104,0x108, 

/*  stimulus  mask  for  4  choice  case  */ 

0x201,0x204,0x208, 

0x401,0x402,0  .408, 

0x801,0x802,0x804; 

int 

min_delay  —  3500; 

int 

max_dclay  =  3500; 

/*  min  &  max  ISI  delays  * / 

int 

wmask(3]  =  1500,2000,2500; 

/*  choice  of  these  wt  foreperiod  */ 

FILE 

*fp; 

/*  file  pointer  for  output  file  */ 

printfj  Kilter  the-  output  parameter  file  name  ?"); 
getp("';t',s",s); 

ifjjfp  fopenjs,  w ’)  )  =  =  ■  NULL){ 

printfj  unable  to  open  output  file  name:  /os  0,s); 
exitj-l); 


s[0]  =-  V; 

while(!(sj0j  -  ’N’  ;;  s[0)  -----  ’n')){ 

printf("Enter  new  value  or  CR  to  keep  old  0); 

printfj  "number  of  channels  —  /7>d  ?",par.nchan); 
gctp("%d",&par.nchan); 

printfj"  number  of  external  interrupt  lines  (0,lor  2)  —  %d  ?",par.external); 
getpj  "f  od",&r  par. external); 

printfj  "number  of  trials  (max  -  %d)  -----  %d  ?",MAX_TRlAL,par.ntrial); 
getpj  "f,  cd",&par.ntrial); 

printfj  "sampling  period  in  usee  (max  —  C4000)  -  %d  ?",par.s_pcr); 
getpj" /'cd  ",&.pa r.s_ per); 

printff'warning  tone  duration  in  msec  =  %d  ?",par.wt_dur); 
getp("f/fd",£'par.wt_dur); 

printfj  "number  of  samples  per  trial  =  %d  ?",par.s_total); 
getp("f<  d  ",4'par.s.  total); 

printfj  "sample  $  after  which  to  turn  stimulus  0  on  =  %d  ?",par.s_stimO); 
get  p("/?d"I&par.s_stimO); 

printfj  "Trial  min, max  delay  (msec)  =  %d,  %d  ?",min„delay,max„dclay ); 
ifjgetlinejs)  !=  0) 

sscanf(s,"/od  %d",&rain_delay ,  &max_delay); 

printfj'  Foreperiod  values  (msec)  =  %d,  %d,  %d  ?",wmask[0],wmask[]],wmask[2]); 
if(getline(s)  !=  0) 

sscanfjs,"%d,%d,%d",&wmask[0],&wmask[l  ],&wmask[2]); 

printfj'Singlc  stimulation  practice  block  (y  or  n)  =  %c  ?",pc); 
gctp( ’’%c",&pc); 


if(pC  =--  y ;;  pc  =--  T’){ 

par. response  —  1;  /*  one  light  so  one  response  * / 

isi  -----  327G5;  /*  second  stim  will  not  come  on  */ 

|  else  ( 

printf(”ISI  in  //  of  sample  pts  (1  pt  ~  %d  usee)  =  %d  ?",par.s_per,isi); 
getp("%d",&isi); 


printff'Two  choice  (2),  mixed  (3)  or  four  choice  (4)  exp.  -  ‘  id  ?  '.nchoice); 
getp("9t,d",&-  nchoice); 

printff’  Oeview  Parameters  Again  ?'); 

getline(s);  /*  return  to  top  of  loop  */ 

1 

printf("Enter  a  seed  for  the  random  number  generator  ?"); 
scanf("%d  ",seed); 

srand(secd);  /*  initialize  random  number  generator  */ 

for (i  =  0;  i<  par.ntria);  i  ++){ 

delayji]  =  min_  delay  +  (long)  rand()  *  (max_delay  -  min_dclay)  /  32 7 C 7 ; 
wtfore[i]  —  wr.iaskji  %  3]; 

} 

for(i=i0;  i<par.ntrial;  i  +  +  ){  /*  randomize  wt_fore  periods  */ 

rnd  —  (long)  rand()  *  par.ntrial  /327 67 ; 
if(  rnd  =•  —  par.ntrial)  rnd--; 
temp  —  wtfore|i'; 
wtforcli]  —  wtforejrndj; 
wtfejrejrndj  —  temp; 

( 

switch  (nchoice)  { 

case  2:  /*  two  trial  case  */ 

for (i  —  0;  i< pa r.n trial;  i++)j 

stim_xvtjij(0j  =  smask2j  i  %  2]; 
stirn  wt(i][lj  =  wtper2; 

> 

break; 
case  3: 

for (i=  0;  i  <  par.ntrial;  i  +=  2){ 
stim  _wt [i] [0]  =  smask2[  (i/2)  %  2]; 
stim_wt[i][l]  =  wtper2; 
stiinjwtji  e  l)[0]  =  smask4[  (i/2)  %  12]; 
stiin_wt|i+l][l]  —  wtperd; 

1 

break; 
case  4: 

for(i~0;  i  <  par.ntrial;  i-M  ){ 

stim_wt(i][0j  —  smask4[  i  %  12]; 
stim_ wtjijjl]  --  wtpcr4; 

I 

break; 

default: 

printf("Error,  nchoice  \-  2,3  or  4  0); 


/*  randomize  stimuli  * / 


for(i"  0;  i-  par.ntrial;  i-r-t)| 

rnd  (long)  rand()  *  par.ntrial  /327G7; 

if (  rnd  —  =  par.ntrial)  rnd—; 

temp  _  stim.  wt |i][0l; 

stim  _ wt  ji j jo]  ~  stim  _wt  [rnd !j()J; 

stim.  wt|rnd]|0]  —  temp; 

temp  =-  stim.wt  jij|l|; 

stiin„wt(i![l]  =  stim  _wt[rnd;jl  i; 

stim  .  wt  [rnd]|l]  ==  temp; 


rnd  =  (pc  —  = 


V)?  OxOOfT  :  Oxffff; 


fry 

r-v 


m 

.v. 


for(i-  0;  i-  par.ntrial;  i  +  -t)(  /*  copy  to  parameter  structure  */ 

par.stims|i]  —  stim  wt|i',0j  A  rnd; 
par.s.  isiji]  =  isi; 
pnr.delays'ij  —  delay li ) ; 
par.wl  pcrsji]  —  st iin _ w t [ij ;  1  ’; 
par.wt.  forejij  —  wtfore[i;; 


/* 


CHECK  RANGE  OF  PARAMETERS  */ 


iffpar. response  >  2)  crror("makcpar:  par. response  >  2"); 
if(par. external  >  2)  error("inakepar:  par. external  >  2  "); 
if(par.ntrial  >  MAX_.TRlAI.)errorf'niakepar:  par.ntrial  >  MAX_TR1AL 
for(i— 0;  i<'par.ntrial;  i  +  +  ){ 

if(par.delays[i)  <  par.wl  dur  •  par.wt  fore (i j) 

error(  "makepar:  par. delays  >  par.wt  dur  t-  par.wt  foreperiod"); 
if(par.s... isiji]  = 0) 

errorf'makepar:  par.s.  isi  *==  0"); 
if(par.s__isi[i]  +  par.s_stim0  —  =  par.sjotal) 

error("makepar:  par.s  _isi  -t  par.s_stim0  =--  par.s_total"); 


writepar(fp); 

fclose(fp); 

ret  urn; 


INCLUDE  FILE  -  I’AKAMASM 


CONTAINS  AI.I,  NECESSARY  DEFINITIONS  FOR  TECMAR  BOARD 

IMPLEMENTED  AS  AN  I/O  DEVICE  AT  0710H  -  071FH 

R/W  AT  THE  END  01  THE  COMMENTS  REFER  TO  READ  &  WRITE 


m 

$] 

•  »t#: 

ts 


TECMAR  EQU 


07 1 0!  I; 


/*  TECMAR  STARTING  ADDRESS  */ 


/*  DIGITAL  /  ANALOG  CONVERTERS  */ 


m 

i 

saf 

j 

■«*!• 

w 

so 

| 

•!»? 

»$■ 

a 


LB..0DA 

EQU0710H;  / 

*  LBYTE  OF  D/A#  0 

W 

7 

HB.ODA 

EQU  071 1H;  / 

*  HBYTK  OF  D/A  M  0 

W 

V 

LB. IDA 

EQU  07 1 2H;  / 

*  LBYTE  OF  D/A  #  1 

w 

V 

MB. IDA 

EQU0713II;  / 

*  HBYTK  OF  D/A  #  1 

w 

V 

;  /* 

GENERAL  CONTROL  */ 

CONTROL  EQU  071 -ill; 

/*  CONTROL  BYTE 

w 

7 

STATUS 

EQU  0714H 

/*  STATUS  BYTE 

R 

V 

;  /* 

ANALOG  /  DIGITAL  CONVERTERS  */ 

CHAN. AD  EQU  0715H; 

/*  A/D  INPUT  CHANNEL  # 

w 

/ 

CONV_AD  EQU  0716H; 

/*  START  A/D  CONVERSION 

w 

V 

LB_AD 

EQU  071511;  / 

*  LOW  BYTE  OF  A/D 

R 

V 

HBAD 

EQU071GH;  / 

*  HIGH  BYTE  OF  A/D 

R 

V 

i  /* 

TIMER  CONTROLS  */ 

TACK 

EQU  0717H 

/*  TIMER  1NT  ACKNOWLEDGE  W 

V 

T_DATA 

EQU  0718H; 

/*  TIMER  DATA 

R/W 

V 

T_CONT 

EQU  0719H; 

/*  TIMER  CONTROL  BYTE 

R/W 

V 

;  r 

PARALLEL  PORT  CONTROLS  */ 

PORTA 

EQU  071  Cl  1 

/*  PORT  A 

R/W  * 

/ 

PORTB 

EQU  07101 1 

/*  PORT  B 

R/W  * 

/ 

PORTC 

EQU  071  El  I 

/*  PORT  C 

R/W  * 

/ 

P_CONT 

EQU  071  FI  I 

/*  PORT  CONTROL  BYTE 

w 

V 

USEFUL  TECMAR  DEFINITIONS  THIS  APPLICATION 


c 


v.,^S'!vx*!5 


SAVES34 

EQU 

Oach 

;T(’ONT:  save  timer  3&  t  counts  in  hold  regs  1&2 

ILCT3 

EQU 

lbh 

;TCONT:  1'  DATA  -  -*  hold  reg  3  &  cycle  bold  reg 

DISAS5 

EQU 

OdOh 

jTCONT:  disable  timer  5 

CLRT5 

EQU 

0e5h 

;TCONT:  clear  -  reset  timer  5  high  (TC  low) 

STEPT5 

EQU 

0f5h 

;TCONT:  increment  timer  5 

UENERAL  DEFINITIONS 


PORTB_lRQ_BIT  EQU  08h  ;PORTB  IRQ3  1th  bit  in  OCWl  of  8259,  check  -  param.c 


-.MAXIMUM  EXT  BUF  SIZE  IS  OFFFO  &  MUST  BE  MULT  OF  4  */ 


BUF^EXTO  EQU  Otffch 

BUF_EXTl  EQU  Offfch 

FJEXTO  EQU  40h 

F _EXT l  EQU  80H 


;MAX  NUMBER  OF  EXTERNAL  0  INTERRUPTS  .=  BUF.EXTO 
;max  number  of  external  1  interrupts  —  BUF.KXTl  /  4 
;PORTC  bit  C6  is  used  for  external  0  input 
-.PORTC  bit  C7  is  used  for  external  1  input 


/  •' 


.*****«««««*«****«*t««***«*«***«*****«******«4*t****«**««*****««**«* 


;  INTERRUPT  CONTROLLER  DEFINITIONS 


EOI  EQU  20h 
1NTA00EQU  20h 
INTA01EQU  2 1  h 
INTBOOEQU  OaOh 
1NTB01EQU  Oalh 
INT  TYPE  EQU 


;end  of  interrupt  command  for  8259 
;lst  8259  interrupt  controller  port  1 
;lst  8259  interrupt  controller  port  2 
;2nd  8259  interrupt  controller  port  1 
;2nd  8259  interrupt  controller  port  2 

70h  jstart  of  interrupt  vectors  for  2nd  8259  (IRQ  8-15  ) 


;  STRUCTURE  DEFINITIONS 


;  IC_STRUC  flag  bit  definition 


FJIESPO 

EQU 

Olh  jwaiting  for  response  0 

F_RESPI 

EQU 

02h  jwaiting  for  response  1 

F_AD  EQU 

10b 

jperform  A/D  conversion 

F_DELAY 

EQU 

40h  ;delay  finished 

F  JSAMP 

EQU 

80h  jsampling  finished 

IC_STRUC 

struc 

;  this  definition  corresponds  to  C  struct  of  same  name 

neban  db 

0 

jnumber  of  channels  to  be  converted 

flags  db 

0 

jcontrol  byte  for  interrupt  routine 

stimO  db 

0 

-.stimulus  lero  to  be  sent  to  parallel  port  out 

sliml  db 

0 

;st i in ii  1  u s  one  to  be  sent  to  parallel  port  out 

s  stimO  dw 

0 

;sa tuple  number  to  turn  stimulus  0  on 

s  stiml  dw 

0 

;samplr  number  to  turn  stimulus  1  on 

s  total  dw 

0 

;total  number  of  samples 

scur  dw 

0 

;current  sample  number  (set  to  0  initially) 

buf.  of!  dw 

0 

r  offset  value  for  next  data  location 

buf  seg  dw 

0 

jbufler  segment  location 

res  ofl  dw 

0 

;respons«*  structure  offset  for  next  data  location 

res_seg  dw 

0 

.response  structure  segment  location 

IC.STRUC 

ends 

;  BYTE  OFFSETS  TO 

MATCH  .  RESPONSE  STRUCTURE  IN  C 

;typedef  struct  { 

STIMEO 

EQI' 

0  ;  unsigned  long  int 

stimeO, 

RT1ME0 

EQI' 

4  ; 

rtimeO, 

RTIMEl 

EQE 

x  ; 

rtimel ; 

RESPO  EQI' 

12 

;  unsigned  char  respO, 

RESI’l  EQI  ’ 

13 

;  respl ;  }  RESP.STRUCT 

MACRO  DEFINITIONS 

ft**************************** 


pushreg  macro 

;push  register  macro 

push 

ax 

push 

di 

push 

ds 

push 

dx 

push 

es 

endm 

popreg  macro 

;pop  register  macro 

pop 

es 

pop 

dx 

pop 

ds 

pop 

di 

pop 

ax 

endm 

macro 

x  off 

niov 

al,  SAVES34 

inov 

dx,  T..CONT 

out 

dx,a! 

niov 

al,  H  CT3 

out 

dx,al 

NA-VAC-'.y 


;get  time  from  timers  3  &  4  al,dx,di,es  used 

;savc  counts  in  hold  reg 

;T  DATA  timerl  hold  reg  &  cycle  hold  reg 


les  di,  dword  ptr  x  of! 

mov  dx,  T  DATA 

ins!> 

insb 

insb 

insb 

endm 


;es:di  -  •  x  ofl':x  .seg 
;dx  -  •  T. DATA  reg 
;cs:di  -  tinier  3  hold  reg,  inc  di 

;es:di  <-  timer  4  hold  reg,  inc  di 
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UNCLASSIFIED 


INCLUDE  FILE  -  PARAM.H 

**************************************************************************** 

includes  all  parameter  and  macro  definitions  for  DAC  software 
*********»*************»************»************«***********»**********»**», 


1***************** *********************************  ************************* 

MACRO  DEFINITIONS 

****************************************************************************  I 

#define  out2b(p,d)  outb(p,d);  outb(p,{d)>>8);  /*  out  2  bytes  ,  low  first  */ 

#define  min(a,b)  (  (a)  >  (b)  ?  (b)  :  (a))  /*  minimum  function  * / 

#define  roax(a,b)  (  (a)  >  (b)  ?  (a)  :  (b))  /*  maximum  function  * / 


/ 


GENERAL  CONSTANT  DEFINITIONS 


#define  MAX_TRIAL  200 
#define  MAXLINE  80 
#  define  HIGH  J'REQ  3000 


#define  BIT0 
# define  BITl 
^define  BIT2 
#define  BIT3 
if  define  BIT4 
#define  BIT5 
#define  BIT6 
if  define  BIT7 


0x01 

0x02 

0x04 

0x08 

0x10 

0x20 

0x40 

0x80 


/*  max  number  trials  -  can  change  */ 

/*  max  number  of  chars/line  -  can  change  */ 

/*  threshold  for  high  freq  line  to  go  high  */ 

/*  bit  selectors  */ 


INTERRUPT  STRUCT  FLAG  DEFINITIONS 


/*  IC_STRUCT  ic.flags  bits  */ 

# define  F_RESP0  B1T0  /*  waiting  for  response  0  */ 

#define  F_RESPl  BITl  /*  wailing  for  response  1  */ 

#define  F_AD  BJT4  /*  perform  a/d  conversion  */ 


#  define  F_DEL  AY  BIT6  /*  delay  finished  */ 

# define  F_SAMP  BIT7  /*  sampling  this  trial  finished  */ 

/********** *************** ********* **************** ************************* 

STRUCTURE  DEFINITIONS 


typedef  struct  { 


/*  INT_STRUCT  (interrupt)  structure  */ 


int  intnum; 

unsigned  Loldseg; 
unsigned  i_oldoff; 
void  (*i_new)(); 

JINT_STRUCT  ; 


/*  interupt  request  line  IRQ  #  * / 

/*  segment  from  old  interupt  vector  * / 
/*  offset  from  old  interupt  vector  */ 

/*  name  of  new  routine  for  interupt  */ 


typedef  struct  {  /*  RESP_STRUCT  response  structure  type  */ 

/*  MUST  MATCH  param.asm  byte  offsets  */ 


unsigned  long  int  stimeO, 
rtimeO, 
rtimel; 

unsigned  char  respO, 
respl; 


/*  time  of  stimulus  1  being  turned  on  * / 
/*  time  of  response  to  stimulus  0  */ 

/*  time  of  response  to  stimulus  1  */ 

/*  response  time  0  =  rtimeO  -  stimes[i]  * / 
/*  response  byte  for  stimulus  0  */ 

/*  response  byte  for  stimulus  1  * / 


}  RESP_STRUCT; 


typedef  struct  { 


/*  IC  (interrupt  control)  structure  type  */ 

/*  must  match  IC  assembler  structure  EXACTLY  */ 


unsigned  char  /*  same  as  define  byte  in  asm  * / 

nchan,  /*  number  of  channels  to  convert  */ 

flags;  /*  control  byte  on  interrupt  entry  */ 

/*  see  above  interrupt  flag  definitions  */ 
unsigned  int  /*  same  as  define  word  in  asm  * / 

stim,  /*  lbyte  =  stimulus  0,  hbyte  =  stimulus  1  * / 

s_stim0,/*  after  samp#  =  s_stim0  -  send  stimulus  0  */ 
s_stiml,/*  after  samp#  =  s_stiml  -  send  stimulus  1  */ 

/*  stimulus  1  =  current  stimulus  |  stim  1  */ 

/*  correct  respO  will  clear  stim  0  * / 

s_total,  /*  total  number  of  sample  points  per  trial  * / 

/*  1st  samp=samp#  0,  stiin  cleared  after  s_total  * / 
/*  s_stim0,  s_stiml,  s_total  must  be  all  diff  */ 
s_cur;  /*  current  sample  number  (initially  0)  */ 


unsigned  char 


*buf; 


/*  (low  word)  offset  of  data  buffer  */ 

/*  (high  word)  segment  of  data  buffer  */ 


RESP_STRUCT 


rcsp; 


/*  (low  word)  offset  of  resp  struct  */ 

/*  (high  word)  segment  of  resp  struct  */ 


}  ICLSTRUCT; 

typedef  struct  {  /*  PAR_STRUCT  (parameter)  structure  type  */ 

unsigned  char 
nchan, 
response, 
external, 
reject, 
int_mask, 

*bufs(2]; 

unsigned  int 

ntrial,  /*  number  of  trial  (stimuli)  to  be  presented  * / 

itrial,  /*  current  trial  * / 

s_per,  /*  sampling  period  in  usee  (FI)  */ 

wt_dur,  /*  warning  tone  duration  in  msec  (F4)  */ 

s_total,  /*  number  of  samples  to  be  taken  */ 

s_stimO,  /*  sample  number  to  turn  stimulus  0  on  * / 

stirns[MAX_TRiAL],  /*  stimulus  -  lbytc  =  stimO,  hbyte=stiml  */ 
s_isi[MAX_TRlAL],  /*  #  of  samples  after  stimO  to  turn  stiml  on  */ 
delays[MAX_TRlAL],  /*  trial  delay  between  sampling  in  msec  (F4)  */ 
wt_fore[MAX_TRlAL],  /*  wt  fore  period  =  ms  after  tone  until  samp  */ 
wt_pers[MAX_TRlAL];  /*  warning  tone  period  in  usee  (1/f)  Fl  */ 
int 

fd_dat,  /*  file  descriptor  for  data  output  */ 

fd_extO,  /*  file  descriptor  for  ext  0  interrupt  times  */ 

fd_extl;  /*  file  descriptor  for  ext  1  interrupt  times  * / 

FILE 

*fp_res,  /*  file  pointer  for  response  file  */ 

*fp_doc;  /*  file  pointer  for  documentation  file  */ 

}PAR_STRUCT  ; 


/*  total  number  of  channels  */ 

/*  #  of  responses  expected,  0,1  or  2  * / 

/*  flag  if  ext  int  on  PORTC  using  PORTBJRQ  * / 
/*  flag  to  reject  previous  trial  * / 

/*  old  8259  interrupt  mask  */ 

/*  data  buffers  for  temporary  storage  */ 


title  RESPONSE  INTERRUPT  HANDLER  -  res.int 


.****************«****•**••*•********•*******•*•*******•******•*•**** 


include  param.asm 

extrn  ic:far 


codeseg  segment  para  public 
assume  cs:codeseg,  ds:seg  ic 
public  res_int_ 


res_int_  proc  far 


;entry  point  for  interrupt  (res_int  in  C) 


checkl: 


pushrcg 

;save  ax,  di,  ds,  dx,  es 

mov 

al,  EOI 

{acknowledge  interrupt  to  8259 

out 

1NTA00,  al 

I 

mov 

ax,  seg  ic 

;load  dataseg 

mov 

ds,  ax 

lcs 

di,dword  ptr  ic.res. 

.off  ;es:di  ->  resp  structure 

mov 

dx,  T  .CONT 

;  dx  ->  timer  control  port 

mov 

al,  SAVES34 

{save  count  in  timer  3  £  4  in  hold  reg 

out 

dx,al 

mov 

al,H_CT3 

;T_DATA  reg  ->  timer  3  hold  reg  £  cycle 

out 

dx,al 

mov 

dx,  PORTA 

;dx  ->  input  port 

in 

al,dx 

;al  =  response  -  reset  8255  IBF  latch 

test 

ic.flags,  F_RESP0 

{determine  if  waiting  for  response  0 

j* 

checkl 

xor 

ic.flags,  F_RESP0 

{clear  waiting  for  respO  flag 

mov 

es:jdij[RESPO],al 

{store  response  0 

add 

di,RTIME0 

;es:di  ->  response  time  0 

cmp 

al,ic  .strimO 

{check  if  right  response 

j* 

clear_stim 

{equal  so  clear  stim  £  save  time 

jmp 

short  save_time 

{save  response  time  only 

test 

ic.flags,  F_RESPl 

{waiting  for  response  1 

j* 

rcs_int_exit 

;not  waiting  for  response 

cmp 

al,es:(di][RESPO] 

;not  valid  if  same  as  response  0 

j* 

res_int_exit 

t 

xor 

ic.flags,  F_RESPl 

{clear  waiting  for  response  1  flag 

mov 

es:(di](RESPl],al 

;s*ve  response 

add 

di.RTIMEl 

;es:di  ->  response  time  1 

;check  if  response  1  =  stimulus  1 
;save  time  only 


emp 

al,  ic.stiml 

jne 

save_time 

clear_stim: 

mov 

ah,al 

not 

ah 

mov 

dx,  PORTB 

in 

al,dx 

and 

al,ah 

out 

dx,al 

ah  =  current  response 

make  clear  mask 

dx  ->  output  port 

at  =  current  stimulus 

clear  response  from  current  stimulus 

CAUTION:  This  *ill  reset  PORTB  OBF  &  INTR 


save_time: 

mov  dx,  T_DATA 

insb 

insb 

insb 

insb 

res_int_exit: 

popreg 

iret 

res_int_endp 


;T_DATA  reg  ->  timer  2  hold  reg 
;store  current  time  from  timer  1 

;store  current  time  from  timer  2 

;restorc  ax,  di,ds,dx,es 


codeseg  ends 


/ 


SUBROUTINE  -  RESETJNT(ip) 


PURPOSE: 

reset  an  interrupt  vector  to  the  value  stored  in  the  interrupt  struct 

INPUT: 

ip  -  INTERRUPT  STRUCTURE  POINTER 

OUTPUT: 

none 

DEPENDENCIES: 

param.h  —  general  parameters 
stdio.h  --  C  standard  input/output 
dos.h  —  C  interface  to  dos  routines 

***********•*•*»***•***•*•*•***•***•****»*»•****»»••*»**•*»»**»**»•*« 

#includc  <stdio.h> 

^include  <dos.h> 

^include  "param.h" 

void  reset_int(ip) 

1NT_STRUCT  *  ip;  /*  interupt  info  structure  */ 

{ 

struct  reg  r;  /*  register  structure  for  DOS  int  */ 

r.r_ds  =  ip->i_oIdseg; 
r.r_dx  =  ip->Loldoff; 

r.r_ax  =  SET1NT  |  ip->intnum;  /*  offset  &  pointer  */ 

intcall  (fcr,  &r,  DOSINT);  /*  redirect  interupt  */ 

return; 

} 


/ 


FUNCTION  SUBROUTINE  SAMPQ 


PURPOSE: 

Perform  all  neccessary  operations  to  start  the  tecmar 
board  sampling  data  for  each  new  trial. 

1)  check  if  last  trial  to  be  rejected 

1)  check  if  all  trials  are  finished 

2)  delay  finished  flag  set  ->  call  erroT  &  terminate 

inter  trial  delay  finished  before  high  level  processing  =  fatal  error 

3)  increment  all  ic  parameter  values 

4)  wait  for  delay  finished  flag  to  be  set 

5)  reset  delay  finished  flag 

6)  start  sampling  in  background 

7)  copy  data  from  previous  trial  to  hard  disk 


INPUT: 

none  -  external  ic,  par 

OUTPUT: 

data  from  previous  buffer  to  output  file 

RETURN  (int): 

1  =  normal  delay  -  sampling  initiated  this  trial  -  stored  previous 
0  =  normal  finish  >  no  sampling  initiated  -  last  two  trials  stored 

DEPENDENCIES: 

tm_int()  -  performs  actual  sampling  in  background 

errorQ  -  error  exit  &  terminate  routine 

ic  -  interrupt  delay  structure  (external) 

par  -  experiment  parameter  structure  (external) 

resps  -  response  structure  Tor  storing  responses  &times  (external) 

teemar.h  -  parameters  for  tecmar  board 

param.h  -  general  parameters 

stdio.h  •  C  standard  input/output 


. "**"***7 

^include  <stdio.h> 

^include  "teemar.h" 

^include  "param.h" 

int  samp() 

{ 

extern  PAR_STRUCT  par; 
extern  IC_STRUCT  ic; 
extern  RESP_STRUCT  respsj); 


/*  experiment  parameter  structure  */ 

/*  interrupt  delay  structure  */ 

/*  storage  for  responses  &  times  */ 


/*  CHECK  REJECT  TRIAL  FLAG  &  IF  FINISHED  -  RETURN  (0)  »/ 


if  (par.reject) 

par.itrial— ; 


/*  reject  -  decrement  trials  * / 


else  if  (par.itrial  ==  par.ntrial){  /*  check  if  done  */ 

write(par.fd_dat,par.bufs|(par.itrial-l)&BlT0),2*par.nchan‘>par.s_total); 
close(par.fd_dat);  /*  close  output  data  file  * / 

return(O);  /*  all  finished  return  0  */ 


I*  CHECK  IF  DELAY  HAS  FINISHED  YET  */ 

if(ic.flags  &  F_DELAY)  error(”startsamp:  delay  finished  early  "); 

/*  INCREMENT  PARAMETERS  IN  INTERRUPT  CONTROL  STRUCTURE  ic  */ 

ic.buf  =  par.bufsfpar.itrial  &  BITO];  /*  next  buffer  * / 

ic.s_stiml  =  par.s_stimO  4-  par.s_isi[par.itrial];  /*  samp  #  for  stiml  * / 

ic.stiin  =  par.stimsjpar.itrial);  /*  stimulus  -word  * / 

ic.s_cur  =  0;  /*  current  sample  #  */ 

ic.resp  —  &resps[par.itrialj;  /*  storage  for  responses  */ 

/*  WAIT  FOR  DELAY  TO  END  &.  THEN  RESET  FLAGS  */ 

while(!(ic.flags  &  F_DELAY));  /*  wait  */ 

ic. flags  =  F.RESPO  !  F_RESPl  |  F^VD  '  (ic.flags  &  *(F_DELAY)); 

/*  clear  delay  flag,  set  a/d  ,  waiting  for  response  flags  */ 

/*  Produce  sync  pulse  on  port  b  (5th  bit)  for  ext  sync  of  scope  */ 
outb(PORTB,  0x10); 

/*  START  SAMPLING  */ 

outb(T_CONT,  T5);  /*  T_DATA  ->  timer  5  mode  reg  */ 

out2b(T_DATA,  TM5us);  /*  set  for  usee  delay  * / 

out2b(T_DATA,  par.s_per);  /*  set  count  down  sampling  freq  divider  */ 
outb(T_CONT,  A_LSX  |  S5);  /*  start  sampling  -  arm  &  load  timer  5  */ 


/*  STORE  DATA  FROM  PREVIOUS  TRIAL  IF  NO  REJECT  FLAG  */ 
if  (par.reject)} 

par.reject  =  0;  /*  reset  reject  flag  */ 

}else{  /*  store  previous  trial  */ 

writc(par.fd.dat,par.bufs((par.itrial-l)&  BIT0j,2*par.nchan*par.s_total); 


/*  RETURN  NORMALLY  */ 
return  (1); 


r-.v  v 
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^**t ***********************************************  **************** 

SUBROUTINE  -  SETJNT(ip) 

******************************************************************* 

PURPOSE: 

redirect  an  interrupt  vector  to  an  interrupt  routine 

INPUT: 

ip  -  INTERRUPT  STRUCTURE  POINTER 

OUTPUT: 

none 

DEPENDENCIES: 

param.h  —  general  parameters 
stdio.h  --  C  standard  input/output 
dos.h  —  C  interface  to  dos  routines 

******************************************************************* j 

^/include  <stdio.h> 

#include  <dos.h> 

#  include  "param.h" 

void  setjnt(ip) 

1NT.STRUCT  ‘ip;  /*  interupt  info  structure  */ 

{ 

struct  reg  r;  /*  register  structure  for  DOS  int  */ 

r.r_ax  =  GETVEC  ]  ip->intnum; 

intcall  (&r,  &r,  DOSINT);  /*  get  old  interupt  vector  */ 

ip->i_oldseg  =  r.r_es;  /*  save  old  int  seg  &  offset  */ 

ip->i_oldoff  =  r.r_bx; 

ptoreg(csreg,r.r_dx,  r.r_ds,  ip->i_new);  /*  convert  procedure  name  to  */ 
r.r_ax  =  SETINT  |  ip->intnum;  /*  offset  &  pointer  */ 

intcall  (&r,  &r,  DOSINT);  /*  redirect  interupt  */ 

return; 


/' 


SUBROUTINE  FUNCTION  -  SETUP() 

****************************************** *************************** 

PURPOSE: 

prompt  user  for  output  file  name  &  parameters 
load  parameter  structure  &  open  output  files 


INPUT: 

prompt  -  fpar:  filename  containing  default  parameters  for  input 
-  fout:  output  file  name  (dir  +  fname  (max  7  char)) 

several  files  will  be  produced  with  various  extensions 
external  -  par:  parameter  structure 

OUTPUT: 

fout0.dat  -  output  binary  data  file 
fout0.doc  -  output  documentation  file 
foutO.par  -  output  parameter  file 
foutO.res  -  output  response  file 
foutO.cxO  -  output  external  0  interrupt  times 
foutO.exl  -  output  external  1  interrupt  times 
external  -  par  (modified) 

DEPENDENCIES: 

getpar.c  •  gets  the  parameters  from  a  parameter  file 

checkpar.c  -  prompts  the  user  for  the  parameter  changes 

writepar.c  -  writes  the  parameter  file  to  disk 

stdio.h 

teemar.h 

parara.h 


/ 


^include  <stdio.h> 
^include  "teemar.h" 
#include  "param.h" 


void  sctupQ 


{ 

extern  PAR_STRUCT  par;  /*  parameter  structure  defined  in  param.h  */ 


char  run, 

stringjMAXLINE], 

fpar|MAXLINE], 

fout[MAXLINE  -7 1; 


/*  run  number  * / 

/*  temp  storage  for  a  string  */ 

/*  parameter  file  name  */ 

/*  parameters  must  be  in  proper  format  */ 
/*  output  file  name  (max  7  char)  */ 

/*  various  extensions  will  be  appended  */ 


a 
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FILE  *fp; 

/*  general  file  pointer 

V 

#define  ECAL 

".cal" 

/*  fout.cal  -  calibration  file 

V 

#  define  ED  AT 

".dat" 

/*  fout.dat  -  binary  data  file 

V 

#define  EDOC 

".doc" 

/*  fout.doc  -  documentation  file 

V 

#define  EPAR 

“.par" 

/*  Tout. par  -  parameter  file 

V 

■#  define  ERES 

".res" 

/*  fout.res  -  response  file 

V 

#define  EEX0 

M.ex0" 

/*  fout.exO  -  external  0  int  file  * 

/ 

#define  EEXl 

".exl" 

/*  fout.exl  -  external  1  int  file  * 

/ 

/*  PROMPT  FOR  FILE  NAME  */ 

printff'Enter  file  name  in  the  following  form  drive:directory\fname  0); 
printf("drive  and  directory  arc  optional  0); 
printf("fname  must  be  <=  7  char  (with  NO  .ext):"); 
while  (getline(fout)  ==  0); 


/*  CHECK  FOR  EXISTENCE  OF  A  CALIBRATION  FILE  */ 

sprintf(string,"%s%s“,fout,ECAL);  /*  make  cal  file  name  (NO  RUN  #)  * / 
if((fp  =  fopen(string,  V))  =  =  NULL) 

fprintf(stderr, "setup:  WARNING  no  cal  file  with  matching  nameO); 

else 

fclose(fp); 

/*  DETERMINE  RUN  NUMBER  -  SET  STRING  ->  OUTPUT  PARAMETER  FILE  •/ 

for(run  =  'O';  run  <=  ’9’  ;  run  —  run  -t  l){  /*  max  of  10  runs  * / 

sprintf(string,"/os%c%s",fout,run,EPAR); 
if((fp  =  fopen(string,  "r"))  ==  NULL)  break; 
fclose(fp); 

} 

if  (run  >  ’9’)  error("setup:  more  than  9  runs  with  this  name"); 

/*  ENTER  DEFAULT  PARAMETER  FILE  OR  PROMPT  FOR  VALUES  */ 

printf("Enter  file  name  for  default  parameters  (or  return):"); 

iffgetline(fpar)  ==  0)  /*  allow  user  to  review  &  change  default  par  */ 

checkparQ; 

else{  /*  get  values  from  default  file  */ 

if((fp  =  fopcn(fpar,  "r")  )  =====  NULL) 

errorf'setup:  unable  to  open  input  parameter  file"); 
gctpar(fp); 

checkpar();  /*  allow  user  to  review  parameters  */ 

fclosc(fp); 

) 


/*  WRITE  OUT  PARAMETER  FILE  */ 
if  ((fp  =  fopen(string,  "w")  )  =====  NULL) 


errorf'sctup:  unable  to  open  output  parameter  file"); 

writepar(fp);  /*  write  out  parameter  file  &  close  */ 

fclose(fp); 

/*  OPEN  OUTPUT  DATA  FILE  (BINARY)  */ 

sprintf(string,"%s%c%s",fout,run,EDAT);  /*  make  data  filename  * / 
if  {(par.fd_dat  =  creat(string,l))  ==  -1) 

error("setup:  unable  to  create  data  file"); 

/*  OPEN  RESPONSE  FILE  IF  NECCESSARY  */ 

if(par.response){ 

sprintf(string,"%s%c%s", foul, run, ERES);  /*  make  response  filename  * / 

if((par.fp_res  =  fopen(string,"w"))  ==  NULL) 
errorf'setup:  unable  to  create  response  file"); 

) 

/*  OPEN  EXTERNAL  FILES  IF  NECCESSARY  */ 
if(par.external  >  0 ){ 

sprintf(string,"%s%c%s",fout,run,EEXO);  /*  make  external  0  filename  */ 

if((par.fd_extO  =  creat(string,l))  =r==  -1) 

crror("setup:  unable  to  create  external  0  file"); 

} 


if(par.external  ~=  2){ 

sprintf(string,"%s%c%s",fout,run,EEXl);  /*  make  external  1  filename  * / 

if((par.fd_extl  =  creat(string,l))  ==  -1) 

error("setup:  unable  to  create  external  1  file"); 

} 

/*  OPEN  DOCUMENTATION  FILE  &  ALLOW  USER  TO  INPUT  DOCUMENTATION  */ 

sprintf(string,"%s%c%s",fout,run,LDOC);  /*  make  documentation  fname  * / 
if((par.fp_doc  =  fopen(string,"w"))  ==  NULL) 

error("setup:  unable  to  create  documentation  file"); 
printff’Onter  documentation  for  this  run  here  (CR  to  end):  0); 
printf("Record  A/D  channels  &  externals  here  also  0); 
whilc(getline(string)  !=  0)fprintf(par.fp_doc,"%s0, string); 
printff'documentation  written  to  %s%c%s  0,  fout,run,EDOC); 

/*  RETURN  */ 


return; 


STORAGE.ASM 


>***«********************************************************************* 

> 

include  param.asm 

bufferO  segmentpara  .  public  'data' 
public  bufO_ 

bufO_  db  Offffh  dup  (?) 

bufferO  ends 

bufferl  segmentpara  public  ’data’ 

public  bufl_ 

bufl_  db  Offffh  dup  (?) 

bufferl  ends 

externalO  segment  para  public  ’data’ 
public  extO,  extO_ 

extO_  label  byte 

extO  db  BUFJEXTO  dup  (0) 

externalO  ends 

externall  segment  para  public  ’data’ 
public  extl ,  extl_ 

cxtl_  label  byte 

extl  db  BUF_£XTl  dup  (0) 

externall  ends 

dataseg  segmentpara  public  ’data’ 

public  ic,  ic_ 
ic_  label  byte 

ic  ICLSTRUC  <> 
dataseg  ends 


;access  this  structure  as  ic  in  C 

{allocate  space  for  id  struc 


/ 


SUBROUTINE  -  STOREQ 


*****************************,************»******************* 

PURPOSE: 

Write  responses,  external  interrupt  times,  documentation  files 
Close  all  files. 


INPUT: 

external  par  —  contains  all  file  pointers 
external  extO,extl  —  buffers  containing  ext  int  times 
external  resps  —  array  of  response  structures 

OUTPUT  FILES: 

fout?.doc  --  documentation  file 
fout?.res  —  response  file  if  selected 
fout?.exO  —  ext  int  line  0  times  if  selected 
foutf.exl  —  ext  int  line  1  times  if  selected 

DEPENDENCIES: 

getline()  -  subroutine  to  get  one  line  at  a  time 
stdio.h 
param.h 

************************  ***********^**f ********  *  ***************************** 


^include  <stdio.h> 

#include  "param.h" 

#  include  "tecmar.h" 

void  store() 

{ 

extern  PAR_STRUCT  par;  /*  external  parameter  structure  */ 

extern  RESP_STRUCT  rcspsf);  /*  external  response  structure  */ 
extern  void  extO(),extl();  /*  external  int  buffers  (NOT  func)  */ 
extern  unsigned  long  int 

*extO_p,*extl_p;  /*  pointers  to  end  of  ext  int  buffer  */ 
long  int  rO,r  1 ;  /*  computed  response  times  */ 

int  i;  /*  temporary  index  */ 

char  string[80];  /*  temporary  string  */ 

/*  DOCUMENTATION  FILE  */ 

printf("Enter  any  additional  documentation  here  (Clt  to  cnd):0); 
while(getline(string)  1=  0)printf(par.fp_doc,"%s0,  string); 
fclose(par.fp_doc); 


/ 

/ 


RESPONSE  FILE  (five  items  per  line)  */ 

respO.  respl,  stimO  ->  respO  time,  sliml  ->  rcspl  time,  stimO  time  */ 


if(par. response)! 

for(i=0;  i<  par.ntrial;  i++){ 

/*  set  resp  time  to  0  if  no  response  */ 
rO  =  (resps(i].rtimeO  =  =  0)?  0  :resps[i].rtimeO  -  resps[i].stimeO; 
rl  =  (resps|i].rtimel  =  =  0)?  0  :resps[i].rtimel  -  respsjiJjstimeO 
-  (long  int)  par.s_isi[i]  *  par.s_per/100; 
fprintf(par.fp_res,"%x  %x  %ld  %\d  %luO, 
respsji).respO,respsji).respl,  rO,  rl,  resps[i]jtimeO); 

} 

fclose(par.fp_rcs); 

} 

I*  EXTERNAL  INTERRUPT  LINE  TIMES  FILES  (binary  -  long  int)  */ 

if(par.external){ 

write(par.fd_extO,  extO, 

sizeof(unsigncd  long  int)  *  (extO_p  -  (unsigned  long  int  *)extO)); 
close(par.fd_extO); 
if(par.external  ==  2){ 
write(par.fd_extl,  cxtl, 

siseof(unsigned  long  int)  *  (extl_p  -  (unsigned  long  int  *)extl)); 
close(par.fd_extl); 

) 

) 


return; 


/' 


SUBROUTINE  -  TEC_DONE() 


PURPOSE: 

deallocate  tecmar  board  &  reset  interrupts 

INPUT: 

none  -  external  par 

OUTPUT: 

none 

DEPENDENCIES: 

par - experiment  parameter  structure  (external) 

reset_int()  —  reset  the  interrupt  vector  to  its  original  value 

error))  --  error  &  terminate  routine 

tecmar. h  -  parameters  for  tecmar  board 

param.h  —  general  parameters 

stdio.h  —  C  standard  input/output 

********•«*********»«•*•**•»******•*»«•***«•***«•»»*»*•*»***•»***»*/ 

#includc  <stdio.h> 

#include  "param.h" 

#include  "teemar.h" 

void  tcc_done() 

{ 

extern  PAR_STRUCT  par;  /*  experiment  parameter  structure  */ 

extern  INT_STRUCT  tm,  res,  ext;  /*  interrupt  vector  structures  */ 

outb(T_CONT,  TIM.CLR);  /*  reset  all  timers  */ 

outb(CONTROL,  D1SAUTO);  /*  disable  tecmar  timer  interrupts  */ 

reset_int(&tm);  /*  reset  timer  interrupt  vector  */ 

if(par.response){ 

reset_int(&res);'  /*  reset  response  interrupt  vector  * / 

outb(P_CONT,PC_RSET|  INTE.AIN);  /*  reset  PORTA  interrupt  enable  */ 

} 

if(par.external){ 

reset_int(&ext);  /*  reset  externa)  interrupt  vector  */ 

outb(P_CONT,PC_.RSET!  INTE_B);  /*  reset  PORTB  interrupt  enable  */ 

1 

outb(!NTA01,  pnr.int.mask);  /*  reset  interrupt  mask  */ 

printff'  SAMPLING  DONE  -  TECMAR  &  INTERRUPTS  RESET  0); 
return; 


»  ■  ,j-i  »'  i’,  »c >•  o-tV »’  *  .  »*  •*.  »*..  ♦*.  »*.,« ,  4'  «■  •  ,.»•  ♦*  ■  > 


/' 


yf 


:: 


SUBROUTINE  -  TEC _INIT() 


PURPOSE: 

initialise  the  tecmar  board  &  start  initial  delay 
Next  step  after  this  call  should  be  the  sampling  loop 

Tecmar  Board  Jumpers  are  assumed  to  be  set  at  the  following 
timer  5  to  generate  interupts  (J7  pin  2  to  pin  3) 
timer  interupt  to  IRQ5  (J9  pin  14  to  pin  5) 
parallel  port  A  interupt  to  IRQ4  (J9  pin  12  to  pin  4) 
parallel  port  B  interupt  to  IRQ3  (J9  pin  13  to  pin  3) 

timer  1  is  used  to  generate  the  warning  tone  gating  for  timer  2 
timer  2  is  used  to  generate  the  warning  tone  frequency 
timers  3&4  are.  used  to  form  a  32  bit  real  time  clock 
timer  5  is  used  to  produce  interrupts  for  sampling  &  delays 

PORT  A  is  used  as  the  input  port  to  read  switches  which  cause  interupt 
PORT  B  is  used  as  the  output  port  to  turn  lights  or  devices  on 
PORT  C  is  used  for  handshaking  on  A&B  and  2  input  lines  which 
generate  an  "external"  interrupt  on  Port  B 


INPUT: 


external  -  ic,  par,  bufO,  bufl 


OUTPUT: 


DEPENDENCIES: 

set_int()  -  set  interrupt  vector  by  redirecting  old  vector 
error()  —  error  &  terminate  routine 

tecmar. h  •  parameters  for  tecmar  board 
param.h  —  general  parameters 
stdio.h  —  C  standard  input/output 
dos.h  —  C  library  for  dos  interface 


^include  <stdio.h> 
^include  <dos.h> 

^include  "param.h" 
^include  "tecmar. h" 

# define  !NIT_DELAY  5000 


/*  program  initial  delay  in  msec  */ 


void  tec_init() 


*  \  k*,y  »vy  v  ■> >  v  v  v 


extern  PAft_STRUCT  par;  /*  parameter  structure  */ 

extern  IC .STRUCT  ic;  /*  interrupt  control  structure  * / 

extern  INT.STRUCT  tm,  res,  ext;  /*  interrupt  vector  structures  */ 

extern  void  bufO(),  bufl();  /*  buffers  address  (NOT  functions)  * / 

extern  void  extO(),  extl();  /*  ext  int  buffer  address  (NOT  tunc)  */ 

extern  unsigned  long  int 

*cxt0.p,*extl_p;  /*  pointers  to  ext  int  buffers  */ 

extern  char  ext.f;  /*  ext  int  Qag  =  #  interrupts  */ 

extern  void  tm_int(),  /*  timer  interrupt  handler  */ 

res.  int(),  /*  response  interrupt  handler  * / 

ext.intQ;  /*  external  interrupt  handler  */ 

unsigned  char  ncw.inask;  /*  new  interrupt  mask  tor  8259  */ 

int  temp; 

/***  SET  PARAMETERS  IN  PARAMETER  &  INTERRUPT  CONTROL  STRUCTURE 

par.itrial  =  0;  /*  first  trial  is  number  0  */ 

ic.nchan  —  par.nchan;  /*  number  of  chanels  */ 

ic.s.stimO  par.s.stimO;  /*  sample  to  turn  stim  0  on  */ 

ic.s.total  =  par.s.total;  /*  number  of  samples  */ 

par.bufs;0|  =  bufO;  /*  set  bufer  pointers  */ 

par.bufsllj  —  bufl; 


/’****  SET  UP  MASTER  MODE  REGISTER  &  TIMERS  &  WARNING  TONE  DURATION 

outb(T_CONT,  TIM.CLR);  /*  master  clear  of  all  timers  */ 

outb(T_CONT,MM  REG);  /*  TJ5ATA  ->  MM  reg*/ 
out2b(TJDATA,MM_SET);  /*  enable  data  pointer  auto  increment  */ 

outb(T_CONT,Tl); 

out2b(T  J)ATA,TMl);  /*  set  Timer  1  mode  reg  for  warning  gating  */ 

outb(T_CONT,HOLDTX  i  f  1);  /*  T.DATA  reg  ->  timer  1  hold  reg  */ 

out2b(T_DATA,  par.wt_dur); 

outb(T_CONT,T2);  /*  set  Timer  2  mode  reg  for  warning  freq  */ 

out2b(T_DATA,TM2); 

outb(T_CONT,T3);  /*  set  timers  3  &  4  to  form  32  bit  .1  msec  */ 

out2b(T_DATA,TM3);  /*  counter  for  all  program  timing  */ 

outb(T_CONT,T4); 
out2b(TJDATA,TM4); 

outb(T_CONT,T5);  /*  set  Timer  5  mode  reg  for  delay  interrupts  * / 

out2b(T_DATA,TM5ms);  /*  initially  set  for  milli  sec  range  */ 

/*  CONFIGURE  TECMAR  BOARD  */ 

outb(CONTROL,DlSAUTO  ]  TIME.INT);  /*  timer  intcrupt  &  no  auto  inc  */ 
outb(P_CONT,  PORT.SET);  /*  set  up  port  control  */ 


/*  REDIRECT  INTERRUPT  VECTORS  &  MAKE  NEW  8259  INTERRUPT  MASK  */ 


par.int.mask  =  inb(lNTAOI);  /*  save  old  interrupt  mask  */ 

/*  new  mask  =  old  +  timer  int  -  IRQO  (sys  clock)  */ 
new.mask  =  ((*(1  <<  TM_IRQ)  &  par.int_mask)  [  BITO); 
tm.intnum  =  8  +  TM_IRQ;  /*  Tech.  Ref.  page  5-5  */ 

tm.Lnew  =  tmjnt;  /*  tm.Lnew  ->  timer  int  handler  */ 

set_int(&tro);  /*  redirect  int  vector  &  store  old  */ 

if(par.response)(  /*  responses  expected  */ 

res.intnum  =  8  +  PORTA_IRQ;  /*  Tech.  Ref.  page  5-5  */ 

res.Lnew  =  res_int;  /*  rcs.i_new  ->  response  int  handler  */ 

set_int(&res);  /*  redirect  int  vector  &  store  old  */ 

new__mask  &=  *(1  <<  PORTA_IRQ);  /*  set  mask  to  allow  PORTA_IRQ  */ 
outb(P_CONT,PC_SET  j  INTE_AIN);  /*  enable  input  PORTA  Tor  interrupt  */ 

} 


if(par.external){  /*  ext  interrupts  using  PORTB  IRQ  */ 

ext.intnum  =  8  +  PORTB_IRQ;  /*  IBM/AT  Tech.  Ref.  page  5-5  */ 

ext.Lnew  =  ext  Jnt;  /*  ext.Lnew  ->  external  int  handler  */ 

set_int(&ext);  /*  redirect  int  vector  &  store  old  * / 

new_mask  &—  '(1  <<  PORTB_IRQ);  /*  set  mask  to  allow  PORTBJRQ  */ 
outb(P_CONT,  PC..SET  ]  INTEJB);  /*  enable  output  PORTB  for  interrupt  * / 
extO_p  =  extO;  /*  set  pointer  to  ext  int  buffer  */ 

extl_p  =  extl; 

ext_f  =  par. external;  /*  ext_flag  =  #  of  external  lines  * / 

) 


/•***  RESET  FLAGS  ***7 


inb(HB^AD); 

outb(PORTB,  0); 
inb(PORTA); 

outb(T_CONT,  CLRTX  j  T5); 
outb(T_ACK,  0); 


/*  reset  a/d  done  teemar  flags  */ 

/*  clear  port  b  * / 

/*  reset  port  a  */ 

/*  set  T5  output  high  (low  TC)  */ 
/*  clear  timer  interrupt  */ 


/*  PROMPT  USER  FOR  START:  START  DELAY  &  WARNINC  TONES  &  ENABLE  8259  */ 

outb(T_CONT,  LOADTX  j  T5);  /*  T_DATA  ->  timer  5  load  reg  */ 

out2b(T_DATA,INIT_DELAY);  /*  set  delay  in  load  reg  * / 

printf("tec_init:  waiting  for  CR  to  start  clock,  timers  &  first  delay:"); 
getchar(); 

out(T_CONT,A_LSX  [  S3  [  S4  [  S5);  /*  load  and  arm  timers  */ 

outb(INTA01,  new_mask);  /*  set  8259  for  new  interrupt  mask  */ 


if  (par.wt.dur  !=  0){  /*  produce  warning  tones  */ 

outb(T_CONT,LOADTX  J  Tl);  /*  T_DATA  ->  timer  1  load  reg  »/ 
temp  =  1NIT_DELAY  -  par.wt.dur  -  par.wt_fore|0)  + 

(par  j_per  *  par  j_itimO)  /  1000  ; 
out2b(TJDATA,  temp); 

outb(T_CONT,LOADTX  [  T2);  /*  T_DATA  ->  timer  2  load  reg  */ 
out2b(T_DATA,  par.wt_pers|par.itrial|); 

outb(T_CONT,  A  LSX  j  Si  [  S2);  /*  start  warning  tone  generator  */ 

} 

/•  SET  PARAMETERS  FOR  INITIAL  ENTRY  INTO  START  SAMPLING  ROUTINE  •/ 

par.reject  =1;  /*  reject  last  trial  since  this  is  first  */ 

par.itrial  =1;  /*  par.itrial  will  be  decremented  on  reject  */ 

return; 

) 


INCLUDE  FILE  -  TECMAR.H 


I  . . . . . 

I 

!  CONTAINS  ALL  NECESSARY  DEFINITIONS  FOR  TECMAR  BOARD 

t 

IMPLEMENTED  AS  AN  I/O  DEVICE  AT  0x0710  -  0x07 IF 

R/W  AT  THE  END  OF  THE  ADDRESS  COMMENTS  REFER  TO  READ  &  WRITE 


/ 


# define  TECMAR  0x0710  /*  TECMAR  STARTING  ADDRESS  */ 

/ . * . 

ADDRESS  OF  ALL  INTERNAL  TECMAR  REGISTERS 


/ 


/*  DIGITAL  /  ANALOG  CONVERTERS  •/ 


#  define 

LB0DA 

0x0710 

/»  LBYTE  OF  D/A  #  0 

Write 

# define 

HB0DA 

0x0711 

/*  KBYTE  OF  D/A  #  0 

W  « 

^define 

LB_1DA 

0x0712 

/*  LBYTE  OF  D/A  #  1 

W  • 

£  define 

hb_ida 

0x0713 

/*  KBYTE  OF  D/A  #  1 

w  • 

/*  GENERAL  CONTROL  */ 


#define  CONTROL  0x0714  /*  CONTROL  BYTE  W  •/ 

#define  STATUS  0x0714  /*  STATUS  BYTE  Read  */ 

I*  ANALOG  /  DIGITAL  CONVERTERS  •/ 


#  define  CHAN^AI)  0x0715 

#  define  CONVJU1  0x0715 
#define  LB_AD  0x0715 

#  define  HB_AD  0x0716 


/•  A/D  INPUT  CHANNEL  #  W  */ 
/*  START  A/D  CONVERSION  W  */ 
/*  LOW  BYTE  OF  A/D  R  •/ 

/*  HIGH  BYTE  OF  A/D  R  */ 


/*  TIMER  CONTROLS  •/ 


#  define  T.ACK  0x0717 
# define  T  DATA  0x0718 

#  define  TCONT  0x0718 
# define  T  STAT  0x0719 


/*  TIMER  INT  ACKNOWLEDGE  W  •/ 
/*  TIMER  DATA  R/W  •/ 

/*  TIMER  CONTROL  BYTE  W  */ 
/*  TIMER  STATUS  BYTE  R  •/ 


/*  PARALLEL  PORT  CONTROLS 


/ 


//define  PORTA  0x07 1C 

//define  PORTB  0x07  ID 

//define  PORTC  0x07  IE 

#  define  P_CONT  0x07 1 F 


/*  PORT  A  R/W  */ 

/*  PORT  B  R/W  •/ 

/*  PORT  C  R/W  */ 

/•  PORT  CONTROL  BYTE  W 


USEFULL  VALUES  FOR  ABOVE  REGISTERS 


/*  TECMAR  CONTROL  -  MAY  BE  |  TOGETHER 


//define  DISAUTO  0x80/*  CONTROL:  disable 
//define  AD  .INT  0x40  /*  CONTROL: 

//define  OVRNJNT  0x20  /*  CONTROL: 

//define  TIMEJNT  0x10  /*  CONTROL: 

//  define  PORT JNT  0x08  /*  CONTROL: 

//define  EXTCONV  0x04  /*  CONTROL: 


auto  increment  option  */ 
enable  interrupt  by  a/d  done  */ 
enable  interrupt  by  a/d  overrun  */ 
enable  interrupt  by  timer  */ 
enable  interrupt  by  parallel  port  */ 
enable  external  start  a/d  conversion  */ 


/*  TIMER  CONTROLS  -  GENERAL 


^define  MM  REG 
//define  ALARMl 
//define  ALARM2 


//define 
//define 
//define 
//define 
//define 
//define 
//define 
//define 
//define 
# define 
//define 


LOADTX 

HOLDTX 

H&CTX 

SETTX0xe8 

CLRTX 

STEPTX 


//define  Si  0x01 

//  define  S2  0x02 

//define  S3  0x04 

//define  S4  0x08 

//define  S5  0x10 

//define  ARMSX 
//define  LOADSX 
//define  A  LSX  0x80 


0x17  /*  T.CONT:  T  DATA  ->  master  mode  register  (MM)  */ 

0x07  /*  T_CONT:  T_DATA  ->  alarm  1  register  */ 

OxOf  /*  T.CONT:  T  DATA  ->  alarm  2  register  */ 

/*  T_CONT:  T^DATA  ->  timer  1  mode  register  */ 

/*  T_CONT:  T_DATA  -  >  timer  2  mode  register  */ 

/*  T_CONT:  T_DATA  ->  timer  3  mode  register  */ 

/  T_CONT:  T_DATA  -~>  timer  4  mode  register  */ 

/*  T_CONT:  T_DATA  -  •  timer  5  mode  register  */ 

0x08  /*  T_CONT—  LOADTX  |  Tx;  T_DATA  ->  x  load  reg  */ 

0x10  /*  T_CONT=  HOLDTX  J  Tx;  T_DATA  ->  x  hold  reg  */ 

0x18  /*  T_CONT=  H&CTX  ;  Tx;  T.DATA  ->  x  hold  reg  &  cycle 

/*  T_CONT=  SETTX  ;  Tx;  set  output  of  timer  x  high  */ 

OxeO  /*  T^CONT=  CLRTX  |  Tx;  set  output  of  timer  x  low  */ 

OxfO  /*  T_CONT=  STEPTX  \  Tx;  increment  (or  dec)  timer  x  */ 

/*  select  timer  1  */ 

/*  select  timer  2  */ 

/*  select  timer  3  */ 

/*  select  timer  4  */ 

/*  select  timer  5  */ 

0x20  /*  T.CONT  =  ARMS  ,Sx  ;  arm  timer  x  */ 

0x40  /•  T_CONT--  LOADSX  \  Sx;  load  timer  x  from  x  load  reg  */ 


/  T  CONT—  A&LSX  j  Sx;  arm  &  load  timer  x  */ 


#define  DISASX  OxcO  /*  T_CONT=  DISASX  J  Sx;  disarm  timer  x  from  counting  ’/ 

# define  SAVESX  OxaO  /*  T_CONT=  SAVESX  ]  Sx;  save  timer  x  count  in  hold  reg  x  */ 

^define  D&SVSX  0x80  /*  T_CONT  —  D&SVSX  j  Sx;  disarm  &  save  count  of  timer  x  */ 

#define  DDPS  0xe8  /*  T_CONT:  disable  data  pointer  sequencing  (bit  MM14)  */ 

#define  EDPS  OxeO  /*  T_CONT:  enable  data  pointer  sequencing  (bit  MM14)  */ 

^define  E16BDB  Oxef  /*  T_CONT:  enable  a  16  bit  data  bus  (bit  MM13)  */ 

#define  E8BDB  0xe7  /*  T_CONT:  enable  a  8  bit  data  bus  (bit  MM13)  */ 

#define  ONFOUT  0xe8  /*  T_CONT:  gate  on  FOUT  (bit  MM12)  */ 

#define  OFFFOUT  Oxee  /*  T_CONT:  gate  off  FOUT  (bit  MM12)  */ 

#define  TIM_CLR  Oxff  /*  T_CONT:  master  reset  to  clear  timer  regs  *  / 


/*  TIMER  CONTROLS  -  THIS  APPLICATION  */ 


/***•  MASTER  MODE  SPECIFICATIONS  *'  **/ 

^define  MM_SET  0x9000  /*  T_DATA:  compare  1,2  disabled;  TOD  disabled  * / 
/*  enable  data  pointer  inc,  FOUT  off,  BCD  */ 

^define  Fl  ((long)  1000  *  1000)  /*  Fl  system  clock  is  1  MHi  */ 

#define  F2  ((long)  1000  *  100)  /*  F2  system  clock  is  100  KHi  »/ 

#define  F3  ((long)  1000  *  10)  /*  F3  system  clock  is  10  KHx  * / 

#de6ne  F4  1000  /*  Fl  system  clock  is  1  KHs  */ 

^define  F5  100  /*  F5  system  dock  is  100  Hz  */ 


/•****  TIMER  l  SPECIFICATIONS  ♦****/ 

/•  USED  TO  GENERATE  GATING  FOR  WARNING  TONE  */ 

r  SET  LOAD  REG  l  to  DELAY  -  WT_FORE  -  WT_DUR  */ 

/*  SET  HOLD  REG  1  TO  WTJ)UR  in  msec  */ 

#define  TM1  OxOe42  /*  T_DATA:  timer  1  -  disable  special  gate,  reload  */ 

/*  from  hold  Sc  load  reg,  count  once,  binary  */ 

/*  count  down,  TC  toggled  (set  low  initialy)  */ 

/*  no  gating,  count  on  rising  edge  F4=  lKHs  */ 

/*•*•*  TIMER  2  SPECIFICATIONS  ****•/ 

/'  OUTPUT  WILL  BE  WARNING  TONE  7 

/*  SET  LOAD  REG  2  TO  1/2  WARNING  TONE  PERIOD  7 

#de6ne  TM2  0x2b22  /*  TJDATA:  tim  2  —  disable  special  gate,  reload  from  */ 

/*  from  load  reg,  count  rep,  binary,  count  */ 

/*  down,  TC  toggled  (for  sq.  wave  out)  */ 

/*  gated  from  TC  of  timer  1,  count  rising  Fl  * / 


/**•*  TIMER  3  SPECIFICATIONS  •••••*/ 

/*  USED  AS  COUNTER  WITH  TIMER  4  TO  KEEP  TIME  */ 

/•  COUNTS  F3  (lOKIIz  =  100  usee  period)  */ 


# define  TM3  0x0d29  /*  T.DATA:  timer  3  =  no  gating  ,  rising  edge  of  F3  */ 

/*  no  special  gate,  reload  from  load,  binary  */ 

/*  count  up  rep,  reload  from  load,  high  TC  */ 

/***•  TIMER  4  SPECIFICATIONS  »••♦**/ 

/•  TRIGGERED  BY  TIMER  3  TO  KEEP  TIME  */ 

/•  COUNTS  TC  OF  TIMER  3  TO  MAKE  32  BIT  COUNTER  */ 

# define  TM4  0x002a  /*  T  .DATA:  timer  4  —  no  gating  ,  rising  edge  of  TC  3  */ 

/*  no  special  gate,  reload  from  load,  binary  */ 

/*  count  up  rep,  reload  from  load,  toggle  TC  * / 

/••♦«  TIMER  5  SPECIFICATIONS  «•****/ 

/•  USED  TO  GENERATE  INTERRUPTS  AFTER  A  FIXED  DELAY 

/•  USES  Fl  (TM5us  lMHs  clock  =  micro  second  counter)  */ 

/•  or  F4  (TM5ms  =  1KH*  clock  =  milli  second  counter)  */ 

/'  load  reg  5  <’-  amount  of  delay  (us  or  ms)  */ 

/*  pin  2  to  pin  3  on  Tecmar  J7  for  timer  5  ->  interrupts  * / 

irdefine  TM5us  0x0b25  /*  T  DATA:  timer  5  —  disable  special  gate,  reload  */ 

/*  from  load  reg,  count  rep,  binary,  down  */ 

/*  low  terminal  count  (TC)  pulse  */ 

/*  no  gating,  rising  edge  of  Fl  (1  MHs)  */ 

^define  TM5msOxOe25  /*  same  as  above  but  uses  F4  (lKH*  clock)  */ 


/*  PORT  CONTROLS  -  GENERAL  */ 


^define  PA.OM  0x80 
^define  PA.lM  OxaO 
/^define  PA.2M  OxcO 


^define  PA.IN  0x10 
#define  PB.lM  0x04 
#  define  PB.IN  0x02 
#define  PCU.IN 
^define  PCLJN 


/*  P.CONT:  port  A  set  for  mode  0  operation  */ 

/*  P.CONT:  port  A  set  for  mode  1  operation  */ 

/*  P.CONT:  port  A  set  for  mode  2  operation  */ 

/*  the  following  6  options  may  be  or’d  to  above  */ 

/*  eg:  P_CONT  <-  PA_1M  ;  PA_IN  |  PB_1M  */ 

/*  port  A  set  for  input  (default  =  output)  */ 

/*  port  B  set  for  mode  1  (default  =  mode  0)  */ 

/*  port  B  set  for  input  (default  =  output)  */ 

0x08  /*•  port  C  upper  nybble  =  input  (default  =  out  */ 

0x01  /*  port  C  lower  nybble  =  input  (default  =  out  */ 


^define  PC_SET  0x01  /*  P.CONT:  set  bit  0  in  port  c  to  1  */ 

^define  PC _RSET  0x00  /*  P.CONT:  reset  bit  0  in  port  c  to  0  */ 

/*  bit  §  is  selected  by  bits  d3,d2,dl  in  P.CONT  */ 

/*  ex:  P.CONT  <-  PC.SET  |  INTE.AIN  to  enable  int  */ 

# define  INTEJMN  0x08  /*  set/rset  PORTC  bit  4  -  input  PORTA  interrupt  */ 

/^define  INTE^AOUT  0x0c  /*  set/rset  PORTC  bit  6  -  output  PORTA  interrupt  */ 

^define  INTE_B  0x04  /*  set/rset  PORTC  bit  2  -  PORTB  interrupts  */ 


/ 


PORT  CONTROLS  -  THIS  APPLICATION  */ 


#dcGne  PORT_SET  PA^lM  \  PA_1N  \  PB_lM  j  PCUJN 
/*  mode  1  for  ports  A  &  B  */ 

/*  A  =  input,  B  =  output,  C  upper  2  bits  =  in  */ 

/*  this  corresponds  to  #5  in  tecmar  8255  manual  p  6  */ 

/*  INTERRUPT  VECTORING  -  THIS  APPLICATION  */ 

# define  TM_IRQ  5  /*  timer  interrupt  on  IRQ5  -  low  memory  =  34  */ 

# define  PORTA_IRQ  4  /*  port  A  interrupt  on  IRQ4  -  low  memory  =  30  */ 

#definc  PORTB_IRQ  3  /*  port  B  interrupt  on  [RQ3  -  low  memory  =  2C  */ 

#define  SYS_CLK  0  /*  system  clock  to  be  masked  off  to  avoid  interrupts  */ 

/*  see  page  5-5  in  IBM/AT  Technical  Reference  */ 

#define  INTA00  0x20  /*  1st  8259  int  controller  port  1  -  command  byte  */ 

#define  INTA010x21  /*  1st  8259  int  controller  port  2  -  int  mask  */ 

^define  EOl  0x20  /*  end  of  interrupt  command  -  send  to  INTA00  */ 


title  TIMER  INTERRUPT  HANDLER  -  tmjnt 

.ttt***************************************************************** 


include  param.asm 

extrn  ic:far 


codescg  segment  para  public  ’code1 
assume  cs:codeseg 
public  tm_int_ 


proc 

far 

;entry  point  for  interrupt  (tm_int  in  C) 

pushreg 

;save  ax,  di,  ds,  dx,  es 

rnov 

al.EOl 

;acknowledge  interrupt  to  8259 

out 

INTAOO,  al 

sti 

interrupts  back  on 

mov 

ax,  seg  ic 

;ds  points  to  data  segment 

mov 

ds,  ax 

test 

ic  .flags,  F_AD 

;check  F^AD  flag  bit  if  sampling 

jnr 

sample 

i 

;no  sampling  -  delay  only 

or 

ic.flags,  FJDELAY  ;sct  delay  finished  flag 

mov 

dx,  T_CONT 

;disarm  timer 

mov 

al,  D1SASS 

out 

dx,al 

mov 

al,  CLRT5 

;set  timer  5  to  clear  =  high  (low  TC) 

out 

dx,al 

mov 

al,  STEPT5 

;step  timer  to  force  immediate  clear 

out 

dx,al 

jmp 

tm_int_exit 

;go  to  exit  routine 

sampling  so  set  other  registers 

les 

di,  dword  ptr  ic 

■buf_off  ;  set  es.’di  ->  buffer 

cld 

;set  direction  for  insb  command 

mov 

ax,  0 

;set  channel  number  to  0 

mov 

dx,  CHAN  .AD 

;dx  points  to  channel  register 

out 

dx,  al 

inc 

dx 

;dx  point  to  start  conv  reg 

out 

dx,al 

;start  conv  of  chan  0 

mov 

dx,  STATUS 

;dx  points  to  status  register 

in 

al,dx 

;get  status  register 

test 

al,  8 OH 

{check  if  a/d  done  bit  set 

j* 

stat 

;not  done  yet 

store: 

and 

al,  OfH 

;al  =  channel  number  converted 

inc 

dx 

;dx  =  lbyte  register  of  a/d 

insb 

;es:di  =  lbyte,  di++ 

inc 

dx 

;dx  =  hbyte  register  of  a/d 

insb 

;es:di  =  hbyte,  di-t-+ 

inc 

al 

;al  =  next  channel  number  to  convert 

cmp 

al,ic.nchan 

1 

je 

cdone 

;all  channels  done  ? 

dec 

dx 

;dx  points  to  channel  number  reg 

out 

dx,al 

inc 

dx 

out 

dx,al 

1 

jmp 

short  adc 

;start  conv  of  next  chan 

;all  channels  have  been  converted 

mov 

ic.buf_off,  di 

;save  offset  of  next  free  byte  in  buffer 

mov 

ax,  ic.s_cur 

;ax  =  previous  sample  number 

inc 

ax 

;increment  &  store  new  sample  number 

mov 

ic.s_cur,  ax 

cmp 

ax,  ic.s_stimO 

je 

sO_on 

;timc  to  turn  stimulus  0  on  ?? 

cmp 

ax,  ic.s.stiml 

je 

sl_on 

;timc  to  turn  stimulus  1  on  ?? 

cmp 

ax,  ic.s_total 

jne 

tm_int_exit 

{sampling  not  done  so  exit  &  wait  for  next  int 

or 

ic .flags,  F_SAMP  ;set  sampling  done  flag 

mov 

dx,  T_CONT 

{disarm  timer 

cdone: 


s_done: 


mov 

al,  D1SAS5 

out 

dx,al 

mov 

al,  0 

jmp 

short  stimulus 

sO_on: 


mov 

di,  ST1ME0 

getime 

ic.res_off[dij 

mov 

al,  ic .stimO 

jmp 

short  stimulus 

sl_on: 


mov 

dx,  PORTB 

in 

al,dx 

or 

al,  ic.stiml 

;clear  stimulus 


;di  <-  offset  of  STIME  ptr  in  resp  stru 
;save  time  stimO  is  turned  on 
;al  =  stimulus  0 


;dx  ->  output  port 

;gct  current  stimulus 

;new  stim  =  old  j  stimulus  1 


■*i#l 

V»> 


VJ' 

i 


stimulus: 

mov  dx,  PORTB 
out  dx,al 

tm_int_exit: 

mov  dx,  T_ACK 
out  dx,al 

popreg 
iret 

tm_int_endp 

codeseg  ends 
end 


;dx  points  to  parallel  port  a 
;output  new  stimulus 

;acknowledge  interrupt  to  tecmar 
;restore  ax,  di,  ds,  dx,  es 


— 

1 

/ . * . . . * . 

1 

SUBROUTINE  -  WRITEPAR(fp) 

.1; 

1 

PURPOSE: 

Write  the  parameter  struct  to  the  output  parameter  file 

'■ 

If 

INPUT: 

fp  -  file  pointer  to  output  parameter  file 

external  -  par:  parameter  structure  defined  in  param.h 

'! 

t 

OUTPUT: 

file  pointed  to  by  fp 

•** 

3 

DEPENDENCIES: 

stdio.li 

1 

param.h 

V 

L 

»***»*<»*»*»**»**»*»*»»»**«*»«********»«*»****««**»**»****•»»»**•**»*  j 

4 

y 

#  include  <stdio.h  • 

^include  "param.h" 

1 

1 

:! 

void  writcpar(fp) 

»» 

FILE  *  fp;  /*  parameter  file  pointer  */ 

•* 

*s 

\ 

t> 

1 

extern  PAR  STRUCT  par;  /*  parameter  structure  */ 

int  i;  /*  temporary  index  */ 

h 

-tl 

? 

«' 

* 

f 

»• 

l 

fprintf(fp,”%d  49d  49d  O.par.nchan, par. response, par.external); 
fprintf(fp,"%d  49d  49d  0,par.ntria!,par.s_per,par.wt_dur); 
fprintf(fp,”%d  49d  0,par.s_total,par.s_stimO); 
for(i=0;  i<  par.ntrial;  i++) 

fprintf(fp,"%x  %x  %d  %d  %d  %d0, 

par.stimsji]  &  OxOOff,(par.stims[i]  >>  8)  Sr.  OxOOff,par.s_isi[i], 

par.dclays{ij,par.wt_pers[i),par.wt_forc|i]); 

return; 

• 

*r 

i. 

/ 

r 

A. 

,« 

V 

% 

1 

; 

1 

'1 

'1 

* 

i 
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SUBROUTINE  -  WRITEPAR(fp) 

PURPOSE: 

Write  the  parameter  struct  to  the  output  parameter  file 


INPUT: 

fp  -  file  pointer  to  output  parameter  file 

external  -  par:  parameter  structure  defined  in  param.h 

OUTPUT: 

file  pointed  to  by  fp 


DEPENDENCIES: 

stdio.h 

param.h 

************************************ ****** *************************** 

#include  <stdio.h> 

#include  "param.h" 

void  writepar(fp) 

FILE  *fp;  /*  parameter  file  pointer  */ 

{ 

extern  PAR_STRUCT  par;  /*  parameter  structure  * / 

int  i;  /*  temporary  index  */ 

fprintf(fp,"%d  49d  49d  O,par.nchan,par.response,par.external); 
fprintf(fp,"%d  49d  49d  0,par.ntrial,par.s_per,par.wt_dur); 
fprintf(fp,"%d  49d  0,par.s_total,par.s_stimO); 
for(i=0;  i<  par.ntrial;  i++) 

fprintf(fp,"%x  %x  %d  %d  %d  %d0, 

par.stims(i)  &  OxOOff,(par.stims(i]  >>  8)  &  OxOOff,par.s_isi|i], 

par.delays[i],par.wt_pers[i],par.wt_fore[i]); 


return; 


} 


